資源描述:
《《matlab函數(shù)學(xué)習(xí)》PPT課件》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、數(shù)字信號(hào)處理與MATLABhttp://www2.nbu.edu.cn/jiangzhidi/mediasystemhttp://bighand.zxq.net/mediasystem主講:蔣志迪寧波大學(xué)信息科學(xué)與工程學(xué)院辦公室:曹光彪信息樓424電話:665929email:bighandjzd@126.com多媒體信息處理系統(tǒng)設(shè)計(jì)Matlab中的信號(hào)表示常見連續(xù)信號(hào)離散信號(hào)離散卷積LTI系統(tǒng)常見連續(xù)信號(hào)指數(shù)信號(hào)Aeaty=A*exp(a*t);指數(shù)序列ak冪運(yùn)算a.^k實(shí)現(xiàn)正弦型信號(hào)??????????
2、????????內(nèi)部函數(shù)cos()和sin()抽樣函數(shù)Sa(t)?????????????sinc(t)矩形脈沖信號(hào)y=rectpuls(t,width)三角波脈沖信號(hào)y=tripuls(t,width,skew)周期信號(hào):正弦信號(hào),周期方波ex1.產(chǎn)生一個(gè)幅度為2,頻率為4Hz,相位為的正弦信號(hào)A=2;f=4;phi=pi/6;w0=2*pi*f;t=0:0.01:1;x=A*sin(w0*t+phi);plot(t,x);6/p非周期信號(hào):指數(shù)信號(hào),階躍信號(hào),取樣函數(shù)Sa(x)等ex2.產(chǎn)生一個(gè)高度為1,
3、寬度為3,延時(shí)為2秒的矩形脈沖信號(hào)t=0:0.02:6;plot(t,rectpuls(t-2-1.5,3));axis([0,6,0,1.5]);注:rectpuls(t-延時(shí)-width/2,width));ex3取樣函數(shù):用sinc(x)命令N=1000;t=-10:20/N:10;x=sinc(t/pi);plot(t,x);gridon隨機(jī)信號(hào):rand產(chǎn)生均勻分布的白噪聲,randn產(chǎn)生高斯分布的白噪聲t=0:0.01:1;y=randn(1,length(t));plot(t,y);gridon
4、;離散信號(hào)周期序列ex4產(chǎn)生一個(gè)幅度為2,頻率為4Hz,相位為的正弦信號(hào)A=2;phi=pi/6;omega=2*pi/12;n=-10:10;x=A*sin(omega*n+phi);stem(n,x,'fill');gridon;單位脈沖序列和單位階躍序列u(n)可以用ones(1,n)和zeros(1,n)來生成單位脈沖序列和單位階躍序列ones(1,n)命令產(chǎn)生1行n列的1值z(mì)eros(1,n)命令產(chǎn)生1行n列的0值ex5產(chǎn)生單位脈沖序列和單位階躍序列u(n-3)n=-2:30;x=[zeros(1,
5、5),1,zeros(1,27)];y=[zeros(1,5),ones(1,28)];subplot(2,1,1);stem(n,x,'fill');gridon;subplot(2,1,2)stem(n,y,'fill');gridon;信號(hào)基本運(yùn)算尺度變換、翻轉(zhuǎn)、時(shí)移、相加、相乘差分與求和、微分與積分離散卷積LTI(LinearTime-Invariant)系統(tǒng)尺度變換、翻轉(zhuǎn)、時(shí)移、相加、相乘t=-3:0.001:3;ft1=tripuls(2*t,4,0.5);subplot(2,1,1)plot(t
6、,ft1)title('f(2t)')ft2=tripuls((2-2*t),4,0.5);subplot(2,1,2)plot(t,ft2)title('f(2-2t)')差分與求和、微分與積分例:已知三角波f(t),畫出其微分與積分的波形%differentiation%integrationh=0.001;t=-3:h:3;t=-3:0.1:3;y1=diff(f2_2(t))*1/h;forx=1:length(t)plot(t(1:length(t)-1),y1)y2(x)=quad('f2_2',
7、-3,t(x));endplot(t,y2)差分y=diff(f);求和y=sum(f(k1:k2));微分y=diff(f)/h;h為數(shù)值計(jì)算所取時(shí)間間隔定積分quad(‘function_name’,a,b);%function_name為被積函數(shù)名,a和b指定積分區(qū)間。其中函數(shù)f2_2(t)為:functionft2=f2_2(t)ft2=tripuls(t,4,0.5);離散卷積matlab中卷積運(yùn)算的指令是c=conv(a,b),其中a和b是有限長(zhǎng)的序列。ex6:計(jì)算下面兩個(gè)序列的卷積a=[-201
8、–13];b=[120-1];兩個(gè)序列的都從0開始的a=[-201-13];b=[120-1];c=conv(a,b);M=length(c)-1;n=0:1:M;stem(n,c,'fill');gridon;xlabel('Timeindexn');系統(tǒng)的時(shí)域分析連續(xù)時(shí)間系統(tǒng)零狀態(tài)響應(yīng)的求解y=lsim(sys,f,t);t:表示計(jì)算系統(tǒng)響應(yīng)的抽樣點(diǎn)向量f:系統(tǒng)輸入信號(hào)向量sys:LTI系統(tǒng)模型