数字信号处理实验全部程序MATLABWord格式.doc
- 文档编号:14624551
- 上传时间:2022-10-23
- 格式:DOC
- 页数:16
- 大小:302.50KB
数字信号处理实验全部程序MATLABWord格式.doc
《数字信号处理实验全部程序MATLABWord格式.doc》由会员分享,可在线阅读,更多相关《数字信号处理实验全部程序MATLABWord格式.doc(16页珍藏版)》请在冰豆网上搜索。
subplot(4,2,1);
stem(n,A,'
fill'
);
xlabel('
时间序列n'
ylabel('
A'
subplot(4,2,2);
stem(n,B,'
时间序列n'
B'
subplot(4,2,3);
stem(n,C,'
A+B'
subplot(4,2,4);
stem(n,D,'
A-B'
subplot(4,2,5);
stem(n,E,'
A.*B'
subplot(4,2,6);
stem(n,F,'
A./B'
subplot(4,2,7);
stem(n,G,'
A.^B'
运行结果:
(2)用MATLAB实现以下序列。
a)x(n)=0.8n0≤n≤15
n=0:
15;
x=0.8.^n;
stem(n,x,'
xlabel('
x(n)=0.8^n'
b)x(n)=e(0.2+3j)n0≤n≤15
x=exp((0.2+3*j)*n);
x(n)=exp((0.2+3*j)*n)'
a)的时间序列 b)的时间序列
c)x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π)0≤n≤15
1:
x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);
xlabel('
x(n)=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)'
d)将c)中的x(n)扩展为以16为周期的函数x16(n)=x(n+16),绘出四个周期
63;
x=3*cos(0.125*pi*rem(n,16)+0.2*pi)+2*sin(0.25*pi*rem(n,16)+0.1*pi);
x16(n)'
e)将c)中的x(n)扩展为以10为周期的函数x10(n)=x(n+10),绘出四个周期
39;
x=3*cos(0.125*pi*rem(n,10)+0.2*pi)+2*sin(0.25*pi*rem(n,10)+0.1*pi);
x10(n)'
d)的时间序列 e)的时间序列
(3)x(n)=[1,-1,3,5],产生并绘出下列序列的样本。
a)x1(n)=2x(n+2)-x(n-1)-2x(n)
3;
x=[1-135];
x1=circshift(x,[0-2]);
x2=circshift(x,[01]);
x3=2*x1-x2-2*x;
stem(x3,'
x1(n)=2x(n+2)-x(n-1)-2x(n)'
b)
x1=circshift(x,[01]);
x2=circshift(x,[02]);
x3=circshift(x,[03]);
x4=circshift(x,[04]);
x5=circshift(x,[05]);
xn=1*x1+2*x2+3*x3+4*x4+5*x5;
stem(xn,'
x2(n)=x(n-1)+2x(n-2)+3x(n-3)+4x(n-4)+5x(n-5)'
a)的时间序列 b)的时间序列
(4)绘出时间函数的图形,对x轴、y轴图形上方均须加上适当的标注。
a)x(t)=sin(2πt)0≤t≤10sb)x(t)=cos(100πt)sin(πt)0≤t≤4s
clc;
t1=0:
0.001:
10;
t2=0:
0.01:
xa=sin(2*pi*t1);
xb=cos(100*pi*t2).*sin(pi*t2);
subplot(2,1,1);
plot(t1,xa);
t'
x(t)'
title('
x(t)=sin(2*pi*t)'
subplot(2,1,2);
plot(t2,xb);
x(t)=cos(100*pi*t2).*sin(pi*t2)'
(5)编写函数stepshift(n0,n1,n2)实现u(n-n0),n1<
n0<
n2,绘出该函数的图形,起点为n1,终点为n2。
n1=input('
请输入起点:
'
n2=input('
请输入终点'
n0=input('
请输入阶跃位置'
n=n1:
n2;
x=[n-n0>
=0];
xlable('
ylable('
u(n-n0)'
2
请输入终点:
8
请输入阶跃位置:
6
(5)运行结果 (6)运行结果
(6)给一定因果系统求出并绘制H(z)的幅频响应与相频响应。
a=[1-0.670.9];
b=[1sqrt
(2)1];
[hw]=freqz(b,a);
fp=20*log(abs(h));
plot(w,fp);
xlabel('
时间序列t'
幅频特性'
xp=angle(h);
plot(w,xp);
相频特性'
(右上图)
常用典型序列
单位采样序列
function[x,n]=impseq(n1,n2,n0)
n=[n1:
n2];
x=[(n-n0)==0];
[x,n]=impseq(-2,8,2);
stem(n,x);
电信1201'
)
n0=-2;
n=[-10:
10];
nc=length(n);
x=zeros(1,nc);
fori=1:
nc
ifn(i)==n0
x(i)=1
end
end
stem(n,x);
电信1201采样序列第二种方法'
单位阶跃序列
function[x,n]=stepseq(n1,n2,n0)
x=[(n-n0)>
[x,n]=stepseq(-2,8,2);
实数指数
n=[0:
x=0.9.^n;
复数指数序列
alpha=-0.1+0.3*j;
x=exp(alpha*n);
real_x=real(x);
image_x=imag(x);
mag_x=abs(x);
phase_x=angle(x);
subplot(2,2,1);
stem(n,real_x);
subplot(2,2,2);
stem(n,image_x);
subplot(2,2,3);
stem(n,mag_x);
subplot(2,2,4);
stem(n,phase_x);
正余弦
x=3*cos(0.1*pi*n+pi/3);
例:
求出下列波形
x1(n)=2x(n-5)-3x(n+4)
function[y,n]=sigadd(x1,n1,x2,n2)
n=min(min(n1),min(n2)):
max(max(n1),max(n2));
y1=zeros(1,length(n));
y2=y1;
y1(find((n>
=min(n1))&
(n<
=max(n1))==1))=x1;
y2(find((n>
=min(n2))&
=max(n2))==1))=x2;
y=y1+y2;
function[y,n]=sigshift(x,m,n0)
n=m+n0;
y=x;
n=[-2:
x=[1:
7,6:
-1:
1];
[x11,n11]=sigshift(x,n,5);
[x12,n12]=sigshift(x,n,-4);
[x1,n1]=sigadd(2*x11,n11,-3*x12,n12);
stem(n1,x1);
已知某一系统方程为:
y[n]-y[n-1]+0.9y[n-2]=x[n]计算并画出脉冲响应h(n),n=(-20,100)
n=(-20:
100);
num=
(1);
den=[1-10.9];
x=impseq(-20,100,0);
h=filter(num,den,x);
stem(n,h)
时间序号N'
脉冲响应h'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 实验 全部 程序 MATLAB
![提示](https://static.bdocx.com/images/bang_tan.gif)