數(shù)字信號(hào)處理實(shí)驗(yàn)教案

（1112341用） 數(shù)字信號(hào)解決實(shí)驗(yàn)規(guī) 預(yù)備實(shí)驗(yàn)用于信號(hào)解決的MATLAB基本操 實(shí)驗(yàn) 系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定 實(shí)驗(yàn) 時(shí)域采樣與頻域采 實(shí)驗(yàn) 用FFT對(duì)信號(hào)作頻譜分 實(shí)驗(yàn) IIR數(shù)字濾波器設(shè)計(jì)及軟件實(shí) 實(shí)驗(yàn) FIR數(shù)字濾波器設(shè)計(jì)與軟件實(shí) 附錄1：MATLAB編程及繪圖辦 附錄2：MATLAB矩陣及矩陣的運(yùn) 實(shí)驗(yàn)基本規(guī)定（以軟件設(shè)計(jì)為重要體現(xiàn)形式預(yù)備實(shí)驗(yàn)MATLAB（3實(shí)驗(yàn) 系統(tǒng)響應(yīng)及系統(tǒng)穩(wěn)定性（3學(xué)時(shí)實(shí)驗(yàn) 時(shí)域采樣與頻域采樣（3學(xué)時(shí)實(shí)驗(yàn) 用FFT對(duì)信號(hào)作頻譜分析（3學(xué)時(shí)實(shí)驗(yàn) IIR數(shù)字濾波器設(shè)計(jì)及軟件實(shí)現(xiàn)（3學(xué)時(shí)實(shí)驗(yàn) FIR數(shù)字濾波器設(shè)計(jì)與軟件實(shí)現(xiàn)（3學(xué)時(shí)姓名 學(xué)號(hào) 實(shí)驗(yàn)日期 [1]高西全，丁玉美.數(shù)字信號(hào)解決（第三版）西電出版社張志勇等.MATLABRa.北京航空航天大學(xué)出版社.MATLAB[3]程佩清.數(shù)字信號(hào)解決教程[M].清華大學(xué)出版社.預(yù)備實(shí)驗(yàn)用于信號(hào)解決的MATLABMATLAB②MATLAB①熟悉MATLAB運(yùn)用冒號(hào)（：）X1=[1234X2=[1.0001.5002.000X3=[5432C）分別輸入矩陣123456789D）分別計(jì)算A+B，A-B，A+3，A-4，A*B，A.*B，C=inv（A）,A/B,A./BE）分別計(jì)算sin(x1),cos(x1),exp(x1),lg(x2),sqrt(x2)運(yùn)用冒號(hào)（：）X=j:kX=[j,j+1,j+2,…,k-I>0j<kI<0j>k,則生成向量I<0j>k,則生成向量例X1=[1234X2=[1.0001.5002.0002.500X3=5：-432 0B= n*n0B= m*n0B=Zeros(Size(A)) A0注：m,n二矩陣的算術(shù)運(yùn)算A-對(duì)于矩陣和標(biāo)量（一種數(shù)）A-

Sin(A),cos(A),exp(A),log(A)( 式中A能夠是標(biāo)量也能夠是矩陣練習(xí)xlabel('0-四．繪圖（二維繪圖1plot(y);%繪持續(xù)曲線圖。stem(y);%繪離散序列波形yyy中元2plot(x,y);繪持續(xù)曲線圖。stem(x,y);%繪離散序列波形x,yxy元素的值為縱坐標(biāo)，一一MATLAB會(huì)自動(dòng)設(shè)定所畫曲線的顏色和線型?？勺詣?dòng)設(shè)立，在輸入變量組后，加一叉標(biāo)，-.點(diǎn)劃線，下面給出窗口分割語句x，yx*ysubplot(3,2,2)3*21MATLAB2MATLAB1：MATLABMATLABMATLABC語言：以字母開頭，由字母、下31個(gè)。MATLAB慣用運(yùn)算符（1＋—矩陣右除(矩陣元素右除矩陣左除（矩陣元素左除矩陣乘（矩陣元素乘∧（.矩陣求冪（矩陣元素求冪[%‘><&|~MATLABlog(x)、log2(x0、e（2、10）[1,2,…]=函數(shù)名（1MATLAB3n1*n2title(‘字符串xlable(‘字符串xxlable(‘字符串y1yyy中元2x,yxy元素的值為縱坐標(biāo)，一一MATLAB>>Y=[0.7,0.5,0.2,->>Y=[0.7,0.5,0.2,- %畫連線圖, %畫散點(diǎn)圖,圖1.1持續(xù)波形 圖1.2散點(diǎn)波形;..]holdon2holdholdon3plot 4plotyy命令，設(shè)有兩個(gè)縱坐標(biāo)1：figure（n）：MATLABplotplotfigure命令打開新的圖形窗口。格式2：clf：去除現(xiàn)在圖形窗的內(nèi)容3：close:4：closeall:subplot(2,2,4);fill(t,y,’r’);subplot(2,2,4);fill(t,y,’r’);MATLAB與大多數(shù)計(jì)算機(jī)語言同樣，MATLABMATLAB的次序構(gòu)造就是由體現(xiàn)式構(gòu)成的語句，語句之間用逗號(hào)或分號(hào)隔開，這些MATLABforwhilefor-endforfor循環(huán)變量＝1212313之間，如果條件滿（DFT

jNX(k)x(n)Wnk NN

k

,N fork=0:N-1forn=0:N-while-endwhilewhile3if-else-endififswitchswitchcase體現(xiàn)式的值＝1為真時(shí)的語句組case值2體現(xiàn)式的值＝23forM文獻(xiàn)的運(yùn)行時(shí)間,>>elapsed_time>>elapsed_time=00.38forfork=1:100在命令窗口運(yùn)行,>>elapsed_time >>elapsed_time

2：MATLAB矩陣須用中括號(hào)[]>>x=[12345678];%一維1×8>>x=0:0.02:1%以起始值=0、增量值=0.02、終止值=1的矩陣(陣列)%linspace函數(shù)，起始值=0,終止值=1,以等間隔（1-0）/100100A（i,j）Aij例1．A>>>>>>A162738495MATLABA（u,v）A的子矩陣，u,vA([1,3],:)%A1ansA(:,[2,4,5]A2,4,5ans A（L1，：）表達(dá)提取A的L1指定的行、全部列；A（：，L2A([1,3],:)%A1ansA(:,[2,4,5]A2,4,5ans 123457930.5>>rand(>>rand(‘seed’,0);>>x=rand(1,10);%1X10>>xColumns1through Columns7throughx4A2,4,6>>>>>>A(:,[2,4,6])=[]A= 5>>A 35>>A 35079000000使用，這時(shí)使用*.mat>>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>load save途徑\文獻(xiàn)名 >>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>loadMATLABMATLAB命令窗口定義的變量、運(yùn)算成果和有關(guān)數(shù)據(jù)，此內(nèi)存空間就是workspaceworkspace*.matFilesaveworkspaceas….項(xiàng)即可矩陣 矩陣 a-ABAB矩陣的除法從逆矩陣的概念引申，