[理学]中南大学现代信号处理课程设计报告1.doc

现代信号处理课程设计报告 姓名： 班级： 学号： 指导老师： 赵亚湘、郭丽梅 2012年10月15 一、课程设计目的： 1．所学知识 3．提高综合运用所学知识独立分析和解决问题的能力； 给定模拟信号： 1）选择采样频率Fs = 5000Hz和合适的信号长度，采样得到序列 x1(n)。求并画出x1(n)及其序列傅里叶变换 |X1(ejw)|。 2）选择采样频率Fs = 1000Hz和合适的信号长度，采样得到序列 x2(n)。求并画出x2(n)及其需列傅里叶变换 |X2(ejw)|。 结果展示： 源程序： function q1 figure(); t=-0.005:0.00005:0.005; %模拟信号 xa=exp(-1000*abs(t)); Ts=(1/5000);n=-25:25; %离散时间信号 Fs=5KHz x=exp(-1000*abs(n*Ts)); Ts1=(1/1000);n1=-5:1:5; %离散时间信号 Fs=1KHz x1=exp(-1000*abs(n1*Ts1)); K=500;k1=0:1:K; %离散时间傅里叶变换(Fs=5kHz) w1=pi*k1/K; X=x*exp(-j*n*w1); %离散时间傅里叶变换 X=real(X); w1=[-fliplr(w1),w1(2:501)]; X=[fliplr(X),X(2:501)]; subplot(2,2,1);plot(t*1000,xa,:);ylabel(X_1); title(Discrete Signal (Fs=5000Hz));hold on;stem(n*Ts*1000,x,k*);hold off; subplot(2,2,2);plot(w1/pi,X);ylabel(|X_1(e^j^\omega)|); title(Discrete-time Fourier Transform (Fs=5000Hz)); X1=x1*exp(-j*n1*w1); %离散时间傅里叶变换(Fs=1kHz) X1=real(X1); w1=[-fliplr(w1),w1(2:K+1)]; X1=[fliplr(X1),X1(2:K+1)]; subplot(2,2,3);plot(t*1000,xa,:);ylabel(x_2); title(Discrete Signal (Fs=1000Hz));hold on;stem(n1*Ts1*1000,x1,k*);hold off; subplot(2,2,4);plot(w1/pi,X1);ylabel(|X_2(e^j^\omega)|); title(Discrete-time Fourier Transform (Fs=1000Hz)); 3） 说明|X1(ejw)|与|X2(ejw)|间的区别，为什么？ 答：前者比后者频谱更加精确，因为采样频率越大信号频谱范围越大即分辨率越好。 2. 已知两系统分别用下面差分方程描述： 试分别写出它们的传输函数，并分别打印曲线 A=1;B=[1,1]; [H,w]=freqz(B,A,whole); figure; subplot(2,1,1);plot(w/pi,abs(H)); xlabel(\omega/\pi); ylabel(|H(e^j^\omega)|); title(y1(n)=x(n)+x(n-1) 幅频响应特性); axis([0,2,0,2.5]); A=1;B=[1,-1]; [H,w]=freqz(B,A,whole); subplot(2,1,2);plot(w/pi,abs(H)); xlabel(\omega/\pi); ylabel(|H(e^j^\omega)|); title(y2(n)=x(n)-x(n-1) 幅频响应特性); axis([0,2,0,2.5]); 3. 已知已调信号，其中调制信号，载波为。 1） 选择合适的采样频率及信号长度，使用FFT分析该已调信号的频谱并画出其幅频特性和相频特性曲线图 结果展示： 源程序： function q3 figure(1); fs=90; N=90; %采样频率和采样点数 Tp=N/fs; n=0:N-1; T=1/fs; t=n*T; x=cos(pi*t).*cos(9*