點(diǎn)擊上方藍(lán)字關(guān)注“公眾號(hào)”

基本信號(hào)在MATLAB中的表示

Part.1

當(dāng)初學(xué)信號(hào)與系統(tǒng)，天書(shū)睡過(guò)去啊，一到考試心中mmp。理論的枯燥，有了MATLAB就生動(dòng)形象啦。

本節(jié)掌握：

1．?學(xué)會(huì)用MATLAB表示常用連續(xù)信號(hào)的方法；

Part.2

連續(xù)信號(hào)的MATLAB表示

MATLAB提供了大量的生成基本信號(hào)的函數(shù)，例如指數(shù)信號(hào)、正余弦信號(hào)。

表示連續(xù)時(shí)間信號(hào)有兩種方法，一是數(shù)值法，二是符號(hào)法。

數(shù)值法是定義某一時(shí)間范圍和取樣時(shí)間間隔，然后調(diào)用該函數(shù)計(jì)算這些點(diǎn)的函數(shù)值，得到兩組數(shù)值矢量，可用繪圖語(yǔ)句畫(huà)出其波形。

符號(hào)法是利用MATLAB的符號(hào)運(yùn)算功能，需定義符號(hào)變量和符號(hào)函數(shù)，運(yùn)算結(jié)果是符號(hào)表達(dá)的解析式，也可用繪圖語(yǔ)句畫(huà)出其波形圖。

指數(shù)信號(hào)

指數(shù)信號(hào)在MATLAB中用exp函數(shù)表示。

A=1; a=-0.4;
t=0:0.01:10;?????? %定義時(shí)間點(diǎn)
ft=A*exp(a*t);???? %計(jì)算這些點(diǎn)的函數(shù)值
plot(t,ft);????????? %畫(huà)圖命令，用直線段連接函數(shù)值表示曲線
grid on;?????????? %在圖上畫(huà)方格

正弦信號(hào)

正弦信號(hào)在MATLAB中用 sin 函數(shù)表示。

調(diào)用格式為 ?ft=A*sin(w*t+phi)

A=1; w=2*pi; phi=pi/6;

t=0:0.01:8; ???????????????????%定義時(shí)間點(diǎn)

ft=A*sin(w*t+phi); ????????????%計(jì)算這些點(diǎn)的函數(shù)值

plot(t,ft); ?????????????????????%畫(huà)圖命令

grid on; ??????????????????????%在圖上畫(huà)方格

異抽樣信號(hào)

抽樣信號(hào)Sa(t)=sin(t)/t在MATLAB中用 sinc 函數(shù)表示。

t=-3*pi:pi/100:3*pi;

?????ft=sinc(t/pi);

?????plot(t,ft);

?????grid on;

?????axis([-10,10,-0.5,1.2]); ???%定義畫(huà)圖范圍，橫軸，縱軸

?????title('抽樣信號(hào)') ????????%定義圖的標(biāo)題名字

三角信號(hào)

三角信號(hào)在MATLAB中用 tripuls 函數(shù)表示。

調(diào)用格式為 ft=tripuls(t,width,skew)，產(chǎn)生幅度為1，寬度為width，且以0為中心左右各展開(kāi)width/2大小，斜度為skew的三角波。

width的默認(rèn)值是1，skew的取值范圍是-1~+1之間。一般最大幅度1出現(xiàn)在t=(width/2)*skew的橫坐標(biāo)位置。

t=-3:0.01:3;

??ft=tripuls(t,4,0.5);

?????plot(t,ft); ?grid on;

?????axis([-3,3,-0.5,1.5]);

復(fù)指數(shù)信號(hào)

調(diào)用格式是f=exp((j*w)*t)

?????? t=0:0.01:15;

w=pi/4;

X=exp(j*w*t);

Xr=real(X); ???????%取實(shí)部 ??

Xi=imag(X); ??????%取虛部

Xa=abs(X); ???????%取模

Xn=angle(X); ?????%取相位

subplot(2,2,1),plot(t,Xr),axis([0,15,-(max(Xa)+0.5),max(Xa)+0.5]),

title('實(shí)部');

subplot(2,2,3),plot(t,Xi),axis([0,15,-(max(Xa)+0.5),max(Xa)+0.5]),

title('虛部');

subplot(2,2,2), plot(t,Xa),axis([0,15,0,max(Xa)+1]),title('模');

subplot(2,2,4),plot(t,Xn),axis([0,15,-(max(Xn)+1),max(Xn)+1]),title('相角');

?????%subplot(m,n,i) 命令是建立m行n列畫(huà)圖窗口，并指定畫(huà)圖位置i

矩形脈沖信號(hào)

調(diào)用格式是f=exp((a+j*b)*t)

單位階躍信號(hào)

單位階躍信號(hào)u(t)用“t>=0”產(chǎn)生，調(diào)用格式為ft=(t>=0)

正弦信號(hào)符號(hào)算法

syms t ????????????????????????%定義符號(hào)變量t

y=sin(pi/4*t) ???????????????????%符號(hào)函數(shù)表達(dá)式

ezplot(y,[-16,16]) ???????????????%符號(hào)函數(shù)畫(huà)圖命令

或者

f=sym('sin(pi/4*t)') ??????????????%定義符號(hào)函數(shù)表達(dá)式

ezplot(f,[-16,16])

單位階躍信號(hào)

MATTLAB符號(hào)數(shù)學(xué)函數(shù)Heaviside表示階躍信號(hào)，但要畫(huà)圖需在工作目錄創(chuàng)建Heaviside的M文件

function f=Heaviside(t)

f=(t>0); ?

保存，文件名是Heaviside ，調(diào)用該函數(shù)即可畫(huà)圖，例

t=-1:0.01:3;

f=heaviside(t);

plot(t,f)

axis([-1,3,-0.2,1.2]) ?

或者

y=sym('Heaviside(t)');

ezplot(y,[-1,5]);grid on ??

MATLAB各種版本，關(guān)注公眾號(hào)回復(fù)相應(yīng)版本即可。

掃碼關(guān)注

不迷路