matlab练习有程序超全

《matlab练习有程序超全》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1000-rOQ-e-e-Q5000G>Q0102030405060理想采样信号序列的相位谱5(!!11-_50102030405060理想采样信号序列1.-1'11111'0102030405060理想采样信号序列的邮度谱5单位冲击信号序列0.5005101520253035404550单位冲击信号的幅度谱50102030405060单位冲击信号序列特定冲击串单位冲击信号序列42I14k•-4ft•<4^••4^••4^•4^-••4^•4b'♦♦4*-•4^…・4^•••♦♦><•■05101520253

2、035404550卩位冲击倍号的幅度谱卷积系统hb[n]输出信号y[n]500输入信号的幅度谱GQJ：0输HI信号的相位谱0501004°0°0101214161416O5oo24oU1ocn400OO程序一fora=1:3kk(a)=4-a;endfora=4:7kk(a)=a-3;endfora=8:15kk(a)=0;endcloseall;subplot(2,1,1);stem(kk);subplot(2,1,2);stem(abs(fft(kk,16)))程序二n=0:15;%定义序列的长度是15a

3、=0.1;f=0.0625;x=exp(-a*n).*sin(2*pi*f*n);closeall;subplot(2,1,1);stem(x);subplot(2,1,2);stem(abs(fft(x)))程序三n=0:15;%定义序列的长度是15P=8；q=2;x=exp(-1*(n-p).A2/q);closeall;subplot(3,1,1);stem(abs(fft(x)))p=8;q=4;x=exp(-1*(n-p).A2/q);subplot(3,1,2);stem(abs(fft(x)))

4、p=8;q=8;x=exp(-1*(n-p).A2/q);subplot(3,1,3);stem(abs(fft(x)))程序四n=1:50;%定义序列的长度是50hb=zeros(1,50);%注意：MATLAB中数组下标从1开始hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;closeall;subplot(3,1,1);stem(hb);title('系统hb[n]*);m=1:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi

5、;T=0.001;%采样率w0=50*sqrt(2.0)*pi;x=A*exp(-a*m*T).*sin(wO*m*T);%pi是MATLAB定义的ti,信号乘可采用"广"subplot(3,1,2);stem(x);title('输入信号x[n]*);y=conv(x,hb);subplot(3,1,3);stem(y);title('输出信号y[n]');k=-25:25;X=x*(exp(-j*pi/12.5)).A(n,*k);magX=abs(X);%绘制x(n)的幅度谱subplot(3,2,1)

6、;stem(magX);title('输入信号的幅度谱');angX=angle(X);%绘制x(n)的相位谱subplot(3,2,2);stem(angX);title(输入信号的相位谱')Hb=hb*(exp(-j*pi/12.5)).A(nf*k);(»归舉gDIP»B@)①壬二(A6UPE3S工9ze)10_dqnsB0舉SUM昴的％HAo_6ueHA6ue-(fflsusDIPesw)皂HpbeluE①迈(gzco)10_dqns«$駆蹙S(UM旋的％工产・u)<.((gCML、一dj)dx①)j

7、ha06二龙06二HU(»归ffigla葺煤咲)①_七二(qgueE①迈(寸zco)10_dqns®归舉gsqJZ昴繃％SH)①一6ueHqH6ue二»超蹙gla區^^)2ssH6eEE①迈(cozco)10_dqns»®ffis(u)qn旋編％sH)sqeHqH6elu(帥亘舉3即恤甬任迫<)①_七1二x6ueE3SHco・Lo)lo_dqns®心ffi*g(u)x軽的％MXr6ueHX6ue-(tin鯨sgnlps+B-e-迫>Q_£(x6eEEaSQ」co)10_dqnsfflgfMgsx昴的％工X)s

8、qeux6euj^*・u<・((g&L、_dJ)dxa)*XHX-(・「_心Ifrnlne+la-s-迫