东南大学4系短学期matlab,信号与系统MATLAB数字信号处理东南大学短学期作业一...-CSDN博客

信号与系统MATLAB数字信号处理东南大学短学期作业一

实验一二、1. >> A=[1,2,3;4,5,6;7,8,9] A = 1 2 3 4 5 6 7 8 9 3. >> B=[9,8,7;6,5,4;3,2,1] B = 9 8 7 6 5 4 3 2 1 >> C=[4,5,6;7,8,9;1,2,3] C = 4 5 6 7 8 9 1 2 3 4. >> A A = 1 2 3 4 5 6 7 8 9 >> B B = 9 8 7 6 5 4 3 2 1 >> C C = 4 5 6 7 8 9 1 2 3 三、1. >> A=[1,2,3];B=[4,5,6]; >> C=A+B C = 5 7 9 >> D=A-B D = -3 -3 -3 >> E=A.*B E = 4 10 18 >> F=A./B F = 0.2500 0.4000 0.5000 >> G=A.^B G = 1 32 729 >> stem(A) >> stem(B) >> stem(C) >> stem(D) >> stem(E) >> stem(F) >> stem(G) 2. (1) >> for i=1:10000 t(i)=(i-1)/1000; f(i)=3-exp(-t(i)); end plot(t,f) (2) >> for i=1:10000 t(i)=(i-1)/1000; f(i)=5*exp(-t(i))+3*exp(-2*t(i)); end plot(t,f) (3) >> for i=1:3000 t(i)=(i-1)/1000; f(i)=exp(-t(i))*sin(2*pi*t(i)); end plot(t,f) (4) >> for i=1:3000 t(i)=(i-1)/1000; f(i)=sin(3*t(i))/(3*t(i)); end plot(t,f) (5) >> clear >> for i=1:6 k(i)=i; f(i)=(-2)^(-k(i)); end >> stem(k,f) (6) >> clear >> for i=1:5 k(i)=i-1; f(i)=exp(k(i)); end >> stem(k,f) (7) >> clear >> for i=1:99 k(i)=i; f(i)= k(i); end >> stem(k,f) 四、 1. >> clear A=[2,3,1;1,1,1;3,-1,-1];B=[11;6;-2]; >> x=A\B x = 1.0000 2.0000 3.0000 2. A=[1,1,1;1,-2,1;1,2,3];B=[2;-1;-1]; >> x=A\B x = 3 1 -2 3. A=[1,1,0;0,1,1;1,0,1];B=[27;33;30]; >> x=A\B x = 12 15 18 实验二二、 1. %定义阶跃函数 function [r] = pulse( k ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here if k> clear >> for i=1:10 k(i)=i-4; f(i)=sin(k(i))*pulse(k(i)); end >> stem(k,f) 2. >> for i=1:10 k(i)=i-5; f(i)=sin(k(i))+sin(pi*k(i)); end >> stem(k,f) 3. >> clear for i=1:50 k(i)=i; f(i)=k(i)*sin(k(i))*pulse(k(i)-3); end stem(k,f) 4. clear for i=1:20 k(i)=i-10; f(i)=(-1)^k(i)+(0.5)^k(i)*pulse(k(i)); end >> stem(k,f) 七、 %直接调用conv函数 1. f1=[1,1,1,1];f2=[1,0.5,0.25,0.125,0.0625]; >> conv(f1,f2) ans = 1.0000 1.5000 1.7500 1.8750 0.9375 0.4375 0.1875 0.0625 2. f1=[1,1,1,1];f2=[3,2,1]; >> conv(f1,f2) ans = 3 5 6 6 3 1 %自编卷积函数 function [y] = my_conv( a,b) %UNTITLED Summary of this function goes here % Detailed explanation goes here lengtha=length(a); lengthb=length(b); len=lengtha+lengthb-1; y=linspace(0,0,len); for i=lengtha:-1:1 for j=lengthb:-1:1 y(i+j-1)=a(i)*b(j)+y(i+j-1); end end end 1. f1=[1,1,1,1];f2=[1,0.5,0.25,0.125,0.0625]; >> my_conv(f1,f2) ans = 1.0000 1.5000 1.7500 1.8750 0.9375 0.4375 0.1875 0.0625 2. f1=[1,1,1,1];f2=[3,2,1]; >> my_conv(f1,f2) ans = 3 5 6 6 3 1 实验三只用DFT计算频谱二、已知离散时间序列 x(k)={1,2,3,4,5,6,6,5,4,3,2,1}用DFTF、DFT计算频谱 1.DTFT法计算频谱 x=[1,2,3,4,5,6,6,5,4,3,2,1]; dtft=zeros(100); for i=1:100 w(i)=(i-50)/10; for k=1:12 dtft(i)=dtft(i)+x(k)*exp(-j*(k-1)*w(i)); end end subplot(1,2,1); plot(w,abs(dtft), r- );xlabel( w ); ylabel( DTFT ); title( 幅频特性 ); subplot(1,2,2); plot(w,angle(dtft), r- );xlabel( w ); ylabel( DTFT );title( 相频特性 ); 2.DFT法计算频谱 %法一：直接调用fft函数语句 X=[1,2,3,4,5,6,6,5,4,3,2,1]; F=fftshift(fft(X,1000)); k=[-500:499]*2*pi/1000; subplot(1,2,1); plot(k,abs(F), b- ) xlabel( w );ylabel( DFT );title( 幅频特性 ); subplot(1,2,2); plot(k,angle(F) , b- ) xlabel( w );ylabel( DFT );title( 相频特性 ); %法二：自编程序计算 xk=[1,2,3,4,5,6,6,5,4,3,2,1]; N=12; n=0:N-1; k=0:N-1; Xn=xk*exp(-j*2*pi/N).^(n *k);%计算DFT X=(Xn*exp(j*2*pi/N).^(n *k))/N;%计算IDFT subplot(1,2,2);stem(k,abs(Xn));title( |X(n)| );%画图 axis([-1,N,1.1*min(abs(Xn)), 1.1*max(abs(Xn))]);%加坐标 subplot(1,2,1);stem(n,xk);title( x(k) ); axis([-1,N,1.1*min(xk),1.1*max(xk)]); 序列x(n)及其DFT变换 3.DTFT法与DFT法的比较 x=[1,2,3,4,5,6,6,5,4,3,2,1]; dtft=zeros(100); for i=1:100 w(i)=(i-50)/10; for k=1:12 dtft(i)=dtft(i)+x(k)*exp(-j*(k-1)*w(i)); end end subplot(1,2,1); plot(w,abs(dtft), r- );xlabel( w ); title( 幅频特性 ); hold; subplot(1,2,2); plot(w,angle(dtft), r- );xlabel( w );title( 相频特性 ); hold; F=fftshift(fft(X,1000)); k=[-500:499]*2*pi/1000; subplot(1,2,1); plot(k,abs(F), b- ) hold; subplot(1,2,2); plot(k,angle(F) , b- ) 三、用DFT法，计算下列信号的频谱 1. %定义函数 function [ y ] = tripi( x ) %UNTITLED Summary of this function goes here % Detailed explanation goes here if((-pi0 else y=0 end %----------------------------------------- for i=1:1000 g(i)=ee(2/1000*i-1); w(i)=(i-1)*0.2*pi; end for i=1001:10000 g(i)=0; w(i)=(i-1)*0.2*pi; end G=fft(g)/1000; subplot(1,2,1); plot(w(1:50),abs(G(1:50))); xlabel( w );ylabel( G );title( 幅度频谱 ); subplot(1,2,2); plot(w(1:50),angle(G(1:50))) xlabel( w );ylabel( Fi );title( 相位频谱 );