资源描述：
《MATLAB信号频谱分析资料报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、题1－周期信号频谱分析题目：给定输入信号x(n)为一周期信号，含基波及若干次谐波分量。要求：1、给出信号基波频率及谐波个数2、基波及各个谐波信号的相位和幅值clearclcA=[0.250.50.25];f=[124];fi=[00.5*pi0.75*pi];N=59;fs=10;t=(1:N)./fs;x_n=A(1).*sin(2*pi*f(1)*t+fi(1))+A(2).*sin(2*pi.*f(2)*t+fi(2))+A(3).*sin(2pi.*f(3)*t+fi(3));save('Signal.mat','x_n')程序如下fs=10;N=6
2、0;n=1:N;t=n/fs;x=0.25*sin(2*pi*t)+0.5*sin(2*pi*2*t)+0.25*sin(2*pi*4*t);figure(1)subplot(211);plot(x);grid;y=fft(x,N);mag=abs(y);k=0:length(y)-1;f=fs/N*k;subplot(212);plot(f,mag);xlabel('Frequence(Hz)');ylabel('Magnitude');grid;title('对周期信号的一个周期进行频谱分析')figure(2)x=[xxxxx];figure(2)su
3、bplot(211);plot(x);grid;N=length(x)y=fft(x,N);mag=abs(y);k=0:length(y)-1;f=fs/N*ksubplot(212);plot(f,mag);xlabel('Frequence(Hz)');ylabel('Magnitude');grid;title('对周期信号的多个周期进行频谱分析')题2－数字滤波器设计题目：已知输入信号x(t)=0.6sin(200πt)+sin(400πt)+0.3sin(800πt)。要求:通过MATLAB编程设计一1)低通滤波器,输出信号x(t)的最低频率成份
4、.2)高通滤波器,输出信号x(t)的最高频率成份.3)带通滤波器,输出信号x(t)的中间频率成份.4)带阻滤波器,输出信号x(t)的最低+最高频率成份.(Rp=1dB,Rs=16dB,Fs=3倍最高频率,截止频率、滤波器设计方法及类型自己选择.)解：编写程序如下：Fs=1200;N=1024;n=0:N;t=n/Fs;x=0.6*sin(200*pi*t)+sin(400*pi*t)+0.3*sin(800*pi*t);y=fft(x,N);mag=abs(y);k=0:length(y)-1;f=Fs/N*k;figure(1)plot(f,mag);xl
5、abel('Frequence(Hz)');ylabel('Magnitude');title('N=1024');grid;Wp=0.12*pi;Ws=0.2*pi;rp=1;rs=16;Fs=1200;%数字滤波器指标wp=Wp*Fs;ws=Ws*Fs;%模拟技术指标[n,wn]=buttord(wp,ws,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[b,a]=lp2lp(b,a,wn);figure(2)h1=freqs(b,a);freqs(b,a);title('模拟低通的幅频和相频特性图');
6、[bz,az]=impinvar(b,a,Fs);figure(3)freqz(bz,az,256);title('数字底通的幅频和相频特性图');y1=filter(bz,az,x);y=fft(y1,N);mag=abs(y);k=0:length(y)-1;f=Fs/N*k;figure(10)plot(f,mag);grid;title('输出信号x(t)的最低频率成份')wp=.4*pi;ws=.35*pi;rp=1;rs=16;Fs=1200;%数字滤波器指标wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2);[n,wn]
7、=buttord(wap,was,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[bt,at]=lp2hp(b,a,wap);w=0:pi:4000*pi;figure(4)freqs(bt,at,w);title('模拟高通的幅频和相频特性图');[bz,az]=bilinear(bt,at,Fs);figure(7)freqz(bz,az,256,1);title('数字高通的幅频和相频特性图');y2=filter(bz,az,x);y=fft(y2,N);mag=abs(y);k=0:length(y
8、)-1;f=Fs/N*k;figure(11)plo
