DSP实验报告――离散信号的产生及运算 精
dsp数字信号处理实验报告华工
![dsp数字信号处理实验报告华工](https://img.taocdn.com/s3/m/3ea77c62eff9aef8941e06f1.png)
实验1 常见离散信号产生与实现一、实验目得1、加深对常用离散信号得理解;2、熟悉使用MATLAB在时域中产生一些基本得离散时间信号。
二、实验原理1、单位抽样序列在MATLAB中可以利用函数实现。
2、单位阶越序列在MATLAB中可以利用函数实现:3、正弦序列在MATLAB中实现过程如下:4、复指数序列在MATLAB中实现过程如下:5、指数序列在MATLAB中实现过程如下:三、预习要求1、预先阅读实验讲义(MATLAB基础介绍);2、讨论正弦序列、复指数序列得性质。
A.绘出信号,当、时、、时得信号实部与虚部图;当时呢?此时信号周期为多少?程序如下:titlez1=-1/12+j*pi/6;titlez2=1/12+j*pi/6;z3=1/12;z4=2+j*pi/6;z5=j*pi/6;n=0:20;x1=exp(titlez1*n);x2=exp(titlez2*n);x3=exp(z3*n);x4=exp(z4*n);x5=exp(z5*n);subplot(5,2,1);stem(n,real(x1));xlabel('n');ylabel('real(x1)'); title('z1=-1/12+j*pi/6时') subplot(5,2,2);stem(n,imag(x1));xlabel('n');ylabel('imag(x1)'); title('z1=-1/12+j*pi/6时') subplot(5,2,3);stem(n,real(x2));xlabel('n');ylabel('real(x2)'); title('z2=1/12+j*pi/6时')subplot(5,2,4);stem(n,imag(x2));xlabel('n');ylabel('image(x2)'); title('z2=1/12+j*pi/6时')subplot(5,2,5);stem(n,real(x3));xlabel('n');ylabel('real(x3)'); title('z3=1/12时')subplot(5,2,6);stem(n,imag(x3));xlabel('n');ylabel('image(x3)'); title('z3=1/12时')subplot(5,2,7);stem(n,real(x4));xlabel('n');ylabel('real(x4)'); title('z4=2+j*pi/6时')subplot(5,2,8);stem(n,imag(x4));xlabel('n');ylabel('image(x4)'); title('z4=2+j*pi/6时')subplot(5,2,9);stem(n,real(x5));xlabel('n');ylabel('real(x5)'); title('z5=j*pi/6时')subplot(5,2,10);stem(n,imag(x5));xlabel('n');ylabel('image(x5)'); title('z5=j*pi/6时')运行结果如下:当Z=pi/6时,序列周期为12。
dsp 研究报告
![dsp 研究报告](https://img.taocdn.com/s3/m/e20b57858ad63186bceb19e8b8f67c1cfbd6ee5e.png)
dsp 研究报告DSP(数字信号处理)研究报告一、引言数字信号处理(DSP)是一种将连续信号转换为离散信号并利用数值计算机技术对其进行处理的领域。
随着计算机和通信技术的迅速发展,DSP在各个领域的应用也变得越来越广泛。
本报告将介绍DSP的基本原理、应用领域以及未来的发展趋势。
二、基本原理DSP的基本原理是对离散信号进行数字化处理。
主要包括信号采样、量化、编码和数值计算等环节。
通过这些处理步骤,可以实现对信号的滤波、变换、压缩和识别等操作。
其中,采样是将连续信号转换为离散信号的过程,量化是测量离散信号幅度的过程,编码是将量化结果转化为二进制数的过程,数值计算是在计算机上对二进制数进行运算和处理的过程。
三、应用领域1. 通信领域:DSP广泛应用于通信系统中的调制解调、编码解码、信道均衡、差错控制等方面。
通过DSP技术,可以实现高效率和高质量的信号传输,提高通信系统的性能。
2. 音频领域:DSP在音频处理方面的应用也非常广泛。
例如,音频信号的降噪、混响、均衡等处理,以及音频压缩、编码、解码等技术都离不开DSP的支持。
3. 映像领域:DSP在映像处理中可以实现图像增强、去噪、边缘检测、图像压缩、图像识别等功能。
这些技术在医学影像、监控系统、数字摄像等方面有重要应用。
4. 传感器信号处理:传感器信号是一些外界环境的模拟信号,通过DSP技术可以对其进行预处理、滤波、增强和识别等操作,获得有用的信息。
5. 音视频编解码:DSP技术在音视频编解码方面有着重要作用。
通过DSP算法,可以将高位率的音视频信号压缩为低位率的信号,实现高效的传输和存储。
四、未来发展趋势随着计算机和通信技术的不断发展,DSP技术也在不断完善和演进。
未来的发展趋势主要包括以下几个方面:1. 高性能和低功耗:随着芯片制造工艺的进步,DSP芯片将实现更高的性能和更低的功耗。
这将推动DSP技术在各个领域的应用向更广泛、更深入的方向发展。
2. 多核并行计算:为了满足大规模信号处理的需求,DSP芯片将趋向于多核并行计算的方向。
DSP实验报告――离散信号的产生及运算 精
![DSP实验报告――离散信号的产生及运算 精](https://img.taocdn.com/s3/m/145bf90f2e3f5727a5e962b4.png)
离散信号的产生及运算报告一、实验目的:1、复习和巩固数字信号处理中离散信号的产生和运算2、学习和掌握用MATLAB产生离散信号的方法3、学习和掌握用MATLAB对离散信号进行运算二、实验原理:1.用MATLAB函数产生离散信号信号是数字信号处理的最基本内容。
没有信号,数字信号处理就没了工作对象。
MATLAB7.0内部提供了大量的函数,用来产生常用的信号波形。
例如,三角函数(sin,cos),指数函数(exp),锯齿波函数(sawtooth), 随机数函数(rand)等。
1 产生被噪声污染的正弦信号用随机数函数产生污染的正弦信号。
2 产生单位脉冲序列和单位阶跃序列按定义,单位脉冲序列为单位阶跃序列为。
3 矩形脉冲信号:在MATLAB 中用rectpuls 函数来表示,其调用形式为:y=rectpuls(t,width,用以产生一个幅值为1,宽度为width,相对于t=0 点左右对称的矩形波信号,该函数的横坐标范围由向量t 决定,是以t=0 为中心向左右各展开width/2 的范围,width 的默认值为1。
例:以t=2T(即t-2×T=0为对称中心的矩形脉冲信号的MATLAB 源程序如下:(取T=1)t=0:0.001:4;T=1;ft=rectpuls(t-2*T,2*T;plot(t,ft;grid on; axis([0 4 –0.5 1.5];4 周期性矩形波(方波)信号在MATLAB 中用square 函数来表示,其调用形式为:y=square(t,DUTY,用以产生一个周期为2π、幅值为±1 的周期性方波信号,其中的DUTY参数表示占空比,即在信号的一个周期中正值所占的百分比。
例如频率为30Hz 的周期性方波信号的MATLAB 参考程序如下:t=-0.0625:0.0001:0.0625;y=square(2*pi*30*t,75;plot(t,y;axis([-0.0625 0.0625 –1.5 1.5];grid on ;2、 MATLAB中信号的运算乘法和加法:离散信号之间的乘法和加法,是指它的同序号的序列值逐项对应相乘和相加。
DSP实验报告--离散时间信号与系统的时、频域表示-离散傅立叶变换和z变换-数字滤波器的频域分析和实现-数字
![DSP实验报告--离散时间信号与系统的时、频域表示-离散傅立叶变换和z变换-数字滤波器的频域分析和实现-数字](https://img.taocdn.com/s3/m/6f3fff72b307e87101f69650.png)
南京邮电大学实验报告实验名称:离散时间信号与系统的时、频域表示离散傅立叶变换和z变换数字滤波器的频域分析和实现数字滤波器的设计课程名称数字信号处理A(双语) 班级学号B13011025姓名陈志豪开课时间2015/2016学年,第1学期实验名称:离散时间信号与系统的时、频域表示实验目的和任务:熟悉Matlab基本命令,理解和掌握离散时间信号与系统的时、频域表示及简单应用。
在Matlab环境中,按照要求产生序列,对序列进行基本运算;对简单离散时间系统进行仿真,计算线性时不变(LTI)系统的冲激响应和卷积输出;计算和观察序列的离散时间傅立叶变换(DTFT)幅度谱和相位谱。
实验内容:基本序列产生和运算:Q1.1~1.3,Q1.23,Q1.30~1.33离散时间系统仿真:Q2.1~2.3LTI系统:Q2.19,Q2.21,Q2.28DTFT:Q3.1,Q3.2,Q3.4实验过程与结果分析:Q1.1运行程序P1.1,以产生单位样本序列u[n]并显示它。
clf;n = -10:20;u = [zeros(1,10) 1 zeros(1,20)];stem(n,u);xlabel('Time index n');ylabel('Amplitude');title('Unit Sample Sequence');axis([-10 20 0 1.2]);Q1.2 命令clf,axis,title,xlabel和ylabel命令的作用是什么?答:clf命令的作用:清除图形窗口上的图形;axis命令的作用:设置坐标轴的范围和显示方式;title命令的作用:给当前图片命名;xlabel命令的作用:添加x坐标标注;ylabel c命令的作用:添加y坐标标注;Q1.3修改程序P1.1,以产生带有延时11个样本的延迟单位样本序列ud[n]。
运行修改的程序并显示产生的序列。
clf;n = -10:20;u = [zeros(1,21) 1 zeros(1,9)];stem(n,u);xlabel('Time index n');ylabel('Amplitude');title('Unit Sample Sequence');axis([-10 20 0 1.2]);Q1.23修改上述程序,以产生长度为50、频率为0.08、振幅为2.5、相移为90度的一个正弦序列并显示它。
数字信号处理实验报告
![数字信号处理实验报告](https://img.taocdn.com/s3/m/fe47a38f32d4b14e852458fb770bf78a64293a7d.png)
一、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握离散时间信号的基本运算和变换方法。
3. 熟悉数字滤波器的设计和实现。
4. 培养实验操作能力和数据分析能力。
二、实验原理数字信号处理(Digital Signal Processing,DSP)是利用计算机对信号进行采样、量化、处理和分析的一种技术。
本实验主要涉及以下内容:1. 离散时间信号:离散时间信号是指时间上离散的信号,通常用序列表示。
2. 离散时间系统的时域分析:分析离散时间系统的时域特性,如稳定性、因果性、线性等。
3. 离散时间信号的变换:包括离散时间傅里叶变换(DTFT)、离散傅里叶变换(DFT)和快速傅里叶变换(FFT)等。
4. 数字滤波器:设计、实现和分析数字滤波器,如低通、高通、带通、带阻滤波器等。
三、实验内容1. 离散时间信号的时域运算(1)实验目的:掌握离散时间信号的时域运算方法。
(2)实验步骤:a. 使用MATLAB生成两个离散时间信号;b. 进行时域运算,如加、减、乘、除等;c. 绘制运算结果的时域波形图。
2. 离散时间信号的变换(1)实验目的:掌握离散时间信号的变换方法。
(2)实验步骤:a. 使用MATLAB生成一个离散时间信号;b. 进行DTFT、DFT和FFT变换;c. 绘制变换结果的频域波形图。
3. 数字滤波器的设计和实现(1)实验目的:掌握数字滤波器的设计和实现方法。
(2)实验步骤:a. 设计一个低通滤波器,如巴特沃斯滤波器、切比雪夫滤波器等;b. 使用MATLAB实现滤波器;c. 使用MATLAB对滤波器进行时域和频域分析。
4. 数字滤波器的应用(1)实验目的:掌握数字滤波器的应用。
(2)实验步骤:a. 采集一段语音信号;b. 使用数字滤波器对语音信号进行降噪处理;c. 比较降噪前后的语音信号,分析滤波器的效果。
四、实验结果与分析1. 离散时间信号的时域运算实验结果显示,通过MATLAB可以方便地进行离散时间信号的时域运算,并绘制出运算结果的时域波形图。
DSP(数字信号处理)实验报告2
![DSP(数字信号处理)实验报告2](https://img.taocdn.com/s3/m/0b76af3d31126edb6f1a1088.png)
本科学生实验报告学号124090314 姓名何胜金学院物电学院专业、班级12电子实验课程名称数字信号处理(实验)教师及职称杨卫平开课学期第三至第四学年下学期填报时间2015 年 3 月 1 9 日云南师范大学教务处编印2.产生幅度调制信号x[t]=cos(2t)cos(200t),推导其频率特性,确定抽样频率,并会出波形。
程序: clc,clear,close all t=[0:0.01:5];x=cos(2*pi*t).*cos(200*pi*t); plot(t,x);clc,clear,close allt0=0:0.001:0.1;x0=0.5*(cos(202*pi*t0)+cos(198*pi*t0)); plot(t0,x0,'r') hold on fs=202;t=0:1/fs:0.1;x=0.5*(cos(202*pi*t)+cos(198*pi*t)); stem(t,x);3.对连续信号x[t]=cos(4t)进行抽样以得到离散序列,并进行重建。
(1)生成信号x(t),时间为t=0:0.001:4,画出x(t)的波形。
程序clc,clear,close all t0=0:0.001:3; x0=cos(4*pi*t0); plot(t0 ,x0,'r');(2)以faam=10HZ对信号进行抽样,画出在0≤t≤1范围内的抽样序列,x[k],利用抽样内插函数恢复连续时间信号,画出重逢信号的波形。
程序:clc,clear,close all t0=0:0.001:3; x0=cos(4*pi*t0); plot(t0,x0); hold onfs=10;t=0:1/fs:3; x=cos(4*pi*t); stem(t,x);4.若x[k]是对连续信号x(t)=cos(0.5t)以samf=2Hz抽样得到的离散序列,如何通过在抽样点之间内插,恢复原连续时间信号x(t)?程序:clc,clear,close all t=0:0.0001:4; x=cos(0.5*pi*t); plot(t,x); Figure1:clc,clear,close allt=0:0.0001:4; x=cos(0.5*pi*t); subplot(2,1,1); plot(t,x);t0=0:0.5:4;x0=cos(0.5*pi*t0); subplot(2,1,2); stem(t0,x0);5.已知序列x[k]={1,3,2,-5;k=0,1,2,3},分别取N=2,3,4,5对其频谱X(e j)进行抽样,再由频域抽样点恢复时域序列,观察时域序列是否存在混叠,有何规律?k=[0,1,2,3]; x=[1,3,2,-5]; n=100;omega=[0:n-1]*2*pi/n;X0=1+3*exp(-j*omega)+2*exp(-2*j*omega)-5*exp(-3*j*omega); subplot(3,4,1);stem(k,x);title('原序列');subplot(3,4,2);plot(omega./pi,abs(X0));title('序列的频谱 N=100');N=2;omega=[0:N-1]*2*pi/N;X1=1+3*exp(-j*omega)+2*exp(-2*j*omega)-5*exp(-3*j*omega); subplot(3,4,5);stem(omega./pi,abs(X1));title('频域抽样 N=2');rx1=real(ifft(X1)); subplot(3,4,9);stem(rx1);title('时域恢复');N=3;omega=[0:N-1]*2*pi/N;X2=1+3*exp(-j*omega)+2*exp(-2*j*omega)-5*exp(-3*j*omega); subplot(3,4,6);stem(omega./pi,abs(X2));title('频域抽样 N=3');rx2=real(ifft(X2)); subplot(3,4,10);stem(rx2);title('时域恢复');N=4;omega=[0:N-1]*2*pi/N;X3=1+3*exp(-j*omega)+2*exp(-2*j*omega)-5*exp(-3*j*omega); subplot(3,4,7);stem(omega./pi,abs(X3));title('频域抽样 N=4');rx3=real(ifft(X3)); subplot(3,4,11);stem(rx3);title('时域恢复');。
离散信号分析实验报告
![离散信号分析实验报告](https://img.taocdn.com/s3/m/be139918bf23482fb4daa58da0116c175e0e1e5c.png)
离散信号分析实验报告离散信号分析实验报告引言离散信号分析是一门重要的信号处理技术,广泛应用于通信、图像处理、音频处理等领域。
本实验旨在通过实际操作,探索离散信号分析的基本原理和方法,并通过实验结果验证理论知识的正确性。
实验一:离散信号采样与重构在离散信号分析中,采样是将连续时间信号转换为离散时间信号的过程。
首先,我们使用示波器对连续时间信号进行采样,得到一组离散时间信号。
然后,通过重构技术,将离散时间信号恢复为连续时间信号。
实验中,我们选择了一个正弦信号作为输入信号,通过改变采样频率和重构方法,观察信号的失真情况。
实验结果表明,当采样频率低于信号频率的两倍时,会发生混叠现象,导致信号失真。
而当采样频率高于信号频率的两倍时,信号可以被完全恢复。
此外,使用不同的重构方法也会对信号的失真程度产生影响。
通过实验,我们深入理解了采样和重构的原理,并了解到了如何选择合适的采样频率和重构方法。
实验二:离散信号频谱分析频谱分析是离散信号分析的重要内容之一。
在实验中,我们使用FFT算法对离散信号进行频谱分析,并观察信号在频域上的特征。
通过改变输入信号的频率、幅度和相位,我们可以观察到频谱分析结果的变化。
实验结果表明,在频域上,信号的频谱图呈现出明显的峰值,对应着信号的频率成分。
当输入信号为单频信号时,频谱图上只有一个峰值;而当输入信号为复合信号时,频谱图上会有多个峰值。
此外,改变信号的幅度和相位也会对频谱图产生影响。
通过实验,我们进一步理解了离散信号在频域上的特性,为后续的信号处理工作奠定了基础。
实验三:离散信号滤波滤波是离散信号处理中常用的技术之一。
在实验中,我们使用FIR和IIR两种滤波器对输入信号进行滤波,并比较它们的性能差异。
通过观察输出信号的波形和频谱,我们可以评估滤波器的效果。
实验结果表明,FIR滤波器具有线性相位特性,能够实现较好的频率响应;而IIR滤波器则具有较窄的带宽和较快的响应速度。
根据不同的应用需求,我们可以选择合适的滤波器类型。
DSP实验报告一
![DSP实验报告一](https://img.taocdn.com/s3/m/d9f95c243968011ca30091db.png)
A = 1.5;
arg = 2*pi*f*n - phase;
x = A*cos(arg);
clf; % Clear old graph
stem(n,x); % Plot the generated sequence
axis([0 40 -2 2]);
grid;
title('Sinusoidal Sequence');
subplot(2,1,1);
plot(m,d','r-',m,s,'g--',m,x,'b-.');
xlabel('Time index n');ylabel('Amplitude');
legend('d[n] ','s[n] ','x[n] ');
x1 = [0 0 x];x2 = [0 x 0];x3 = [x 0 0];
1.3修改程序P1.1,以产生带有延时11个样本的延迟单位样本序列ud[n]。运行修改的程序并显示产生的序列。
1.6运行程序P1.2,以产生复数值得指数序列。
1.7哪个参数控制该序列的增长或衰减率?哪个参数控制该序列的振幅?
1.17运行程序P1.4,以产生正弦序列并显示它。
1.18该序列的频率是多少?怎样可以改变它?哪个参数控制该序列的相位?哪个参数控制该序列的振幅?该序列的周期是多少?
title('Real part');
subplot(2,1,2);
stem(n,imag(x));
xlabel('Time index n');ylabel('Amplitude');
离散信号分析实验报告
![离散信号分析实验报告](https://img.taocdn.com/s3/m/307dacfcb1717fd5360cba1aa8114431b90d8eb0.png)
一、实验名称离散信号分析实验二、实验目的1. 理解离散信号的基本概念和特点。
2. 掌握离散信号的表示方法,包括时域和频域表示。
3. 熟悉离散信号的基本运算,如加、减、乘、除等。
4. 理解离散系统响应的概念,并学会使用MATLAB进行离散信号与系统分析。
三、实验原理离散信号是指只在离散时刻上有定义的信号,其特点是时域上的不连续性。
离散信号可以通过时域采样和频域变换进行分析。
四、实验内容1. 离散信号的生成与表示使用MATLAB生成以下离散信号:- 单位脉冲序列:δ[n]- 单位阶跃序列:u[n]- 单位斜坡序列:r[n]- 正弦信号:sin(nω0)- 指数信号:e^(αn)并分别绘制这些信号的时域波形图。
2. 离散信号的运算对上述生成的信号进行以下运算:- 加法运算:δ[n] + u[n]- 乘法运算:δ[n] e^(αn)- 移位运算:δ[n - 1]- 反褶运算:δ[-n]绘制运算结果的时域波形图。
3. 离散系统响应假设离散系统由以下差分方程描述:y[n] = x[n] + x[n - 1] - y[n - 1]使用MATLAB编写程序,对输入信号x[n] = δ[n] 进行仿真,并绘制系统响应y[n] 的时域波形图。
4. 离散信号的频域分析对上述生成的信号进行傅里叶变换,得到其频域表示。
绘制信号的频谱图,并分析信号的频率特性。
五、实验步骤1. 使用MATLAB编写程序,生成上述离散信号。
2. 绘制信号的时域波形图。
3. 对信号进行运算,绘制运算结果的时域波形图。
4. 使用MATLAB编写程序,对输入信号进行仿真,并绘制系统响应的时域波形图。
5. 对信号进行傅里叶变换,绘制信号的频谱图。
六、实验结果与分析1. 离散信号的生成与表示通过实验,我们成功生成了上述离散信号,并绘制了它们的时域波形图。
可以看出,这些离散信号在时域上是不连续的。
2. 离散信号的运算通过实验,我们验证了离散信号的基本运算规律,如加法、乘法、移位和反褶等。
离散信号的实验报告
![离散信号的实验报告](https://img.taocdn.com/s3/m/a1cfc30ae418964bcf84b9d528ea81c758f52e24.png)
离散信号的实验报告实验名称:离散信号的生成与观测实验实验目的:1. 理解离散信号的概念和性质;2. 掌握离散信号的生成和观测方法;3. 练习使用计算机进行信号的生成和观测。
实验原理:离散信号是定义在离散时间上的信号。
在离散时间上,信号的取样是有限的,所以可以用数字表示信号的取值。
离散信号可以看作是连续信号在时间维度上的抽样。
离散信号的生成可以通过以下步骤进行:1. 选择采样周期T和采样点数N;2. 定义连续信号的表达式;3. 选择合适的采样频率,定时采样连续信号的值,得到离散信号。
离散信号的观测可以通过以下步骤进行:1. 将离散信号输入计算机,利用编程语言读取数据;2. 对离散信号进行处理和分析,如时域分析、频域分析等;3. 输出观测结果,进行可视化展示。
实验步骤:1. 设置采样周期T为1 ms,采样点数N为1000;2. 选择正弦信号作为连续信号的表达式,频率f为10 Hz,振幅为1;3. 选择合适的采样频率,如采样频率为100 Hz;4. 在计算机中编写程序,采样并观测离散信号。
实验结果:通过实验,我们得到了离散信号的采样值。
利用计算机对采样值进行处理和分析,得到了时域和频域上的信号特性。
在时域上,我们观察到离散信号在时间轴上的抽样点,可以根据抽样点的密度和分布情况分析信号的周期性、稳定性等特性。
在频域上,我们观察到离散信号的频谱图,可以根据频谱上的峰值和分布情况分析信号的频率成分、能量分布等。
实验讨论:离散信号的生成和观测是数字信号处理的基础,也是许多实际应用中常用的方法。
生成离散信号时,采样频率的选择要满足奈奎斯特采样定理,以避免采样失真。
观测离散信号时,可以利用计算机进行各种信号处理和分析,如滤波、谱分析等。
离散信号的特性和应用都十分广泛。
在通信系统中,离散信号的传输和处理是实现音视频编码和数据传输的重要手段。
在信号处理领域,离散信号的分析和处理常用于故障诊断、数据压缩、图像处理等方面。
离散信号的运算实验报告
![离散信号的运算实验报告](https://img.taocdn.com/s3/m/7215cf5ae97101f69e3143323968011ca300f7ff.png)
一、实验目的1. 理解离散信号的基本概念及其运算规则。
2. 掌握MATLAB在离散信号运算中的应用。
3. 通过实验,验证离散信号运算的基本原理和规律。
二、实验原理离散信号是指在一定时间间隔上取值的信号,其特点是时间离散、幅度连续。
在数字信号处理中,离散信号运算主要包括信号的时域运算、频域运算和变换运算。
三、实验设备与软件1. 实验设备:计算机、MATLAB软件2. 实验数据:常用离散信号数据四、实验内容与步骤1. 信号生成(1)利用MATLAB内置函数生成常用离散信号,如单位脉冲序列、单位阶跃序列、正弦序列等。
(2)绘制信号的波形图,观察信号的时域特性。
2. 时域运算(1)信号相加与相减:将两个离散信号进行相加或相减,观察运算结果。
(2)信号移位:将离散信号进行左移或右移,观察运算结果。
(3)信号反转:将离散信号进行反转,观察运算结果。
(4)信号尺度变换:将离散信号进行尺度变换,观察运算结果。
3. 频域运算(1)快速傅里叶变换(FFT):将离散信号进行FFT变换,观察频谱特性。
(2)频域相乘与相加:将两个离散信号的频谱进行相乘或相加,观察运算结果。
4. 变换运算(1)离散余弦变换(DCT):将离散信号进行DCT变换,观察变换结果。
(2)离散正弦变换(DST):将离散信号进行DST变换,观察变换结果。
五、实验结果与分析1. 信号生成(1)通过MATLAB生成单位脉冲序列、单位阶跃序列、正弦序列等,绘制波形图,观察信号特性。
(2)分析不同信号的特点,如单位脉冲序列的冲击特性、正弦序列的周期特性等。
2. 时域运算(1)信号相加与相减:将两个离散信号进行相加或相减,观察运算结果,验证运算规则。
(2)信号移位:将离散信号进行左移或右移,观察运算结果,验证移位规则。
(3)信号反转:将离散信号进行反转,观察运算结果,验证反转规则。
(4)信号尺度变换:将离散信号进行尺度变换,观察运算结果,验证尺度变换规则。
3. 频域运算(1)快速傅里叶变换(FFT):将离散信号进行FFT变换,观察频谱特性,验证FFT原理。
dsp-软件实验报告(精)
![dsp-软件实验报告(精)](https://img.taocdn.com/s3/m/faf716c77f1922791788e80e.png)
数字信号处理MATLAB仿真实验报告学院:电子工程学院班级: 2011211203学号: 2011210876姓名:孙月鹏班内序号: 04一、实验一:数字信号的 FFT 分析、实验内容及要求(1 离散信号的频谱分析:设信号此信号的0.3pi 和 0.302pi两根谱线相距很近,谱线 0.45pi 的幅度很小,请选择合适的序列长度 N 和窗函数,用 DFT 分析其频谱,要求得到清楚的三根谱线。
(2 DTMF 信号频谱分析用计算机声卡采用一段通信系统中电话双音多频(DTMF)拨号数字 0~9的数据,采用快速傅立叶变换(FFT)分析这10个号码DTMF拨号时的频谱。
2、实验结果x(n的时域图与频谱:得到三根清晰的谱线号码9的频谱号码8的频谱号码7的频谱号码6的频谱、实现代码及分析(1第一小题:k=1000; %DFT点数n=[1:1:k]; %对时域信号进行采样x=0.001*cos(0.45*n*pi+sin(0.3*n*pi-cos(0.302*n*pi-pi/4;subplot(2,1,1;stem(n,x,'.'; %用.画出时域图title('时域序列';xlabel('n';ylabel('x(n';xk=fft(x,k; %进行K点DFT变换w=2*pi/k*[0:1:k-1]; %数字角频率subplot(2,1,2;stem(w/pi,abs(xk; %画出频谱图axis([0.2,0.5,0,2]; %设置窗函数的宽度与限幅title('1000点dft';xlabel('数字频率';ylabel('|xk(k|';% 此题关键在于DFT点数N的确定。
经过计算和实验,当N=1000时能满足题目要求,看到3条清晰地谱线(2)第二小题clear;close all;f=[941 1336;697 1209;697 1336;697 1477;770 1209;770 1336;770 1477;852 1209;852 1336;852 1477] %0-9的频率n=1:400;fs=4000; %取样频率为4000hzfprintf('请输入数字(0 to 9:\n'k=input (''f1=f(k+1,1; %因为从0开始计算,+1得输f2=f(k+1,2; %入数字的两个频率N=400;x1=sin(2*pi*f1*n/fs+sin(2*pi*f2*n/fs;%DTMF的输入信号时域xn=[x1,zeros(1,400]; %补零subplot(2,1,1;plot(xn %画出时域图xlabel('n'ylabel('xn'subplot(2,1,2;fn=fs*n/N; %取样点的频率plot(fn,abs(fft(xn(1:400; %400点fft变换,画出频谱图axis([0,4000,0,300]xlabel('f'ylabel('FFT'二、实验二: DTMF 信号的编码1、实验内容及要求1)把您的联系电话号码通过DTMF 编码生成为一个 .wav 文件。
DSP实验报告
![DSP实验报告](https://img.taocdn.com/s3/m/ec05359dcd22bcd126fff705cc17552706225e50.png)
DSP实验报告实验⼀信号、系统及响应⼀、实验⽬的1.熟悉理想采样的性质,了解信号采样前后的频谱变化,加深对采样定理的理解。
2.熟悉离散信号和系统的时域特性。
3.熟悉线性卷积的计算编程⽅法:利⽤卷积的⽅法,观察分析系统响应的时域特性。
4.掌握序列傅⽒变换的计算机实现⽅法,利⽤序列的傅⽒变换对离散信号,系统及系统响应进⾏频域分析。
⼆、实验内容在编制以上各部分程序以后,编制主程序调⽤各个功能模块实现对信号,系统和系统响应的时域,频域分析,完成以下实验内容。
1.分析理想采样信号序列的特性。
产⽣理想采样信号序列()a x n ,使A =444.128,α=,0Ω= a.⾸先选⽤采样频率为1000Hz,T =1/1000,观察所得理想采样信号的幅频特性,在折叠频率以内和给定的理想幅频特性⽆明显差异,并做记录。
b.改变采样频率为300Hz,T=1/300,观察所得到的幅频特性曲线的变化,并做记录c.进⼀步减少采样频率为200Hz,T=1/200,观察频谱“混淆”现象是否明显存在,说明原因,并记录幅频特性曲线。
2.离散信号,系统和系统响应的分析(1).观察信号()b x n 和系统()b h n 的时域和幅频特性,利⽤线性卷积求信号通过系统以后的响应。
⽐较系统响应和系统()b h n 的时域和幅频特性,注意它们之间有⽆差异,并绘出曲线。
(2).观察信号()c x n 和系统()a h n 的时域和幅频特性,利⽤线性卷积求系统响应。
判断响应序列图形及序列⾮零值长度是否与理论结果⼀致,说出⼀种定性判断响应序列图形正确与否的⽅法(提⽰:10()()()c a x n h n R n ==)。
利⽤系统的傅⽴叶变换数值计算⼦程序求出()k j Y e ω,观察响应序列的幅频特性。
定性判断结果正确与否,改变信号()c x n 的脉冲宽度,使N =5,重复以上动作,观察变化,记录改变参数前后的差异。
(3)将实验步骤2-(2)中信号变换为()a x n ,其中A=1,α =0.4,0Ω=2.0734,T=1。
DSP实验报告.实验一 .离散信号的分析
![DSP实验报告.实验一 .离散信号的分析](https://img.taocdn.com/s3/m/fab5ed95b9d528ea80c7790a.png)
实验 步骤 及源 程序
n=-4:2;x=[1,-2,4,6,-5,8,10]; [x21,n21]=sigshift(5*x,n,-5); [x22,n22]=sigshift(4*x,n,-4); [x23,n23]=sigadd(x21,n21,x22,n22); [x2,n2]=sigadd(x23,n23,3*x,n); subplot(2,1,1);stem(n2,x2); xlabel('n');ylabel('x_1(n)');
DSP 实验报告书
姓名 学号 实验时间
实验 名称 实验 目的
实验一 离散信号的分析
1.观察常用离散时间信号的图形,掌握离散时间信号的基本序列运算。 2.理解离散时间系统的时域特性,加深对离散时间系统差分方程的理解。 3.熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解。
1.在给出的区间上产生并画出下面序列:
问题 讨论
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 程序—— function [y,n]=sigmult(x1,n1,x2,n2); m=[min(n1),min(n2)]; p=[max(n1),max(n2)]; n=min(m):1:max(p); 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; 函数程序: n=-4:2;x=[1,-2,4,6,-5,8,10];n4=-10:10;x3=2*exp(0.5*n4);x4=cos(0.1*pi); [x31,n31]=sigmult(x3,n4,x,n); [x32,n32]=sigshift(x,n,-2); [x33,n33]=sigmult(x4,n4,x32,n32); [x34,n34]=sigadd(x31,n31,x33,n33); subplot(2,1,1);stem(n34,x34); xlabel('n');ylabel('x_1(n)');
DSP实验1离散时间信号产生及频谱分析
![DSP实验1离散时间信号产生及频谱分析](https://img.taocdn.com/s3/m/5efffcc78bd63186bcebbc38.png)
实验1 离散时间信号产生及频谱分析一、实验目的㈠ 掌握MATLAB 产生常用离散时间信号的产生方法。
㈡ 掌握离散时间信号频谱的分析方法并观察其特点。
二、实验原理㈠ 常用的离散时间信号在 MATLAB 语言主要是研究离散信号的。
常用的离散信号有:1.单位取样序列⎩⎨⎧≠==0001)(n n n δ 2.单位阶跃序列⎩⎨⎧<≥=0001)(n n n u 3.实指数序列R a n a n x n∈∀=;)( 4.复指数序列n en x n j ∀=+)(0)(ωσ 5.正(余)弦序列)cos()(0θω+=n n x n ∀6.周期序列n N n x n x ∀+=)()(㈡ 离散信号的产生离散信号的图形显示使用stem 指令,连续信号的图形显示使用plot 指令。
㈢ 离散信号的频谱MATLAB 语言为我们提供了一些分析信号频谱的有效工具,其中使用FFT 子函数是一种方便快捷的方法。
三、实验任务任务1:编写MATLAB 程序来产生下列基本脉冲序列。
(1)单位脉冲序列:起点no ,终点nf ,在ns 处有一单位脉冲(no ≤ns ≤nf)。
(2)单位阶跃序列: 起点no ,终点nf ,在ns 前为0,在ns 处及以后为l(no ≤ns ≤nf)。
(3)实数指数序列:n ).(x 7503=(4)复数指数序列:n j e x )7.02.0(4+-=任务2:已知一时域连续正弦信号的频率为1Hz ,振幅值幅度为1V ,在窗口上显示一个周期的信号波形,对其进行32点采样后,进行32点的FFT ,观察其采样后的信号及信号频谱。
任务3:已知一个8点的矩形序列,用N=8点和N=32进行FFT 变换,作其时域信号图及信号频谱图。
并利用N=32的FFT 结果画其连续频谱。
DSP第一次实验报告
![DSP第一次实验报告](https://img.taocdn.com/s3/m/5aad437ff5335a8102d2209e.png)
四、实验要求
1、实验前,自己编制 C 语言程序。 2、实验中,自己调试程序,并准备验收时回答老师提出的各方面相关问题。 3、实验后,撰写实验报告和心得体会。 4、参考《数字信号处理实验教程》第四章(基础理论实验) 。
代码:
#include"stdarg.h" #include"graphics.h" #include"math.h" #include"stdio.h" int flag = 0; void conv(float x[],int M,float h[],int N,float y[]); void DO1(float m[],int n); void DO2(float m[],int n); void DO3(float m[],int n); void figureone(float x[],int m); void figuretwo(float x[],int m); void figurethree(float x[],int m); /*figure 函数初始化符号位*/
实验一 离散卷积的 C 语言编程实验 实验报告
一、实验性质
综合设计性实验
二、实验目的
1、认识常用的各种信号; 2、掌握卷积的定义和计算方法; 3、掌握在计算机中生成以及绘制信号序列图的方法。
三、实验内容
1、用 C 语言编程产生以下时间序列信号 ① x1(n)=n u(n) 0≤n≤31 ② x2(n)=an u(n) 0≤n≤31 0<a<1 ③ x3(n)= G(n) 0≤n≤31 ④ 产生自己欲实现的一个任意序列 (参考《数字信号处理实验教程》69 页序列 xa(n)=Ae-anTsin(w0nT)u(n), 0≤n≤31 ) 。 2、用 C 语言编程实现 y(n)=x1(n)*x3(n)和 y(n)=x2(n)*x3(n)并在 显示器上观察到x1(n) 、x2(n) 、x3(n)和 y(n)几个相应序列。
离散时间信号的产生及信号的卷积和运算实验报告2
![离散时间信号的产生及信号的卷积和运算实验报告2](https://img.taocdn.com/s3/m/b40adc330722192e4536f637.png)
离散时间信号的产生及信号的卷积和运算实验报告班级:___________ 姓名:__________ 学号:____________一、实验目的和原理实验原理:(一)DTFT 和DFT 的定义及其相互关系:序列x[n] 的DTFT 定义:∑=∞-∞=-n jn ωj ωx[n]e )X(e它是关于自变量ω的复函数,且是以π2为周期的连续函数。
)X(e j ω可以表示为:)(e jX )(e X )X(e j ωim j ωre j ω+=其中,)(eX j ωre 和)(e X j ωim 分别是)X(e j ω的实部和虚部;还可以表示为:)(ωj j ωj ωe )X(e )X(e θ=其中,)X(ej ω和}arg{)()X(e j ω=ωθ分别是)X(e j ω的幅度函数和相位函数;它们都是ω的实函数,也是以π2为周期的周期函数。
序列x[n]的N 点DFT 定义:∑∑-=-=-===10122][][)(][N n knNN n kn Njk NjW n x en x eX k X ππ][k X 是周期为N 的序列。
)X(e j ω与][k X 的关系:][k X 是对)X(e j ω在一个周期中的谱的等间隔N点采样,即:k Nj ω)X(e k X πω2|][==,而)X(e j ω可以通过对][k X 内插获得,即:]2/)1)][(/2([1)22sin()22sin(][1----=⋅--=∑N N k j N k j ωe Nk N kN k X N)X(e πωπωπω(二) 线性时不变离散时间系统的变换域表示:LTI 离散时间系统的时域差分方程为:∑∑==-=-Mk k Nk kk n x p k n y d)()((1) 传递函数:对上面的差分方程两边求z 变换,得:∑∑∑∑=-=-=-=-=⇒=Nk kkMk kkMk k k Nk kk z dzp z X z Y z p z X zd z Y 000)()()()(我们定义LTI 离散时间系统的输出的Z 变换Y(z)与输入的Z 变换X(z)的比值为该系统的传递函数,即)()()(z X z Y z H =为系统的传递函数。
数字信号处理实验一 时域离散信号的产生与基本运算
![数字信号处理实验一 时域离散信号的产生与基本运算](https://img.taocdn.com/s3/m/50383a621eb91a37f1115c8d.png)
实验一 时域离散信号的产生与基本运算一、实验目的1、了解常用的时域离散信号及其特点;2、掌握Matlab 产生常用时域离散信号的方法;3、掌握时域离散信号简单的基本运算方法;二、实验内容1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列;2、自己设定参数,分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、 信号和、信号积、信号能量。
3、已知信号⎪⎩⎪⎨⎧≤≤-≤≤-+=其他040614452)(n n n n x(1) 描绘)(n x 序列的波形;(2) 用延迟的单位脉冲序列及其加权和表示)(n x 序列; (3) 描绘一下序列的波形)2()(),2(2)(),2(2)(321n x n x n x n x n x n x -=+=-=三、实现步骤1、自己设定参数,分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、 实指数序列、随机序列。
输出图形如图1所示。
x=zeros(1,10); x(2)=1;subplot(3,2,1); stem(x,'filled')axis([0,10,-0.2,1]); title('单位抽样序列');N=10;u=ones(1,N); subplot(3,2,2); stem(u,'filled') axis([-10,10,0,1]); title('单位阶跃序列');x=-20:1:20;y=sin(0.2*pi.*x+0.5*pi); subplot(3,2,3);stem(x,y,'filled'); axis([-20,20,-2,2]); title('正弦序列');n=0:10; a1=1/2; y1=a1.^n;subplot(3,2,4);stem(n,y1,'filled');axis([0,10,0,1]);title('实指数序列,a=1/2');n=0:10;a2=2;y2=a2.^n;subplot(3,2,5);stem(n,y2,'filled');title('实指数序列,a=2');y=rand(1,20);subplot(3,2,6);stem(y,'filled');title('随机序列');图12、自己设定参数,分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、信号和、信号积、信号能量。
DSP数字信号处理实验报告(精品)
![DSP数字信号处理实验报告(精品)](https://img.taocdn.com/s3/m/e136c425376baf1ffd4fad12.png)
FFT 频谱分析一、 实验目的a) 进一步加深DFT 算法原理和基本性质的理解b) 熟悉FFT 算法原理和FFT 程序的应用c) 学习用FFT 对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确的应用FFT二、 实验原理a) 离散傅里叶变换(DFT ):离散傅里叶变换在作为有限长序列的傅里叶变换表示法在理论上相当重要;由于存在着计算离散傅里叶变换的快速算法(FFT ),从而离散傅里叶变换在各种数字信号处理的算法中起到了核心的作用。
其对应的离散傅里叶变换对为:X (K )=DFT [x (n )]=∑x (n )W N nk 0≤k ≤N −1N−1n=0x (n )=IDFT [X (K )]=1N ∑x (n )W N nk 0≤n ≤N −1N−1n=0需要注意:有限长序列的离散傅里叶变换及周期序列的离散傅里叶级数之间的关系是:它们仅仅是n 、k 的取值不同,DFT 只取主值区间。
X(n)、X(k)是一个有限长序列的离散傅里叶变换对,已知其中一个序列,就可以唯一确定另一个序列,这是因为x(n)、X(k)都是长为N 的序列,都有N 个独立值,所以信息量相同。
b) DFT 计算量:长度为N 的DFT 的计算量是N 个复数乘法和N-1个复数加法(4N 个实数乘法和4N-2个复数加法)c) FFT 的计算量:长度为N 的FFT 的计算量是 N 2log 2N 个复数乘法和N log 2N 个复数加法。
三、 实验步骤a) 复习DFT 的定义、性质和用DFT 做谱分析的有关内容b) 复习FFT 算法原理与编程思想,熟悉DIT-FFT 运算流图c) 编制信号产生程序,产生典型信号尽心谱分析。
d) 进行以下几个信号的谱分析i.x 1(n)=R 4(n ) ii. x 2(n)={n +1 ,0≤n ≤38−n ,4≤n ≤70 ,其他niii. x 3(n)= {4−n ,0≤n ≤3n −3 ,4≤n ≤70 ,其他niv. x 4(n )=cos π4n ,0≤n ≤19v.x 5(n )=sin π8n ,0≤n ≤19 vi.x 6(n )= cos 8πt +cos 16πt +cos 20πt vii.令x 7(n )=x 4(n )+x 5(n) N=8,16 viii. 令x 8(n )=x 4(n )+jx 5(n) N=8,16针对上述信号进行逐一的谱分析,下面给出针对各个信号的FFT 点数N 及对连续信号x 6(n )的采样频率f s ,供实验时参考 :x 1(n ),x 2(n ),x 3(n ),x 4(n ),x 5(n ) N=8,16x 6(n ) f s =64Hz ,N =16,32,64四、 实验内容a)对x1(n)=R4(n)进行谱分析1.编辑代码x1=[1 1 1 1];y11 = fft(x1,8);y12 = fft(x1,16);subplot(2,2,1);stem(0:3,x1);title('函数X1的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');2.谱分析图片b) 对x 2(n)={n +1 ,0≤n ≤38−n ,4≤n ≤70 ,其他n进行谱分析i. 编辑代码x2 = [1 2 3 4 4 3 2 1];y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));c)对x3(n)={4−n ,0≤n≤3n−3 ,4≤n≤70 ,其他n进行谱分析i.谱分析程序x3 = [4 3 2 1 1 2 3 4];y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));n ,0≤n≤19进行谱分析d)对x4(n)=cosπ4i.谱分析程序n = 0:1:19;x2 = cos(0.25*pi*n);y11 = fft(x2,32);y12 = fft(x2,64);subplot(2,2,1);stem(0:19,x2);title('函数X4的图像');subplot(2,2,2);stem(0:31,abs(y11));title('N=32的DFT');subplot(2,2,4);stem(0:63,abs(y12));title('N=64的DFT');ii.谱分析图片n ,0≤n≤19进行谱分析e)对x5(n)=sinπ8i.谱分析程序n = 0:1:19;x2 = sin(0.125*pi*n);y11 = fft(x2,32);y12 = fft(x2,64);subplot(2,2,1);stem(0:19,x2);title('函数X5的图像');subplot(2,2,2);stem(0:31,abs(y11));title('N=32的DFT');subplot(2,2,4);stem(0:63,abs(y12));title('N=64的DFT');ii.谱分析图片f)对x6(n)=cos8πt+cos16πt+cos20πt进行谱分析i.谱分析程序n = 0:1:15;x1 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);n = 0:1:31;x2 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64); n = 0:1:63;x3 = cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);y1 = fft(x1,16);y2 = fft(x2,32);y3 = fft(x3,64);subplot(3,2,1);stem(0:15,x1);title('函数X6 N=16 的图像');subplot(3,2,2);stem(0:15,abs(y1));title('N=16的DFT');subplot(3,2,3);stem(0:31,x2);title('函数X6 N=32 的图像');subplot(3,2,4);stem(0:31,abs(y2));title('N=32的DFT');subplot(3,2,5);stem(0:63,x3);title('函数X6 N=64 的图像');subplot(3,2,6);stem(0:63,abs(y3));title('N=64的DFT');ii.谱分析图片g)对x7(n)=x4(n)+x5(n)进行谱分析i.谱分析程序n = 0:1:19;x2 = j*sin(0.125*pi*n) + cos(0.25*pi*n);y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:19,x2);title('函数X=X4 + X5 的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图片h)对x8(n)=x4(n)+jx5(n)进行谱分析i.谱分析程序n = 0:1:19;x2 = j*sin(0.125*pi*n) + cos(0.25*pi*n);y11 = fft(x2,8);y12 = fft(x2,16);subplot(2,2,1);stem(0:19,x2);title('函数X=X4 + X5 的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图片五、实验分析六、实验结论通过这次利用FFT对信号进行频谱分析的实验,更加深刻的理解了DFT算法的理解和性质的理解;同时也更为熟悉了FFT算法的原理和应用;学会使用FFT 对离散信号和连续信号进行频谱分析,了解了可能出现的分析误差和原因。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
离散信号的产生及运算报告
一、实验目的:
1、复习和巩固数字信号处理中离散信号的产生和运算
2、学习和掌握用MATLAB产生离散信号的方法
3、学习和掌握用MATLAB对离散信号进行运算
二、实验原理:
1.用MATLAB函数产生离散信号
信号是数字信号处理的最基本内容。
没有信号,数字信号处理就没了工作对象。
MATLAB7.0内部提供了大量的函数,用来产生常用的信号波形。
例如,三角函数(sin,cos),指数函数(exp),锯齿波函数(sawtooth), 随机数函数(rand)等。
1 产生被噪声污染的正弦信号
用随机数函数产生污染的正弦信号。
2 产生单位脉冲序列和单位阶跃序列
按定义,单位脉冲序列为
单位阶跃序列为。
3 矩形脉冲信号:在MATLAB 中用rectpuls 函数来表示,其调用形式为:
y=rectpuls(t,width,用以产生一个幅值为1,宽度为width,相对于t=0 点左右对称的矩形波信号,该函数的横坐标范围由向量t 决定,是以t=0 为中心向左右各展开width/2 的范围,width 的默认值为1。
例:以t=2T(即t-2×T=0为对称中心的矩形脉冲信号的MATLAB 源程序如下:(取T=1)
t=0:0.001:4;
T=1;
ft=rectpuls(t-2*T,2*T;
plot(t,ft;grid on; axis([0 4 –0.5 1.5];
4 周期性矩形波(方波)信号在MATLAB 中用square 函数来表示,其调用形
式为:y=square(t,DUTY,用以产生一个周期为2π、幅值为±1 的周期性方
波信号,其中的DUTY参数表示占空比,即在信号的一个周期中正值所占的百
分比。
例如频率为30Hz 的周期性方波信号的MATLAB 参考程序如下:t=-0.0625:0.0001:0.0625;
y=square(2*pi*30*t,75;
plot(t,y;axis([-0.0625 0.0625 –1.5 1.5];grid on ;
2、 MATLAB中信号的运算
乘法和加法:离散信号之间的乘法和加法,是指它的同序号的序列值逐项对应相乘和相加。
三、实验步骤及结果:
1.产生噪声污染的正弦信号
波形如下所示:
2.产生以下信号,并作图
作图结果如下示:
3、信号的基本运算
实验结果如下:
四、实验总结:
通过实验,将随机信号、各种基本信号及信号之间的基本运算通过图像展示出来,从而使我直观的巩固了课堂基础知识。
另外,在使用matlab的过程中,在体会到软件自身的强大性的同时,发现自己的操作还不够熟练,今后在这方面要加强。
附:实验源程序:
1、噪声污染的正弦信号
t=0:0.0001:1;
f=5;
y1=sin(2*pi*f*t;
y2=randn(size(t;
y=y1+y2;
plot(t,y
xlabel('时间t'
ylabel('幅值'
title('噪声污染正弦信号' 2、产生基本信号
t1=0:0.001:1;
A1=5;
f=3;
y1=A1*sin(2*pi*f*t1;
subplot(411
plot(t1,y1
title('正弦信号波形'
t2=-6:0.001:6;
A2=13;
y2=A2*rectpuls(t2,2; subplot(412
plot(t2,y2
title('矩形脉冲信号'
t3=-8:0.001:8;
y3=sawtooth(2*pi/8*t3,0.70; subplot(413
plot(t3,y3
title('三角波脉冲'
t4=-6:0.001:6;
y4=square(2*pi/6*t4,30;
subplot(414
plot(t4,y4
title('方波信号'
3、信号基本运算
①t=-5:0.001:5;
y=sin(pi*t/5+4*cos(pi*t/4;
plot(t,y
title('正弦叠加信号'
xlabel('时间t'
②k=1:1:500;
fs=1000;
x1=sin(2*pi*100*k/fs;
x2=sin(2*pi*10*k/fs;
y=x1.*x2;
subplot(311
stem(k,x1
title('x1的函数图像' xlabel('时间t' subplot(312
stem(k,x2
title('x2的函数图像' xlabel('时间t' subplot(313
stem(k,y
title('y的函数图像' xlabel('时间t'。