中華民國萬萬歲

用MATLAB程式用sinc(x) fcn 實現類比信號 x(t)之重建,取樣序列&重建類比信號如圖...

序列 x[n] 以取樣間格 Ts = 1/ fs = 0.0002 sec 取樣類比信號 x(t) ,

range  為 -0.005 < t < 0.005 ,因此以有限個數之序列 x[n] , -25 < n < 25 代表取樣序列 .

 >> Ts=0.0002;
>> n=[-25:1:25];
>> nTs=n*Ts;
>> fs=1/Ts;
>> xn=exp(-1000*abs(nTs));
>> t=[-0.005:0.00005:0.005];
>> xr=xn*sinc(fs*(ones(length(n),1)*t-nTs'*ones(1,length(t))));
>> plot(t*1000,xr);
>> plot(t*1000,xr);
>> axis([-5,5,0,1,2]);
??? Error using ==> axis>LocSetLimits
Vector must have 4, 6, or 8 elements.

Error in ==> axis at 96
                LocSetLimits(ax(j),cur_arg);

>> axis([-5,5,0,1.2]);
>> xlabel('t in msec.');
>> ylabel('x[n] and wr(t)');
>> title('time domain analysis');
>> hold on
>> title('time domain analysis');hold
Current plot released
>> title('time domain analysis'); hold
Current plot held
>> title('time domain analysis'); hold on
>> stem (n*Ts*1000,xn);
>> error=max(abs(x-exp(-1000*abs(t))))
??? Undefined function or variable 'x'.

>> error=max(abs(x-exp(-1000*abs(t))));
??? Undefined function or variable 'x'.

>> %check error
>> error=max(abs(x-exp(-1000*abs(t))));
??? Undefined function or variable 'x'.

>> legend('x[n]');
>> gtext('x[n]');
>> gtext('重建信號');
>>