matlab程序(杂)

《matlab程序(杂)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、clearclff1=1000;  %f1=1kHzf2=2500;  %f2=2.5kHzf3=3000;  %f3=3kHzfs=10000;%采样频率f0=10kHZT=1/fs;  %采样周期t=(0:100).*T;y=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);plot(t,y);N=20;    %采样点数Nfft=20;  %FFT点数  不足后面补零W0=fft(y(1:N),Nfft);Aw=2*abs(W0)/N;df=fs/Nfft;  %[0,fs）频点间隔%改变采样点数和f

2、ft变换点数比较分析N1=40;    %采样点数Nfft1=128;  %FFT点数  不足后面补零W1=fft(y(1:N1),Nfft1);Aw1=2*abs(W1)/N1;df1=fs/Nfft1;  %[0,fs）频点间隔figure(2);plot((0:Nfft-1)*df,Aw,'-r');holdonplot((0:Nfft1-1)*df1,Aw1);x(n)=cos(0.48*pi*n)+cos(0.52*pi*n),利用MATLAB程序求如下X（exp(jw)）,X(k)。(1)取x(n)的前10点数据，求N=10点的X（e

3、xp(jw)）,X(k)并作图。(2)取（1）中的x(n)补零至100点，求N=100点的X（exp(jw)）,X(k)并作图。(3)取x(n)的前100点数据，求N=100点的X（exp(jw)）,X(k)并作图。(4)取x(n)的前128点数据，求N=128点的X（exp(jw)）,X(k)并作图。(5)取x(n)的前50点数据，求N=100点的X（exp(jw)）,X(k)并作图。(6)讨论以上五种情况的区别。2．(1)matlab程序如下：figure(1)n=[0:1:9];x=cos(0.48*pi*n)+cos(0.52*pi*n)

4、;w=[0:1:500]*2*pi/500;%0-2*pi区域分为501点subplot(3,1,1);stem(n,x);title('x(n)(0<=n<=9');xlabel('n');axis([0,10,-2.5,2.5]);%axis([xminxmaxyminymax])line([0,10],[0,0]);x1=fft(x);magx1=abs(x1(1:1:10));x=x*exp(-j*n'*w);%内部的矩阵维数必须一致magx=abs(x);k1=0:1:9;w1=2*pi/10*k1;subplot(3,1,2);plo

5、t(w/pi,magx);title('DTFT');xlabel('w');axis([0,1,0,10]);subplot(3,1,3);stem(w1/pi,magx1);title('X(K)');xlabel('frequencyinpiunits');axis([0,1,0,10]);图（1）（2）程序如下：n=[0:1:9];y=cos(0.48*pi*n)+cos(0.52*pi*n);n1=[0:1:99];x=[y(1:1:10),zeros(1,90)];subplot(3,1,1);stem(n1,x);title('x(

6、n)(0<=n<=9+90zeros');xlabel('n');axis([0,100,-2.5,2.5]);%axis([xminxmaxyminymax])line([0,100],[0,0]);w=[0:1:500]*2*pi/500;%0-2*pi区域分为501点x1=fft(x);magx1=abs(x1(1:1:51));x=x*exp(-j*n1'*w);magx=abs(x);k1=[0:1:50];w1=2*pi/100*k1;subplot(3,1,2);plot(w/pi,magx);title('DTFT');xlabe

7、l('w');axis([0,1,0,10]);subplot(3,1,3);stem(w1/pi,magx1);title('X(K)');xlabel('frequencyinpiunits');axis([0,1,0,10]);结果如下：图（2）（3）取x(n)的前100点数据，程序如下：n=[0:1:99];x=cos(0.48*pi*n)+cos(0.52*pi*n);subplot(3,1,1);stem(n,x);title('x(n)(0<=n<=99');xlabel('n');axis([0,100,-2.5,2.5]);%a

8、xis([xminxmaxyminymax])line([0,100],[0,0]);w=[0:1:500]*2*pi/500;%0-2