数字信号处理 实习大报告 21.docx
- 文档编号:10267424
- 上传时间:2023-05-24
- 格式:DOCX
- 页数:15
- 大小:189.11KB
数字信号处理 实习大报告 21.docx
《数字信号处理 实习大报告 21.docx》由会员分享,可在线阅读,更多相关《数字信号处理 实习大报告 21.docx(15页珍藏版)》请在冰点文库上搜索。
数字信号处理实习大报告21
一、实习任务:
能够利用Matlab熟悉地画图,内容包括:
X,Y轴上的label,每幅图上的title,绘画多条曲线时的legend,对图形进行适当的标注等。
(1)在一幅图上作多幅小图
源程序:
problem1
2)画出一组二维图形
源程序:
problem2
3)画出一组三维图形
源程序:
problem3
(4)画出复数的实部与虚部。
源程序:
problem4
(5)完成对一个源程序进行详细注释。
%----------------------------------------------------------------------------
%Gibbs_Phenomena_CFST.m
%----------------------------------------------------------------------------
clear;(清零)
closeall;
%改变窗函数的长度,观察Gibbs现象
%
T=2;(定义常量T)
Lambda=0.3;
ifLambda>=T/2
Lambda=T/4;(需要小于T/2,若大于则另其为T/4)
end
Len=[13102050200];(数组的定义)
Dimension=401;维数为401
delta_t=T/(Dimension-1);定义函数为0.005
Y=ones(Dimension,1)*2*Lambda/T;定义(401,1)*2*0.3/2结果为均为0.3的401行1列数组
YY=zeros(Dimension,6);定义YY为0矩阵401行6列
figure
(1);画图1
t0=-T/2;-1
p1=2*pi/T*Lambda;p1的计算公式
p2=2/pi;
p3=2*pi/T;
fork=1:
Len(6)循环从1到Len6=200
t=t0-delta_t;
fori=1:
Dimension
t=t+delta_t;
Y(i)=Y(i)+sin(p1*k)/k*p2*cos(p3*k*t);循环CFST变换算出Y0到Y401
end
ifk==Len
(1)
YY(:
1)=Y;如果K为Len
(1)把Y矩阵赋给YY中的第一列
Elseifk==Len
(2)
YY(:
2)=Y;;如果K为Len
(2)把Y矩阵赋给YY中的第2列
Elseifk==Len(3)
YY(:
3)=Y;;如果K为Len(3)把Y矩阵赋给YY中的第3列
Elsefk==Len(4)
YY(:
4)=Y;;如果K为Len(4)把Y矩阵赋给YY中的第4列
Elseifk==Len(5)
YY(:
5)=Y;;如果K为Len(5)把Y矩阵赋给YY中的第5列
end
end
YY(:
6)=Y;
fork=1:
6(K从一取到6递加1)
subplot(3,2,k);画出子图3行2列每个图标号为K
plot([-T/2-Lambda],[00],':
c','LineWidth',2);(C代表颜色,线宽为2)
holdon;保留
plot([-Lambda-Lambda],[01],':
c','LineWidth',2);
holdon;保留
plot([-LambdaLambda],[11],':
c','LineWidth',2);(C代表颜色,线宽为2)
holdon;保留
plot([LambdaLambda],[10],':
c','LineWidth',2);
holdon;保留
plot([LambdaT/2],[00],':
c','LineWidth',2);
holdon;保留
二、计算普通褶积与循环褶积,分别使用时间域与频率域两种方法进行正、反演计算,指出循环褶积计算时所存在的边界效应现象;编写一个做相关分析的源程序。
(1)普通褶积Problem2-1
结果:
xy0=
5111932461714
yx0=
5111932461714
(2)循环褶积源程序:
problem2-2
结果:
C0=
28
39
49
38
(3)在时间域正演计算循环褶积
源程序Problem2-3
结果;
xy=
55
23
14
6
35
22
27
4)在频率域反演计算循环褶积
源程序problem2-4
结果:
C=
28
39
49
38
5)循环褶积的边界效应
源程序:
problem2-5
计算结果:
K<=1时,h=0
K>1时,h≠0
结论:
当N≥N1+N2时,循环褶积等于线性褶积
当N 6)编写一个做相关分析的源程序 源程序: problem2-6 结果 xy= 100 184 109 30 50 135 三、设计一个病态(矩阵)系统,分析其病态程度;找出对应的解决方法(提示: 添加白噪因子)。 (1)设计一个病态矩阵源程序3-1 计算结果: B= 5 7 5 7 (2)解决方案: 源程序: 3-2 d= 18.0000 -8.0000-2.0000i 2.0000 -8.0000+2.0000i a= 5 4 5 4 X= -0.1105+0.0000i 0.3883-0.0000i -0.1105+0.0000i 0.3883-0.0000i 四、设计一个一维滤波处理程序(1、分别做低通、高通、带通、带阻等理想滤波器进行处理;2、窗函数)。 (1)分别做低通、高通、带通、带阻等理想滤波器进行处 结果: 源程序见lvboqi1 结果: 源程序见lvboqi2 1 (2)窗函数;shichuang 此题在是和同学一起商量完成的 六、设计一个二维滤波处理程序(分别做低通、高通等处理)。 (1)原始信号problem5-1 (2)二维低通滤波 程序ditong.m (3)二维高通滤波 程序gaotong.m (4)二维带阻滤波 程序见daizu.m 六、验证时间域的循环褶积对应的是频率域的乘积;线性褶积则不然。 源程序6-1 结果: xy0= 4 13 32 69 87 98 103 70 xy1= 63.7500 62.9910-8.4280i 68.0519+8.5519i 60.2512-0.3112i 59.5000+15.3292i 45.8271-5.6635i 59.8923-0.3923i 55.7364-9.0861i 结论: 比较xy0和xy1结果不相等,验证了时间域的循线性环褶积对应的不是频率域的乘积 xy2= 4 13 32 69 87 98 103 70 xy3= 4.0000 13.0000 32.0000 69.0000 87.0000 98.0000 103.0000 70.0000 结论: 比较xy2和xy3验证时间域的循环褶积对应的是频率域的乘积 七、请用通俗、易懂的语言说明数字信号处理中的一种性质、一条定理或一个算例(顺便利用Matlab对其进行实现)。 证明线性叠加定理: ax(t)+by(t)的谱等于aX(f)+bY(f) 结果
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号处理 实习大报告 21 数字信号 处理 实习 大报