本科實驗報告
課程名稱:訊號與系統實驗
學院(系):電信學部
專業:電子
班級:學號:
姓名: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頁...