matlab仿真.doc

%题目1clear;clear all;fs=1000;%抽样频率N=256;n=0:N-1;t=n/fs;x=sin(2*pi*10*t)+0.2*sin(2*pi*20*t);figure(1);subplot(211);plot(t,x);xlabel(t);ylabel(y);title(信号时域波形，N=256);grid;y=fft(x,N);%进行fft变换m=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%频率转换figure(1);subplot(212);plot(f,m);axis(0,1000,0,500);xlabel(频率(Hz);ylabel(幅值);title(信号频谱图，N=256);grid;% N=512；clear;clear all;fs=1000;%抽样频率N=512;n=0:N-1;t=n/fs;x=sin(2*pi*10*t)+0.2*sin(2*pi*20*t);figure(1);subplot(211);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(信号时域波形，N=512);grid;%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换m=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%频率转换figure(1);subplot(212);plot(f,m);%频谱图axis(0,1000,0,500);xlabel(频率(Hz);ylabel(幅值);title(信号频谱图，N=512);grid;% N=1024;clear;clear all;fs=1000;%抽样频率N=1024;n=0:N-1;t=n/fs;x=sin(2*pi*10*t)+0.2*sin(2*pi*20*t);figure(1);subplot(211);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(信号时域波形，N=1024);grid;%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换m=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%频率转换figure(1);subplot(212);plot(f,m);%做频谱图axis(0,1000,0,500);xlabel(频率(Hz);ylabel(幅值);title(信号频谱图，N=1024);grid;%题目2：clear;clear all;a=7 10 0 7 1;b=1 2 0 -1;subplot(411);stem(a,fill);title( x(n);grid on;subplot(412);stem(b,fill);title( y(n);grid on;c=conv(a,b)M=length(c)-1;n=0:1:M;subplot(413);stem(n,c,fill);title(线性卷积 x(n)*h(n) );grid on;a2=a zeros(1,4);b2=b zeros(1,5);ak=fft(a2);bk=fft(b2);ck=ak.*bk;c2=i