数字信号处理实验报告汇总.docx
- 文档编号:16008337
- 上传时间:2023-07-09
- 格式:DOCX
- 页数:13
- 大小:179.67KB
数字信号处理实验报告汇总.docx
《数字信号处理实验报告汇总.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验报告汇总.docx(13页珍藏版)》请在冰点文库上搜索。
数字信号处理实验报告汇总
数字信号处理实验报告
实验1:
使用MATLAB实现连续和离散时间信号的表示和运算
1.
程序:
t=0:
0.001:
1;
y=sin(10*pi*t)+1./(cos(pi*t)+2)
plot(t,y)
运行结果显示
2.
程序:
A=2;duty=30;
w=pi/3;
n=-12:
12;
y=A*square(w*n,duty);
stem(n,y);
ylabel('x(n)');xlabel('Timeindex')
title('squaresignal');
axis([-12,12,-2.5,2.5]);
运行结果显示为:
3.
(1).
程序:
方法一:
n=-5:
10;
x=2.^(n-1).*(n>=1);
stem(n,x)
方法二:
n=-5:
10;
u=[zeros(1,6),ones(1,10)];%构造序列u
x=2.^(n-1).*u
stem(n,x)
方法三:
functiony=u(n)
y=(n>=0)
将上述程序在当前工作目录保存为u.m,然后在命令窗口输入如下代码:
n=-5:
10;
x=2.^(n-1).*u(n-1);
stem(n,x)
运行结果显示:
(2).
程序:
n=-30:
30;
x=cos(n*pi/10-pi/5);
stem(n,x)
运行结果显示:
4.
程序:
function[y,n]=sigmult(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))&(n<=max(n2))==1))=x2;
y=y1.*y2;
将上述程序另存为sigmult.m,保存在当前工作目录下,然后在命令窗口中输入如下代码:
n1=-3:
3;
f1=[3,-2,5,1,-1,2,-3];
n2=-1:
4;
f2=[9,-2,0,-2,-4,8];
y=sigmult(f1,n1,f2,n2)
[y,n]=sigmult(f1,n1,f2,n2)
运行结果:
数字信号处理实验报告
实验2:
使用MATLAB实现离散信号的描述和运算
1.
程序:
A=3;
omega=pi/4;
phi=pi/3;
n=-20:
20;
x=A*sin(omega*n+phi);
stem(n,x);xlabel('n');ylabel('x(n)');
运行结果:
2.
程序:
A=0.7;
r=5/7;
n=-10:
10;
x=A*r.^n;
stem(n,x);
xlabel(‘n’);ylabel(‘x(n)’);
运行结果:
3.
程序:
n1=-2:
3;
x1=[211200];
n2=-1:
3;
x2=[12120];
x3=conv(x1,x2);
subplot(311);
stem(n1,x1);axis([-34-13]);xlabel(‘n1’);ylabel(‘x1(n)’);
subplot(312);
stem(n2,x2);axis([-34-13]);xlabel(‘n2’);ylabel(‘x2(n)’);
subplot(313)
n3=-3:
6;
stem(n3,x3);axis([-47-110]);xlabel(‘n3’);ylabel(‘x3(n)’);
运行结果:
X3=[2559744000]
4.
程序:
functiony=cirshift(x,m,N)
iflength(x)>N
error('Nmustbe>=thelengthofx')
end
x=[x,zeros(1,N-length(x))];
n=[0:
1:
N-1];
n=mod(n-m,N);
y=x(n+1);
在当前工作目录下另存为cirshift.m,然后在命令窗口中输入如下代码:
n=0:
11;
x=cos(pi/8*n+pi/3);
y=cirshift(x,7,16);
n=0:
15;
x=[x,zeros(1,4)];
subplot(211);
stem(n,x);
title('Originalsequence');
xlabel('n');ylabel('x(n)');axis([-116-11])
subplot(212);
stem(n,y);
title('Circularlyshiftedsequence,N=16');
xlabel('n');ylabel('y(n)');axis([-116-11])
运行结果:
数字信号处理实验报告
实验三.MATLAB实现滤波器的结构变换
1.
信号流图:
2.
信号流图:
数字信号处理实验报告
实验4.MATLAB实现数字滤波器的设计
1.程序:
Wp=40;Ws=240;Rp=3;Rs=40;
[N,fc]=buttord(Wp,Ws,Rp,Rs,'s');
[B,A]=butter(N,fc,'s');
[hf,f]=freqs(B,A,1024);
plot(f,20*log10(abs(hf)));
axis([0,240,-40,5])
xlabel('f/Hz');ylabel('H(dB)');
运行结果:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 实验 报告 汇总