信号与系统实验报告,（范例大全）【精彩4篇】

【路引】由阿拉题库网美丽的网友为您整理分享的“信号与系统实验报告,（范例大全）【精彩4篇】”文档资料，以供您学习参考之用，希望这篇范文对您有所帮助，喜欢就复制下载支持吧！

信号与系统仿真实验报告【第一篇】

信号与系统

仿真

实 验 报 告

班级： 学号： 姓名： 学院：

实验一

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

5(s25s6)求三阶系统H(s)3的单位阶跃响应，并绘制响应波形图。

s6s210s8

五、解题分析：要知道求单位阶跃响应需知道所用函数，以及产生波形图所需要用到的函数。

六、试验程序：

num=[5 25 30]； den=[1 6 10 8]； step(num,den,10)； title（‘Step response’）

七、实验结果：

实验所得波形图如下：

Step (sec)678910

八、实验心得体会：通过本次试验了解学会了一些新的函数的应用。了解到了N阶系统的单位阶跃响应的计算方法，和系统的响应波形图的函数应用和绘制方法。为后面的实验打下基础，并对信号仿真和《信号与系统》这门课程之间的联系有所增加，对《信号与系统》这门课里的问题也有了更加深入地了解。

九、实验改进想法：无。

实验二

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目:

一个因果线性移不变系统y（n)(n2)x(n)x(n2)，求：（1）H(z)；（2）冲激响应h(n)；（3）单位阶跃响应u(n)；（4）H(ej），并绘出幅频和相频特性。

五、解题分析：离散卷积是数字信号处理中的一个基本运算，MTLAB提供的计算两个离散序列卷积的函数是conv，其调用方式为 y=conv(x,h) 。其中调用参数x,h为卷积运算所需的两个序列，返回值y是卷积结果。

MATLAB函数conv的返回值y中只有卷积的结果，没有y的取值范围。由离散序列卷积的性质可知，当序列x和h的起始点都为k=0时，y的取值范围为k=0至length(x)+length(h)-2。

许多离散LTI都可用如下的线性常系数的差分方程描述

ay[kn]bx[kn]

nnn0n0NN其中x[k]、y[k]分别系统的输入和输出。在已知差分方程的N个初始状态y[k]，和输入x[k]，就可由下式迭代计算出系统的输出

y[k](an/a0)y[kn](bn/b0)x[kn]

n1n0NM利用MATLAB提供的filter函数，可方便地计算出上述差分方程的零状态响应。filter函数调用形式为 y=filter(b,a,x) 。其中 a[a0,a1,。.。，aN]，b[b0,b1,。.。，bM] ，分别表示差分方程系数。X表示输入序列，y表示输出序列。输出序列的长度和序列相同。

当序列的DTFT可写成ej的有理多项式时，可用MATLAB信号处理工具箱提供的freqz函数计算DTFT的抽样值。另外，可用MATLAB提供的abs、angle、real、imag等基本函数计算 DTFT的幅度、相位、实部、虚部。若X（ej）可表示为

b0b1ej。.。bMejMB（ej） X(e)jjjNA(e)a0a1e。.。aNe则freqz的调用形式为 X=freqz(b,a,w) ，其中的b和 a分别是表示前一个

j式子中分子多项式和分母多项式系数的向量，即a[a0,a1,。.。，aN] ，w为抽样的频率点，向量w的长度至少为2。返回值X就是DTFTb[b0,b1,。.。，bM]。在抽样点w上的值。注意一般情况下，函数freqz的返回值X是复数。

六、实验程序:

clc;clear;close; b=[1 0 -1]； a=[1 0 -]； figure(1)； subplot(2,1,1)； dimpulse(b,a,20) subplot(2,1,2)； dstep(b,a,50) w=[0:1:512]*pi/512; figure(2)； freqz(b,a,w)

七、实验结果：

冲击响应图及阶跃响应图:

Impulse (sec)Step (sec)3035404550 100Magnitude (dB) Frequency （ rad/sample） (degrees) Frequency （ rad/sample）

八、实验心得体会：通过实验我们知道了使用Matlab来绘出出一个线性移不变系统的幅频和相频曲线。并知道了在《信号与系统》中得一些差分方程和各种响应，譬如零输入相应、零状态响应、全响应、自由响应、强迫响应、冲击响应、单位阶跃响应等等各种响应在Matlab中的函数表达方式和他们的求法，以及系统的幅频和相频曲线的绘制都有了一定深刻的认识。

九、实验改进想法：无。

实验三

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

模拟信号x(t)2sin(4t)5cos(8t)，求N64的DFT的幅值谱和相位谱。

五、解题分析：在MATLAB信号处理工具箱中，MATLAB提供了4个内部函数用于计算DFT和IDFT，它们分别是：fft(x)，fft(x,N)，ifft(X)，ifft(X,N)。

fft(x) 计算M点的DFT。M是序列x的长度，即M=length(x)。

fft(x,N) 计算N点的DFT。若M>N，则将原序列截短为N点序列，再计算其N点DFT；若M

ifft(X) 计算M点的IDFT。M是序列X的长度。

ifft(X,N) 计算N点IDFT。若M>N，则将原序列截短为N点序列，再计算其N点IDFT；若M

六、实验程序：

clc;clear;close; N=64; n=0:63; t=d*n; q=n*2*pi/N; x=2*sin（4*pi*t)+5*cos(8*pi*t)； y=fft(x,N)； subplot(3,1,1)； plot(t,x)； title（‘source signal’）； subplot(3,1,2)； plot(q,abs(y)）； title（‘magnitude’）； subplot（3,1,3)； plot(q,angle(y)）； title（‘phase’）；

七、实验结果：

***0100806040200|F(k)|05101520Frequency253035

***0100806040200|F(k)|05101520Frequency253035 4321|jW|0-1-2-3-405101520Frequency253035Step Response400020000-2000 Amplitude-4000-6000-8000-10000-12000-1400001234n (samples)5678

八、实验心得体会：通过本次试验我知道了求取模拟信号在N等于一定值时的的DFT的幅值谱和相位谱的求法。通过本次实验，对幅值谱和相位谱有了更深的了解，并与课程《信号与系统》里的一些相关知识连接到了一起，使得学到的只是更加深刻、有意义。

九、实验改进想法：无。

实验四

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

将信号x(t)sin(240t)做离散序列，比较原序列与经过FFT和IFFT变换后的序列，并做出说明。

五、解题分析：此题需要对信号做离散序列，还要做FFT和IFFT变换，然后得到图像进行比较。连续时间函数与离散时间函数在编程中的区别主要体现在如下两个方面：第一，自变量的取值范围不同，离散时间函数的自变量是整数，而连续时间函数的自变量为一定范围内的实数；第二，绘图所用的函数不同，连续函数图形的绘制不止一个。本实验中要求绘制离散时间信号图，可以应用MATLAB中的函数来实现。用MATLAB表示一离散序列，可用两个向量来表示。其中一个向量表示自变量的取值范围，另一个向量表示序列的值。之后画出序列波形。当序列是从0开始时，可以只用一个向量x来表示序列。由于计算机内寸的限制，MATLAB无法表示一个无穷长的序列。对于典型的离散时间信号，可用逻辑表达式来实现不同自变量时的取值。

六、实验程序：

t=0:1/255:1; x=sin（2*pi*120*t)； y=real(ifft(fft(x)）)； subplot(2,1,1)； plot(t,x)； title（‘原波形’）； subplot(2,1,2)； plot(t,y)；

七、实验结果：

原波形恢复的波形

八、实验心得体会：通过对做信号的离散序列以及经FFT和IFFT的变换，了解了相关特性。通过计算机做出的信号波形图，我们能够很直白的看出原波形和经过变换后的波形的差别。

九、实验改进想法：无。

实验五

一、实验者姓名：

二、实验时间：

三、实验地点：

四、实验题目：

2s，激励信号22(s1)100x(t)(1cot)sco1s0(t)0，求（1）带通滤波器的频率响应；（2）输出稳态响应并绘制图形。 已知带通滤波器的系统函数为H(s)

五、解题分析：需要知道求频率响应的方法，并绘制图形。

六、实验程序：

clear; t=linspace（0,2*pi,1001)； w=[99,100,101]； U=[,1,]； b=[2,0]； a=[1,2,10001]； u1=U*cos(w’*t+angle(U’）*ones（1,1001)）； H=polyval（b,j*w)。/polyval(a,j*w)； H=freqs(b,a,w)； subplot(2,1,1)，plot(w,abs(H)），grid; subplot（2,1,2)，plot(w,angle(H)），grid; u21=abs（U(1)*H(1)）*cos（99*t+angle(U(1)*H(1)）)； u22=abs（U(2)*H(2)）*cos（100*t+angle(U(2)*H(2)）)； u23=abs（U(3)*H(3)）*cos（101*t+angle(U(3)*H(3)）)； u2=u21+u22+23; figure(2)； subplot(2,1,1)，plot(t,u1)； subplot(2,1,2)，plot(t,u2)；

七、实验结果：

210-1-***222101234567

八、实验心得体会：通过本次试验，了解了频率响应求法，加深了对输出稳态响应的印象。

九、实验改进想法：无。

信号与系统实验报告【第二篇】

中南大学

信号与系统试验报告

姓名:

学号:

专业班级:自动化

实验一

基本信号得生成1.实验目得

 学会使用 MATLAB 产生各种常见得连续时间信号与离散时间信号；  通过ＭATLAB 中得绘图工具对产生得信号进行观察，加深对常用信号得理解；  熟悉 MＡTＬAB 得基本操作，以及一些基本函数得使用，为以后得实验奠定基础。2 。

实验内容

⑴ 运行以上九个例子程序，掌握一些常用基本信号得特点及其ＭＡＴＬAB 实现方法；改变有关参数，进一步观察信号波形得变化。⑵ 在范围内产生并画出以下信号: a)

； b) ； c) ； d) 。源程序: k= — 10 0 ：1 1 0; ；

f1k=[ze ｒ os(1 ， 10 ）， ， 1 ，zer ｏｓ（1 1 ，1 1 0）］； ；

subplot （2 2 ， 2,1)

stem( ｋ， ， ｆ 1k ）

title( ＇ f1[k ］ ’ ）

f2k=[zer ｏs s （1 ，8)， 1，z z ｅr r ｏs s （1 1 ， 12 ）］； ；

su ｂ plot(2 ， 2,2 ）

ｓ te ｍ（k k ， f2k ）

titl ｅ (“f2 ［ｋ］ ’ ）

f3k ＝ [zeros （1 ，14 ）， ， 1， zer ｏs s （1 1 ，6 6 ）］； ；

su ｂｐｌｏt t （2 2 ，2 2 ， 3)

st ｅ m（ k,f 3 k）

ti ｔ le(”f3[k]“ ）

ｆ 4k= 2＊ｆ2k k- - ｆ3 k;

ｓｕ bpl ｏｔ (2 ，2 2 ， 4)

s s ｔ em （ k,f4k)

t t ｉ tle（ ＇ｆ 4[k]”）

⑶ 在 范围内产生并画出以下信号: a) ； b) ； c) 。

请问这三个信号得基波周期分别就是多少？ 源程序: ｋ= = 0：

31 ；

ｆ1 1 ｋ= = ｓｉｎ（ｐｉ /4* ｋ) ) 、＊ cos （ｐｉ /4*k ）； ；

ｓｕｂｐ lo ｔ （3,1, 1) ）

ｓｔ em（k ，ｆ1 1 ｋ） )

ｔ itle( ＇ f1[k ］“ ” ）

ｆ2 k=（cos(pi/4*k)） 、^ ^ 2; ；

subp ｌ ot(3, 1 ，2)

ｓt t ｅm m （ｋ， f2 ｋ) )

ｔi i ｔｌ e( ＇ｆ2 [k ］“ ” ）

f3 ｋ =s ｉn n （ pi ／4＊ｋ）、＊ｃｏs s （ｐ i/8 ＊k k ）； ；

sub ｐ lot （3，1 1 ，3）

st ｅm m （ｋ ，f3k ）

tit ｌe e （ ’f3 ［k k ］ ’ ）

其中ｆ1[k]得基波周期就是4, f2[k］得基波周期就是4, ｆ3［k］得基波周期就是 16． 实验 二

信号得基本运算1 。

实验目得

 学会使用ＭATLAＢ完成信号得一些基本运算；  了解复杂信号由基本信号通过尺度变换、翻转、平移、相加、相乘、差分、求与、微分及积分等运算来表达得方法；  进一步熟悉 MＡTLAB 得基本操作与编程，掌握其在信号分析中得运用特点与使用方式． 2。

实验内容

⑴ 运行以上三个例题程序，掌握信号基本运算得ＭATＬAＢ实现方法；改变有关参数，考察相应信号运算结果得变化特点与规律。

⑵ 已知信号如下图所示：

a) 用 MAＴLAB 编程复现上图； ％作业题2

a: t=-6:0 、001 ：6; ft1=t ｒip ｕｌs(t, 6,0 、5）； sub ｐｌot( 2,1,1 ）

ｐｌot（t,ft1) ｔ ｔit ｌe（ ’f(t） ’）

b) 画出得波形； %b t= －6:0 、00 1:6; ｆ ｆt 1= ｔripuls(2 ＊(1 —t)，6 ，0 、5) ； %s ｕbp ｌot(1,1,1 ）

pl ｏt(t ，ｆt 1) t ｉt ｌｅ（ ’f （2＊(1- ｔ）“） -4 -3 -2 -1 0 1 2 3 (t)给 定 信 号 f(t)

c) 画出得波形；

%c h=0 、00 1；t= —6: ｈ:6; yt= ｔri ｐu ｌs（ ｔ，6,0 、5） ； y 1= ｄiff(yt ）*1/h; plot(t（ 1:lengt ｈ(t） —1 ），ｙ1 ）

tit ｌe（ ’df(t ）/ ｄt ’）

d) 画出得波形。

％d t＝—6:0、1:6; for x=1:ｌｅngｔh(t)

ｙ2（x)=ｑuａd（’ｔｒipuｌs(ｔ，6，0、5）’，－3,t(x)）；

end ｐlot(ｔ，y2）

tｉtｌe（”iｎtｅgｒaｌ ｏf f(t）“）

实验 三

系统得时域分析1 。

实验目得

 学习并掌握连续时间系统得零状态响应、冲激响应与阶跃响应得 MATLAB求解方法；  学习并掌握离散时间系统得零状态响应、冲激响应与阶跃响应得 MＡTＬAB 求解方法；  进一步深刻理解连续时间系统与离散时间系统得系统函数零极点对系统特性得影响； 学习并掌握卷积得 MＡTLAB 计算方法。

2、实验内容

⑴ 运行以上五个例题程序，掌握求解系统响应得 MＡTＬAＢ分析方法；改变模型参数，考察系统响应得变化特点与规律。

⑵ 设离散系统可由下列差分方程表示：

计算时得系统冲激响应。

源程序：

ｋ＝ — 20:100 ；

a= ［1 1

- - 1 0、9] ；

b=[1]；

h= ｉ mpz （b b ，a a ，k k ）； ；

stem （ k, ｈ）； ；

xla ｂ el（’Ti ｍｅ （sec）’）

y y ｌａｂ el（’y（t）”）

⑶ 设，输入，求系统输出。

（取）

源程序: k=- -1 1 0 :50 ；

u u ｋ= = ［z z ｅ ro ｓ（ 1,1 0 )， ｏ nes（1, 51) ） ］； ；

ｕ 1k=[ ｚ er ｏ s( 1， 20 ），o o ｎ es( 1, ， 41）］； ；

ｈ k=0 、9 9 、^k 、*uk;

fk=u ｋ- - u1k;

yk=co o ｎv v （ hk,f ｋ) ) ；

stem （ 0:length(yk ）- - 1， yk ）； ；

⑷ 已知滤波器得传递函数：

输入信号为为随机信号。试绘出滤波器得输出信号波形。(取) 源程序: R=101 ；

d=rand （1 ，R ）

—0 0 、5; ；

t=0:100 ；

s=2 ＊ si ｎ (0 、05*pi*t) ；

f=s ＋d d ；

su ｂｐ lo ｔ (2,1,1)；

plot（t ，ｄ ，’ ｇ- - 、＇，t t ， s,’b- - — “，t,f ， ”r- - ＇ ）；

ｘl l ａb b ｅl l （“ ” Ｔi i ｍ e in ｄ ex t’ ）； ；

legend （ “d ［t t ］ ”，“ ｓ[ [ ｔ］” “ ， ”f[t ］ ’ ）；

tit t ｌe e （“ ” 处理前得波形＇) )

b=[0 、22 0 ]；a=[ 1

－0 0 、8]；

y=fi ｌｔ er（b ，ａ，ｆ） ) ；

su ｂp p ｌ ot （2 2 ， 1,2)；

pl ｏt t （ｔ ，s,“b —- -” “ ，t t ， y,’r- - ’) ；

xl ａb b ｅ l（’ Ｔi i ｍ e i ｎｄ ex t”）；

leg ｅ nd(“s ［t t ］ ’，’ ｙ [t ］＇）； ；

title （” “ 滤波器输出波形’) )

实验 四

周期信号得频域分析

1、。实验目得

 掌握周期信号傅立叶级数分解与合成得计算公式  掌握利用 MATＬAＢ实现周期信号傅立叶级数分解与综合方法  理解并掌握周期信号频谱特点

2、实验内容1、仿照例程，实现下述周期信号得傅立叶级数分解与合成:

要求:

(ａ）首先，推导出求解，，得公式，计算出前 10 次系数；

(b)利用ＭＡTＬAB 求解，，得值，其中，求解前 10 次系数，并给出利用这些系数合成得信号波形。

（ａ)设周期信号得周期为，角频率，且满足狄里赫利条件，则该周期信号可以展开成傅立叶级数。

（1）三角形式傅立叶级数

dt t n t fTbdt t n t fTadt t fTat n b t n a at b t a t b t a t b t a a t fTT nTT nTTnnnnn n n n           ***1 02 2 2 2 1 1 1 1 0111111sin ) （2cos ） （2) (1) sin（ ） cos(。.。sin cos 。.。sin cos sin cos ） (      

（2）指数形式傅立叶级数

（b）求解，,及合成信号波形所用程序: fｕnction [Ａ_sｙm,B_sym]＝ＣTFShｃhsym ％ 采用符号计算求一个周期内连续时间函数 f 得三角级数展开系数，再用这些 %

展开系数合成连续时间函数ｆ、傅立叶级数 ％ 量值数是就都出输入输得数函ﻩ%

数阶得波谐 6=fNﻩ% Nn

数位确准得据数出输ﻩ% 数系开展项 soc 波谐次、、、3,2,1是就次依素元后其，项流直是就素元 1 第ﻩmyｓ_Aﻩ% B_sym 第 2，3，4,、、、元素依次就是 1，2，3、、、次谐波 siｎ项展开系数 %

taｏ=1

tａo/T=0、2 sｙｍs t n k x

T＝4; tao=T/4; ａ=-1、5; ｉf ｎarｇｉn<4

Nf=10； enｄ if nａrｇｉn＜5

Ｎｎ=32; end

1 -3 -4 5 4 1 O

x＝timｅ＿ｆun_x（t)； A0=iｎt(x,t，a,Ｔ+a）/Ｔ；

％求出三角函数展开系数Ａ0 Ａｓ＝2/T*iｎｔ（x＊coｓ（2*pi*ｎ*t/Ｔ)，t,a，T+ａ)；

％求出三角函数展开系数 As Bｓ=2/Ｔ*ｉnｔ（x*sin（2*pｉ＊n＊t/T)，t，a，T+ａ）；

%求出三角函数展开系数 Bs Ａ＿sym（1)=ｄoubｌe（vｐａ(A0,Ｎn））；

％获取串数组 A0 所对应得 ASC2码数值数组 fｏr k＝1:Nf A＿sym(k＋1)＝ｄouble（ｖpａ(subｓ（Ａs，n，k），Nn））；

％获取串数组Ａ所对应得 ASC2码数值数组 B_sym（k+1)=ｄｏｕble(vpa（subs(Ｂs，ｎ，k），Nｎ））；

%获取串数组Ｂ所对应得 AＳC2 码数值数组 end；

if nａrgoｕｔ==0

c=A＿syｍ； disp(ｃ）；

%输出 c 为三角级数展开系数:第 1 元素就是直流项，其后元素依次就是 1,2,3、、、次谐波ｃos 项展开系数 d=B＿ｓym; disp(ｄ)；

%输出 d 为三角级数展开系数：

第 2,3，4,、、、元素依次就是 1,2,3、、、次谐波ｓｉn 项展开系数

ｔ=—3＊T:0、01:3*T;

ｆ 0= ｃ(1)；

%直流

ｆ 1 ＝ c （ 2) 、＊ co ｓ （ 2* ｐ i* 1 * ｔ /T) ＋ d(2 ）、＊ s ｉ n(2 ＊ pi* 1 * ｔ /Ｔ）；

% 基波

ｆ 2= ｃ （ 3 ） 、* ｃ o ｓ （ 2*pi ＊ 2 ＊ t/T)+d(3 ）、*sin （ 2 ＊ pi ＊ 2* ｔ ／T）；

% 2 次谐波

f3=ｃ(4）、＊ｃos(2＊pi*3＊t／T)+d(4)、＊sin(2*ｐi*3＊t/T)；

% 3次谐波

f4=ｃ（5)、＊cos（2*ｐi*4＊t/T）+d（5)、＊siｎ（2*pi＊4*t/Ｔ）；

% 4次谐波

f5=c（6）、*ｃos(2*ｐi*5*ｔ/T）+d(6）、*ｓｉn(2＊ｐi*5*t／T)；

% 5次谐波

f6=c（7)、*cｏｓ（2*pi*6＊t/Ｔ）+d(7)、*sin（2＊pｉ*6*t/T)；

％ 6 次谐波

ｆ 7=c(8 ）、*cos(2*p ｉ ＊ 7 *t/T ）

+d(8 ）、＊ sin(2 ＊ ｐ ｉ ＊ 7 ＊t/T)；

% 7 次谐波

f8=c(9)、＊ｃos（2*ｐi＊8*ｔ/T)+d(9）、*sin（2*ｐi＊8＊t/T）；

% 8次谐波

f9 ＝ ｃ （1 0 ） 、＊ c ｏ ｓ （ 2 * ｐ i*9 ＊ ｔ /T）+d(10) 、* ｓ in(2 ＊ ｐ i ＊ 9 ＊t/T）；

% 9 次谐波

f 1 0=c （ 11) 、*co ｓ ( 2 ＊ pi*10*t/T ）

+d（1 1 ） 、*s ｉ n(2*pi ＊ 1 0 ＊t/T)；

％ 10次谐波

ｆ11=f0+f1+f2；

% 直流+基波＋2 次谐波

f12=f11+f3;

％ 直流+基波+2 次谐波+3 次谐波

f13＝f12＋f4+f5+f6；

% 直流＋基波+2 次谐波+3次谐波+4次谐波+5 次谐波+6 次

谐波

f14＝f13+f7+f8+ｆ9+f10；

%0～10 次

subｐlot(2,2，1)

ｐlｏt(t，f0＋f1），hoｌｄ ｏn

y＝tiｍe_ｆun＿ｅ(t)；

％调用连续时间函数－周期矩形脉冲

pｌoｔ(t，y，”r：“）

titｌe（”直流+基波’)

axiｓ（[-8，8，-0、5，1、5]）

sｕbploｔ（2,2,2）

plot(t,f12)，hold oｎ

y=time_fun_e（t）；

ｐlot(t,y，’r：’）

tiｔle(“1—3 次谐波+直流”）

axis([—8,8,-0、5,1、5］）

sｕｂploｔ(2,2,3)

ploｔ(t,f13），hoｌd on

y＝ｔime_fun＿e（t)；

plot(ｔ，y,’ｒ:’）

ｔitle（“1—6 次谐波＋直流＇)

axis（［-8，8，－0、5,1、5］）

ｓubpｌot（2,2，4)

plｏt（t,ｆ14），holｄ on

y=timｅ_fun_e（t）；

ploｔ（t，y,”r:’）

tｉtlｅ（’1—10 次谐波+直流“）

axｉs（[－8，8,-0、5，1、5］)

hｏｌd off ｅnd fｕnｃtiｏn y=time_fun_e(t) % 该函数就是 CＴFShｃhsym、m 得子函它由符号函数与表达式写成 a＝1、5； T=4； h＝1； taｏ=Ｔ/4; t=—3*Ｔ:0、01:3*T; ｅ1＝1/2+1/2、＊sｉgｎ（ｔ—0、5＋tａo／2）； e2＝1／2+1/2、*sign(t—0、5—tａo／2)； y=h、＊(e1—e2）；

%连续时间函数—周期矩形脉冲 funcｔion x＝tiｍｅ_fun＿ｘ（t）

% 该函数就是 CTFSｈcｈsym、ｍ得子函数。它由符号变量与表达式写成。h=1;

ｘ1=sｙm（”Heaｖisidｅ（ｔ）“）*h; x＝x1－sym(’Heaｖisidｅ（t-1）’）＊ｈ；

2、已知周期为Ｔ=4 得三角波，在第一周期(-2

functｉon [A＿sym,B_sym]=CTFSｓhbpsyｍ（T，Nｆ） % 采用符号计算求［0，T]内时间函数得三角级数展开系数。

%

ﻩ 函数得输入输出都就是数值量 ％ Ｎｎ

输出数据得准确位数 % mys_Aﻩ 第1元素就是直流项，其后元素依次就是1,2,3、、、次谐波 coｓ项展开系数 ％ Ｂ_ｓyｍ 数系开展项 nis 波谐次、、、3,2,1 是就次依素元、、、，4，3,2 第ﻩ％

T

T=m*taｏ， 信号周期 ％ ﻩ ｆNﻩ 谐波得阶数 ％

m (m=T/tao)周期与脉冲宽度之比，如 m=4，8，16,100等 ％

tao

脉宽：ｔａo＝Ｔ/ｍ

sｙｍs t

n

ｙ

iｆ naｒgin＜3

Nｆ=inｐut（’plｅａse Inｐut 所需展开得最高谐波次数:Nf＝’)； end

T=ｉnpuｔ(’ｐleaｓe Iｎput 信号得周期 T＝”）； if nａrｇｉn〈5

Ｎn=32; ｅnd y=tiｍe_ｆun_s(t)； A0＝2/Ｔ＊int（ｙ，t，0,T／2)； Ａs=2/Ｔ＊iｎｔ（y*cｏs（2*pi*n＊t/Ｔ），ｔ，0,T/2)；

Bｓ＝2/T＊iｎt（ｙ*ｓin(2＊pi*n＊t/Ｔ），t,0,Ｔ／2）； Ａ＿sym(1）=double(vpa(Ａ0，Ｎｎ））； for k=1：Ｎf

Ａ_sym(k+1)＝double(vpa（subs（Ａs,n，k），Nn））；

B＿sym（k+1)=ｄouble（ｖｐa(sｕbs(Bs，n，k），Ｎn)）；

end if nａrgout==0

Ａn=flｉｐlr（Ａ_sym）；

%对 A_sym 阵左右对称交换

An（1，k+1）=A_sｙm(1)；

%Ａ＿sym 得 1*k 阵扩展为 1＊（k+1)阵

An=ｆｌｉplr（Aｎ)；

%对扩展后得 S1阵左右对称交换回原位置

Bn＝fｌiｐlr（B_ｓym)；

%对 B_sym 阵左右对称交换

Ｂn（1,ｋ+1)=0;

%B＿sｙｍ得 1＊k 阵扩展为1*(k+1)阵

Bn=fliｐlr（Bｎ)；

%对扩展后得 S3阵左右对称交换回原位置

FnR=Aｎ/2—i＊Ｂn/2；

% 用三角函数展开系数 A、Ｂ值合成付里叶指数系数

ＦnL=fｌｉplr（FnＲ）；

N=Nf*2*pi/T;

k2=—N:2＊pi／T:N;

Fn＝[FnＬ，FnＲ（2:eｎｄ）]；

%subploｔ（3，3，3）

%x＝tｉmｅ_fuｎ_ｅ(ｔ)；

% 调用连续时间函数-周期矩形脉冲

sｕｂplot(2,1，1)

stｅm(k2,aｂs(Fn)）；

％画出周期矩形脉冲得频谱(T=M*tao）

ｔiｔle（＇连续时间函数周期三角波脉冲得双边幅度谱’）

axｉs([－80，80,0,0、12]）

ｌine（[-80，80］，[0,0]，＇colｏr“，’ｒ”）

ｌｉｎe([0,0]，[0，0、12]，’colｏｒ’，＇r“）

end ｆunctioｎ x=time＿fun_e(t) % 该函数就是ＣTFSｓhbpsｙm、ｍ得子函数。它由符号变量与表达式写成。

％ t

组数间时是就ﻩ％ T 2、0=T/oat=yｔuｄ

期周是就ﻩﻩT=5; t=—2＊T:0、01：2*T； tao＝T/5; ｘ=rectpuｌs（t,tａo)；

％产生一个宽度 tａo=1 得矩形脉冲 sｕbplot(2,2,2）

ｐlｏt(ｔ，x）

ｈoｌｄ on x=rectｐuls(t—5，tao)；

％产生一个宽度ｔao=1 得矩形脉，中心位置在ｔ＝5处 pｌot（t,ｘ)

hoｌd ｏn x=rｅctpｕls(t＋5,ｔaｏ）；

%产生一个宽度ｔao＝1得矩形脉，中心位置在 t=—5处 ploｔ(t，x)

ｔitlｅ（”周期为 T＝5,脉宽 tａo=1得矩形脉冲＇） ａxis（[-10,10,0,1、2］） funｃtｉoｎ y＝time_fｕｎ＿s(t）

syms t y＝1—abs（t)； ｘ1=sym（＇Hｅavｉｓidｅ(t+2)’）； ｘ=x1—sｙｍ（"Heａvisidｅ（t-2）’）； y=ｙ＊x; ｅzplｏt(t，y,[—10,10]） grid

信号与系统实验报告总结【第三篇】

信号与系统实验

实验一常用信号的观察

方波：

正弦波：

三角波：

在观测中，虚拟示波器完全充当实际示波器的作用，在工作台上连接AD1为示波器的输入，输入方波、正弦波、三角波信号时，可在电脑上利用软件观测到相应的波形，其纵轴为幅值可通过设置实现幅值自动调节以观测到最佳大小的波形，其横轴为时间，宜可通过设置实现时间自动调节以观测到最佳宽度的波形。 实验四非正弦周期信号的分解与合成

方波DC信号：

DC信号几乎没有，与理论相符合，原信号没有添加偏移。

方波基波信号：

基波信号为与原方波50Hz信号相对应的频率为50Hz的正弦波信号，是方波分解的一次谐波信号。

方波二次谐波信号:

二次谐波信号频率为100Hz为原方波信号频率的两倍，幅值较一次谐波较为减少。

方波三次谐波信号：

三次谐波信号频率为150Hz为原方波信号的三倍。幅值较一二次谐波大为减少。

方波四次谐波信号:

四次谐波信号的频率为200Hz为原方波信号的四倍。幅值较三次谐波再次减小。

方波五次谐波信号：

五次谐波频率为250Hz为原方波信号的五倍。幅值减少到以内，几乎可以忽略。

综上可知：50Hz方波可以分解为DC信号、基波信号、二次、三次、四次、五次谐波信号…，无偏移时即无DC信号，DC信号幅值为0。分解出来的基波信号即一次谐波信号频率与原方波信号频率相同，幅值接近方波信号的幅值。二次谐波、三次谐波、四次谐波、五次谐波依次频率分别为原方波信号的二、三、四、五倍，且幅值依次衰减，直至五次谐波信号时几乎可以忽略。可知，方波信号可分解为多个谐波。 方波基波加三次谐波信号：

基波叠加上三次谐波信号时，幅值与方波信号接近，形状还有一定差异，但已基本可以看出叠加后逼近了方波信号。

方波基波加三次谐波信号加五次谐波信号：

基波信号、三次谐波信号、五次谐波信号叠加以后，比基波信号、三次谐波信号叠加后的波形更加接近方波信号。

综上所述：方波分解出来的各次谐波以及DC信号，叠加起来以后会逼近方波信号，且叠加的信号越多，越是接近方波信号。说明，方波信号可有多个谐波合成。

三角波DC信号：

三角波基波信号：

三角波二次谐波信号：

三角波三次谐波信号：

三角波四次谐波信号：

三角波五次谐波信号：

三角波基波加三次谐波信号：

三角波基波加三次谐波加五次谐波信号：

三角波信号的分析与方波信号的分析基本一致，可以看出三角波也可以分解为多个谐波，并且相应的多个多次谐波可以合成三角波信号，且参与合成的波形越多，合成波越是逼近三角波信号。

综合两个波形来看，可知任何周期性函数均可分解为相应的傅里叶展开式里所包含的直流分量和各次谐波项。且任何周期性函数均可由锁对应的直流分量和各次谐波项所合成，参与合成的信号越多，结果越逼近周期性函数的图形。

实验思考题

1．什么样的周期性函数没有直流分量和余弦项；

答：无偏移的周期性函数没有直流分量，当周期性函数为奇函数时没有直流分量和余弦项。

2．分析理论合成的波形与实验观测到的合成波形之间误差产生的原因。

答：理论合成的波形不能把所有无限个谐波合成起来，故必然产生误差，且实验设备、实验方法也存在一定的误差。

实验二 零输入、零状态级完全响应

零输入响应下降沿采样：

零输入响应上升沿采样：

可见，零输入响应按照指数形式下降，最终降为零。其规律符合-1tU（t）＝RCc2e。

信号与系统实验报告【第四篇】

长江大学电工电子实验中心实验报告

野外获得的重力数据要作进一步处理和解释才能解决所提出的地质任务﹐主要分3个阶段﹕野外观测数据的处理﹐并绘制各种重力异常图﹔重力异常的分解（应用平均法﹑场的变换﹑频率滤波等方法）﹐即从叠加的异常中分出那些用来解决具体地质问题的异常﹔确定异常体的性质﹑形状﹑产状及其他特征参数。解释

解释分为定性的和定量的两个内容﹐定性解释是根据重力图并与地质资料对比﹐初步查明重力异常性质和获得有关异常源的信息。除某些构造外﹐对一般地质体重力异常的解释可遵循以下的一些原则﹕极大的 正异常说明与围岩比较存在剩馀质量﹔反之﹐极小异常是由质量亏损引起的。靠近质量重心﹐在地表投影处将观测到最大异常。最大的水平梯度异常相应于激发体的边界。延伸异常相应于延伸的异常体﹐而等轴异常相应于等轴物体在地表的投影。对称异常曲线说明质量相对于通过极值点的垂直平面是对称分布的﹔反之﹐非对称曲线是由于质量非对称分布引起的。在平面上出现几个极值的复杂异常轮廓﹐表明存在几个非常接近的激发体。定量解释是根据异常场求激发体的产状要素建立重力模型。一种常用的反演方法是选择法﹐即选择重力模型使计算的重力异常与观测重力异常间的偏差小于要求的误差。

由于重力反演存在多解性﹐因此﹐必须依靠研究地区的地质﹑钻井﹑岩石密度和其他物探资料来减少反演的多解性。应用运用领域

在区域地质调查﹑矿产普查和勘探的各个阶段都可应用重力勘探﹐要根据具体的地质任务设计相应的野外工作方法。应用条件

应用重力勘探的条件是﹕被探测的地质体与围岩的密度存在一定的差别﹔被探测的地质体有足够大的体积和有利的埋藏条件﹔干扰水平低。意义 重力勘探解决以下任务﹕

1、研究地壳深部构造﹔研究区域地质构造﹐划分成矿远景区﹔

2、掩盖区的地质填图﹐包括圈定断裂﹑断块构造﹑侵入体等﹔

3、广泛用于普查与勘探可燃性矿床（石油﹑天然气﹑煤）﹐

4、查明区域构造﹐确定基底起伏﹐发现盐丘﹑背斜等局部构造﹔

5、普查与勘探金属矿床（铁﹑铬﹑铜﹑多金属及其他）﹐主要用于查明与成矿有关的构造和岩体﹐进行间接找矿﹔

6、也常用于寻找大的﹑近地表的高密度矿体﹐并计算矿体的储量﹔工程地质调查﹐如探测岩溶﹐追索断裂破碎带等。