资源描述：
《东南大学短学期matlab实验参考答案》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验二：1（2）%时间t的输入如：0:0.01:5clear;t=0:0.01:5;y=(t-1).*((t>1)-(t>3)).*(t>2);plot(t,y,'r-');3（1）clear;x=[1,-2,4,6,-5,8,10];m=length(x);fori=1:mif(i+2<8)if(i-4>0)x1(i)=3*x(i+2)+x(i-4)-2*x(i);elsex1(i)=3*x(i+2)-2*x(i);endelseif(i-4>0)x1(i)=x(i-4)-2*x(i);elsex1(i)=-2*x(i
2、);endendendn=0:m-1;stem(n,x1);4（3）function[y,n]=sigadd(x1,n1,x2,n2)%实现y(n)=x1(n)+x2(n)%[y,n]=sigadd(x1,n1,x2,n2)%y=在包括n1和n2的n上求和序列%x1长为n1的第一个序列，x2长为n2的第二个序列（n2可与n1不等）if(n1>n2)n=n1;fori=n2+1:n1x2(i)=0;endelsen=n2;if(n2>n1)n=n2;fori=n1+1:n2x1(i)=0;endendendy=x1+x2;
3、实验三：2clear;x=[1,2,3,4,5,6,6,5,4,3,2,1];N=0:11;w=-pi:0.01:pi;m=length(x);n=length(w);fori=1:nF(i)=0;fork=1:mF(i)=F(i)+x(k)*exp(-1j*w(i)*k);endendF=F/10;subplot(2,1,1);plot(w,abs(F),'b-');xlabel('w');ylabel('F');title('幅度频谱');gridsubplot(2,1,2);plot(w,angle(F),'b-'
4、);xlabel('w');X=fftshift(fft(x))/10;subplot(2,1,1);holdon;plot(N*2*pi/12-pi,abs(X),'r.');legend('DIFT算法','DFT算法');subplot(2,1,2);holdon;plot(N*2*pi/12-pi,angle(X),'r.');xlabel('w');ylabel('相位');title('相位频谱');grid3（3）clear;t=-2:0.01:2;N=length(t)-1;f=(t>-1)-2*(t>0
5、)+(t>1);%DFT法fori=1:Nw(i)=(i-1)*0.5*pi;endG=fft(f)/100;subplot(2,1,1);plot(w(1:16),abs(G(1:16)),'r.');xlabel('w');ylabel('F');title('幅度频谱');subplot(2,1,2);plot(w(1:16),angle(G(1:16)),'r.');xlabel('w');ylabel('相位');title('相位频谱');%数值计算法w=0:0.01:25;m=length(w);fori=
6、1:mF(i)=0;fork=1:NF(i)=F(i)+f(k)*exp(-1j*w(i)*t(k));endendF=F/100;subplot(2,1,1);holdon;plot(w,abs(F));xlabel('w');ylabel('F');title('幅度频谱');legend('DFT算法','数值算法');subplot(2,1,2);holdon;plot(w,angle(F));xlabel('w');ylabel('相位');title('相位频谱');实验四：2（3）clear;clcsyms
7、tw;ft=(sin(2*pi*t)/2/pi/t)^2;disp('傅里叶正变换：');F=fourier(ft,w)disp('傅里叶反变换：');f=ifourier(F,t)6（5）clear;clcsymszn;Fz=10*z^2/(z-1)/(z+1);disp('反z变换：');fn=iztrans(Fz)实验六：6（3）clear;b=[1,2,1];a=[1,0,0.61];H=freqz(b,a,128);fori=1:128w(i)=(i-1)/128*pi;endsubplot(2,1,1);pl
8、ot(w,abs(H));xlabel('w');ylabel('幅度');title('幅频曲线');subplot(2,1,2);plot(w,angle(H));xlabel('w');ylabel('rad');title('相频曲线');figure;dbode(b,a,1);6（4）clear;w=0:0
