数字信号处理实验全部程序MATLABWord文件下载.docx
- 文档编号:19904444
- 上传时间:2023-01-12
- 格式:DOCX
- 页数:19
- 大小:166.27KB
数字信号处理实验全部程序MATLABWord文件下载.docx
《数字信号处理实验全部程序MATLABWord文件下载.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验全部程序MATLABWord文件下载.docx(19页珍藏版)》请在冰豆网上搜索。
运行结果:
(2)用实现以下序列。
a)x(n)=0.8n0≤n≤15
0:
150.8.^n;
(,'
('
x(n)=0.8^n'
b)x(n)(0.2+3j)n0≤n≤15
15((0.2+3*j)*n);
x(n)((0.2+3*j)*n)'
a)的时间序列b)的时间序列
c)x(n)=3(0.125π0.2π)+2(0.25π0.1π)0≤n≤15
15;
3*(0.125**0.2*)+2*(0.25**0.1*);
('
x(n)=3*(0.125**0.2*)+2*(0.25**0.1*)'
d)将c)中的x(n)扩展为以16为周期的函数x16(n)(16),绘出四个周期
63;
3*(0.125**(n,16)+0.2*)+2*(0.25**(n,16)+0.1*);
x16(n)'
e)将c)中的x(n)扩展为以10为周期的函数x10(n)(10),绘出四个周期
39;
3*(0.125**(n,10)+0.2*)+2*(0.25**(n,10)+0.1*);
x10(n)'
d)的时间序列e)的时间序列
(3)x(n)=[11,3,5],产生并绘出下列序列的样本。
a)x1(n)=2x
(2)
(1)-2x(n)
3;
[1-135];
x1(x,[0-2])2(x,[01])3=2*x12-2*x;
(x3,'
x1(n)=2x
(2)
(1)-2x(n)'
b)
x1(x,[01])2(x,[02])3(x,[03]);
x4(x,[04])5(x,[05]);
1*x1+2*x2+3*x3+4*x4+5*x5;
x2(n)
(1)+2x
(2)+3x(3)+4x(4)+5x(5)'
a)的时间序列b)的时间序列
(4)绘出时间函数的图形,对x轴、y轴图形上方均须加上适当的标注。
a)x(t)(2πt)0≤t≤10sb)x(t)(100πt)(πt)0≤t≤4s
;
t1=0:
0.001:
102=0:
0.01:
(2**t1)(100**t2).*(*t2);
(2,1,1);
(t1)('
t'
x(t)'
x(t)(2**t)'
(2,1,2);
(t2)('
x(t)(100**t2).*(*t2)'
(5)编写函数(n012)实现u(0)1<
n0<
n2,绘出该函数的图形,起点为n1,终点为n2。
n1('
请输入起点:
n2('
请输入终点'
n0('
请输入阶跃位置'
12;
[0>
=0];
u(0)'
2
请输入终点:
8
请输入阶跃位置:
6
(5)运行结果(6)运行结果
(6)给一定因果系统
求出并绘制H(z)的幅频响应与相频响应。
[1-0.670.9];
[1
(2)1];
[hw]();
20*((h));
()('
时间序列t'
幅频特性'
(h);
相频特性'
(右上图)
常用典型序列
单位采样序列
[](n120)
[n12];
[(0)0];
[](-2,8,2);
();
电信1201'
)
n02;
[-10:
10];
(n);
(1);
1
n(i)0
x(i)=1
();
电信1201采样序列第二种方法'
单位阶跃序列
[(0)>
实数指数
[0:
0.9.^n;
复数指数序列
0.1+0.3*j;
(*n);
(x);
(x);
(2,2,1);
()('
(2,2,2);
()('
(2,2,3);
(2,2,4);
正余弦
3*(0.1**3);
例:
求出下列波形
x1(n)=2x(5)-3x(4)
[](x1122)
((n1)(n2)):
((n1)(n2));
y1(1(n));
y21;
y1(((n>
(n1))&
(n<
(n1))1))1;
y2(((n>
(n2))&
(n2))1))2;
12;
[](0)
0;
;
[-2:
[1:
7,61:
1];
[x1111](,5);
[x1212](4);
[x11](2*x11113*x1212);
(n11);
已知某一系统方程为:
y[n][1]+0.9y[2][n]计算并画出脉冲响应h(n)(-20,100)
(-20:
100);
[1-10.9];
(-20,100,0);
()
时间序号N'
脉冲响应h'
脉冲响应电信1201'
实现两个序列的卷积
[3,11,7,01,4,2];
[2,3,05,2,1];
();
(c);
自己给的数
[0.5,1,1.5,0,0];
[1,1,1];
(c);
试求卷积C(t)1(t)*f2(t),并绘制出f1、f2、及卷积以后的波形。
[]=()
%信号处理的改进卷积程序
(1)
(1);
((x))((h));
(:
);
(x,h);
0.1;
t1=[0:
f11.*(t1>
0);
t2=[-1:
0.1:
2];
f22.*(t2).*(t2>
=0)(t2).*(t2<
[](f1122);
(3,1,1)(t11)('
(3,1,2)(t22)('
(3,1,3);
实验二
()
(0:
11);
(*2*);
'
*k;
.^;
*;
[0,1,2,3];
()'
.^();
*;
分析:
因为x(n)是复指数,它满足周期性,我们将在两个周期中的401个频点上作计算来观
10;
(0.9*(j*3)).^n;
200:
200;
(100)*k;
*((*100)).^(n'
*k);
(X);
检验频移特性
100;
(*2);
100:
(100)*k;
(j**4).*x;
(2,2,1)((X));
([-1,1,0,60]);
幅值电信1201'
以为单位的频率'
绝对值X'
(2,2,2)((X));
相位电信1201'
([-1,11,1]);
(2,2,3)((Y));
幅值电信201'
绝对值Y'
(2,2,4)((Y));
相位电信201'
1;
[10.8];
(0.05**n);
n'
x(n)'
输入序列电信1201'
y(n)'
输出序列电信1201'
实验三
已知信号由15幅值0.5的正弦信号和40幅值2的正弦信号组成,数据采样频率
128;
01;
0.5*(2**15*t)+2*(2**40*t);
(0(y)-1)'
*(y);
(y);
128点电信1201'
已知带有测量噪声信号f1=502=1为均值为零、方差为1的随机信号,采样频率为1000,
0.6;
(2**50*t)(2**120*t);
2*(1(t));
(y,512);
.*(Y)/512;
%求功率
1000*(0:
255)/512;
((1:
256));
对信号进行,对其结果进行,并将的结果和原信号进行比较。
128;
01;
(2**40*t)(2**15*t);
(2,2,1)
电信1201'
(2,2,2)
[0()-1];
(2,2,3)
(2,2,4)
电信201'
用函数和计算同一序列的卷积,比较其计算时间。
5000;
*2-1;
1;
x1=0.5*n;
x2=2*n;
t0;
(x12);
(0)
((x1).*(x2));
实验四
用冲激响应不变法设计低通数字滤波器,通带波纹小于1,阻带在
0.2*;
0.3*;
1150.01128;
[](,'
s'
[](N);
[]2();
[](,1);
(,1)
设计一个高通数字滤波器,通带边界频率为300,阻带边界频率为200,通带波纹小于1,阻带衰减大于20,采样频率为1000。
1000;
300/
(2);
200/
(2);
15;
[]();
设计一个24阶带通滤波器,通带频率为
[0.35,0.65];
24;
1(2*);
(b,1,512)
0.002:
m(1:
201)=1(202:
301)=0;
m(302:
351)=0.5(352:
401)=0(402:
501)=1;
r:
2(64);
[1](b);
f11;
(f1(h))
64电信1201'
2(256);
256电信1201'
针对一个含有5、15和30的混和正弦波信号,设计一个带通滤波器。
1=10;
2=20;
100;
[]([7131723],[010],[0.010.010.01],100);
w1=2*1;
w2=2*2;
%使用窗函数
1(n,[w1w2]);
%使用标准频率响应的
(b,1,512);
%数字滤波器频率响应
t=(0:
s=(2**t*5)(2**t*15)(2**t*30);
=(b,1);
%对信号s进行滤波
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 实验 全部 程序 MATLAB