訊號與系統實驗matlab報告

2021-09-19 12:42:17 字數 4577 閱讀 7629

本科實驗報告

課程名稱:訊號與系統實驗

學院(系):電信學部

專業:電子

班級:學號:

姓名:1. 掌握週期訊號的傅利葉級數展開

2. 掌握週期訊號的有限項傅利葉級數逼近

3. 掌握週期訊號的頻譜分析

4. 掌握連續非週期訊號的傅利葉變換

5. 掌握傅利葉變換的性質

t=-3:0.01:3;

n0=-3;n1=-1;t0=2;

for i=0:2

t1=n0:0.01:n0+t0/2;x1=t1-n0;

t2=n1-t0/2:0.01:n1;x2=-t2+n1;

plot(t1,x1,'r',t2,x2,'r');hold on;

n0=n0+t0;

n1=n1+t0;

endn_max=[1 3 7 15 31];

n=length(n_max);

for k=1:n

n=1; sum=0;

while (n<(n_max(k)+1))

b=4./pi/pi/n/n;

y=b*cos(n*pi*t);

sum=sum+y;

n=n+2;

endfigure;

n0=-3;n1=-1;t0=2;

for i=0:2

t1=n0:0.01:n0+t0/2;x1=t1-n0;

t2=n1-t0/2:0.01:n1;x2=-t2+n1;

plot(t1,x1,'r',t2,x2,'r');hold on;

n0=n0+t0;

n1=n1+t0;

endy=sum+0.5;

plot(t,y,'b');

xlabel('t'),ylabel('wove');

hold off;

axis([-3.01 3.01 -0.01 1.01]);

grid on;

title(['the max=',num2str(n_max(k))])

endfor k=1:3;

n=-30:30;tao=k;t=2*k; w=2*pi/t;

x=n*tao*0.5

fn1=sinc(x/pi); fn=tao*fn1.*fn1;

subplot(3,1,k),stem(n*w,fn);grid on

title(['t=',num2str(2*k)]); axis([-30 30 0 k]);

endft=sym('sin(2*pi*(t-1))/(pi*(t-1))');

fw=fourier(ft);

subplot(2,1,1);

ezplot(abs(fw));

grid on;

title('fudupu');

phase=atan(imag(fw)/real(fw));

subplot(2,1,2);

ezplot(phase);

grid on;

title('xiangweipu');

ft=sym('(sin(pi*t)/(pi*t))^2');

fw=fourier(ft);

subplot(2,1,1);

ezplot(abs(fw));

grid on;

title('fudupu');

phase=atan(imag(fw)/real(fw));

subplot(2,1,2);

ezplot(phase);

grid on;

title('xiangweipu');

(1)syms t

fw=sym('10/(3+i*w)-4/(5+i*w)')

ft=ifourier(fw,t)

ezplot(ft);grid on

(2)syms t

fw=sym('exp(-4*w^2)')

ft=ifourier(fw,t)

ezplot(ft);grid on

dt = 0.01;

t = -0.5:dt:0.5;

ft = uct(t+0.5)-uct(t-0.5);

n = 2000;

k = -n:n;

w = 2*pi*k/((2*n+1)*dt);

f = dt * ft*exp(-j*t'*w);

plot(w,f), grid on

這是第一次訊號上機實驗,在這次實驗中第一次接觸到了matlab這個強大的工程軟體,同時學會了對繪製時訊號的時域波形和對訊號進行頻域分析。非常直觀的看到了,訊號的疊加合成等現象。

1 建立系統的概念

2 掌握連續時間系統的單位衝激響應的求解

3 掌握連續時間系統單位階躍響應的求解

4 掌握連續時間系統零極點的求解

5 分析系統零極點對系統幅頻特性的影響

6 分析零極點對系統穩定性的影響

7 介紹常用訊號處理的matlab工具箱

已知系統的微分方程為,計算該系統的單位衝激響應和單位階躍響應,b=[1,4];

a=[1,3,2];

sys=tf(b,a);

t=0:0.1:10;

y=impulse(sys,t);

subplot(2,1,1);

plot(t,y);

xlabel('±')

ylabel('¤ì')

b=[1];

a=[1,2,1];

sys=tf(b,a);

t=0:0.1:10;

y=step(sys,t);

subplot(2,1,2);

plot(t,y);

xlabel('±')

ylabel('×ì')

1. 實現卷積,其中: p=0.01

nf=0:p:2

f=2*((nf>=0)-(nf>=2

nh=0:p:10

h=exp(-nh

y=conv(f,h);

subplot(3,1,1),stairs(nf,f);

title('f(t)');axis([0 3 0 2.1]);

subplot(3,1,2),stairs(nh,h);

title('h(t)');axis([0 3 0 1.1]);

subplot(3,1,3),plot(y);

title('y(t)=f(t)*h(t)');

2. 已知二階系統方程

對下列情況分別求單位衝激響應,並畫出其波形。

a. b.

c. d.

b=[3];

a=[1,4,3];

sys=tf(b,a);

t=0:0.1:10;

y=impulse(sys,t);

plot(t,y);

xlabel('±')

ylabel('¤ì')

b=[1];

a=[1,2,1];

sys=tf(b,a);

t=0:0.1:10;

y=impulse(sys,t);

plot(t,y);

xlabel('±')

ylabel('¤ì')

b=[1];

a=[1,0,1];

sys=tf(b,a);

t=0:0.1:10;

y=impulse(sys,t);

plot(t,y);

xlabel('±')

ylabel('¤ì')

3. 求下列系統的零極點

(1)(2)

b=[1 0 -4

a=[1 2 -3 2 1

sys=tf(b,a

subplot(3,1,1),pzmap(sys);

subplot(3,1,2),step(b,a

subplot(3,1,3),impulse(b,a);

1)b=[5 20 25 0

a=[1 5 16 30

sys=tf(b,a

subplot(3,1,1),pzmap(sys);

subplot(3,1,2),step(b,a

subplot(3,1,3),impulse(b,a);

4. 對於更多零極點和不同零極點位置的連續系統,作出系統的零極點圖;分析系統是否穩定?若穩定,作出系統的幅頻特性曲線和相頻特性曲線;

(1) 1個極點s=0,增益k=1;

(2) 2個共軛極點,增益k=1;

(3) 零點在s=0.5,極點在,增益k=1

b=[1

a=[1 0

sys=tf(b,a

pzmap(sys);

xi tong wending

b=[1

a=[1 0 5

sys=tf(b,a

pzmap(sys);

lijie wending

b=[1 -0.5

a=[1 -0.2 25.01

sys=tf(b,a

pzmap(sys);

wending

這次實驗學會了用matlab對系統進行一些分析,並求解連續時間系統的單位衝響應和單位階躍響應。對於以後實際問題中對連續時間系統的分析有了很大幫助。

訊號與系統MATLAB實驗報告

一 實驗名稱 matlab對連續訊號與系統的時域分析 頻域分析和s域分析 matlab對離散訊號與系統的時域分析。二 實驗目的 1.學習用matlab描述常用訊號的方法。2.掌握連續時間訊號和離散時間訊號的描述。3.利用matlab計算訊號卷積。4.掌握訊號頻譜的定義,理解非週期訊號頻譜密度的概念。...

訊號與系統實驗報告

實驗報告樣板 王寶忠編寫 適用專業 通訊工程 電子資訊工程 電子資訊科學與技術 其他電類專業 江蘇科技大學電子資訊學院 二零零六年七月 實驗學時 2學時 實驗型別 驗證 實驗要求 必修 一 實驗目的 1 研究動態網路的階躍響應。2 驗證卷積積分。二 實驗內容 1 觀察衝激訊號和階躍訊號的關係。2 觀...

「訊號與系統」實驗報告

實驗報告 課程名稱 訊號與系統實驗 任課教師 實驗名稱 年級 專業 學號姓名 日期 年月日 一 實驗目的 1.2.二 實驗內容 1.2.三 主要演算法與程式 四 實驗結果與分析 五 實驗小結 說明 請注意 在正式報告中把本說明頁刪除。1 標題用小三號宋體加粗,正文用小四號宋體。2 實驗報告不超過5頁...