资源描述：
《现代信号处理课程设计报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、现代信号处理课程设计报告姓名：班级：学号：指导老师：赵亚湘、郭丽梅2012年10月15一、课程设计目的：1．全面复习课程所学理论知识，巩固所学知识重点和难点，将理论与实践很好地结合起来。2.掌握信号分析与处理的基本方法与实现3．提高综合运用所学知识独立分析和解决问题的能力；4．熟练使用一种高级语言进行编程实现。二、课程设计内容1.给定模拟信号：1）选择采样频率Fs=5000Hz和合适的信号长度，采样得到序列x1(n)。求并画出x1(n)及其序列傅里叶变换
2、X1(ejw)
3、。2）选择采样频率Fs=1000Hz和合适的信号长度，采样得到序列x2(n)。求并画出x2(n)及其需列傅里叶变
4、换
5、X2(ejw)
6、。结果展示：源程序：functionq1figure();t=-0.005:0.00005:0.005;%模拟信号xa=exp(-1000*abs(t));Ts=(1/5000);n=-25:25;%离散时间信号Fs=5KHzx=exp(-1000*abs(n*Ts));Ts1=(1/1000);n1=-5:1:5;%离散时间信号Fs=1KHzx1=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
7、=[-fliplr(w1),w1(2:501)];X=[fliplr(X),X(2:501)];subplot(2,2,1);plot(t*1000,xa,':');ylabel('X_1');title('DiscreteSignal(Fs=5000Hz)');holdon;stem(n*Ts*1000,x,'k*');holdoff;subplot(2,2,2);plot(w1/pi,X);ylabel('
8、X_1(e^j^omega)
9、');title('Discrete-timeFourierTransform(Fs=5000Hz)');X1=x1*exp(-j*n1'*
10、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('DiscreteSignal(Fs=1000Hz)');holdon;stem(n1*Ts1*1000,x1,'k*');holdoff;subplot(2,2,4);plot(w1/pi,X1);ylabel('
11、X_2(e^j^omega)
12、');title('Discrete-timeFou
13、rierTransform(Fs=1000Hz)');3）说明
14、X1(ejw)
15、与
16、X2(ejw)
17、间的区别，为什么？答：前者比后者频谱更加精确，因为采样频率越大信号频谱范围越大即分辨率越好。2.已知两系统分别用下面差分方程描述：试分别写出它们的传输函数，并分别打印曲线。说明这两个系统的区别。分析与说明：通过计算得出传输函数，再将得出传输函数的分子分母参数作为条件（A、B），通过程序得出两个函数的幅频响应曲线。结果展示：源程序：functionq2A=1;B=[1,1];[H,w]=freqz(B,A,'whole');figure;subplot(2,1,1);plot(w/pi
18、,abs(H));xlabel('omega/pi');ylabel('
19、H(e^j^omega)
20、');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('
21、H(e^j^omega)
22、');title('y2(n)=x(n)-x(n-1)幅频响应特性');axis([0,2,0,2.5]);3.已知已调信号，其中调制信号，载波
23、为。1）选择合适的采样频率及信号长度，使用FFT分析该已调信号的频谱并画出其幅频特性和相频特性曲线图结果展示：源程序：functionq3figure(1);fs=90;N=90;%采样频率和采样点数Tp=N/fs;n=0:N-1;T=1/fs;t=n*T;x=cos(pi*t).*cos(9*pi*t);%已调信号X=fft(x);mag=abs(X);pha=angle(X);subplot(3,1,1);stem(n,x,'.');title(['采样序列波形图
