《数字信号处理》·实验一
实验题目:离散时间信号的时域分析 教师:杨鉴
专业:电子信息科学与技术 学号:姓名:
Q1.29 运行程序P 1.5, 以产生所有相关的信号
MATLAB程序代码如下:
%程序 p1.5 通过平均的信号平滑 clf; R=51;
d=0.8*(rand(R,1)-0.5); m=0:R-1;
s=2*m.*(0.9.^m); x=s+d';
subplot(2,1,1);
plot(m,d','r-',m,s,'g--',m,x,'b-.'); xlabel('时间序号n'); ylabel('振幅'); legend('d[n]','s[n]','x[n]');
x1=[0 0 x]; x2=[0 x 0]; x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2);
plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]');
xlabel('时间序号n'); ylabel('振幅');
10d[n]s[n]x[n]5振幅0-5051015202530时间序号n3540455086y[n]s[n]振幅420051015202530时间序号n35404550 Q 1.30 未污染的信号s[n]是什么样的形式?加性噪声d[n] 是什么样的形式? 信号s[n] 的MATLAB程序代码如下:
clf; R=51; m=0:R-1;
s=2*m.*(0.9.^m); plot(m,s,'b-.'); xlabel('时间序号n'); ylabel('振幅');
legend('未污染信号s[n]');
7未污染信号s[n]654振幅3210051015202530时间序号n35404550
信号d[n]的MATLAB程序代码如下: clf; R=51;
d=0.8*(rand(R,1)-0.5); m=0:R-1;
plot(m,d','r-'); xlabel('时间序号n'); ylabel('振幅');
legend('加性噪声d[n]');
0.4加性噪声d[n]0.30.20.1振幅0-0.1-0.2-0.3-0.4051015202530时间序号n35404550
Q 1.31使用语句x=s+d 能产生被噪声污染的信号吗?若不能,为什么?
答:不能,因为s和d均是 矩阵,d与s要想相加必须匹配,而矩阵d本身与矩阵s不匹配,必须把d转置变为d'才能与s进行相加。
Q 1.32 信号x1,x2 和x3与信号x 之间的关系是什么? 答:延时关系
Q 1.33legend 命令的作用是什么?
答:legend 命令使得对MATLAB所画的图形进行注释