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

Sin(A),cos(A),exp(A),log(A)( 式中A能夠是標量也能夠是矩陣練習xlabel('0-四．繪圖（二維繪圖1plot(y);%繪持續(xù)曲線圖。stem(y);%繪離散序列波形yyy中元2plot(x,y);繪持續(xù)曲線圖。stem(x,y);%繪離散序列波形x,yxy元素的值為縱坐標，一一MATLAB會自動設定所畫曲線的顏色和線型。可自動設立，在輸入變量組后，加一叉標，-.點劃線，下面給出窗口分割語句x，yx*ysubplot(3,2,2)3*21MATLAB2MATLAB1：MATLABMATLABMATLABC語言：以字母開頭，由字母、下31個。MATLAB慣用運算符（1＋—矩陣右除(矩陣元素右除矩陣左除（矩陣元素左除矩陣乘（矩陣元素乘∧（.矩陣求冪（矩陣元素求冪[%‘><&|~MATLABlog(x)、log2(x0、e（2、10）[1,2,…]=函數(shù)名（1MATLAB3n1*n2title(‘字符串xlable(‘字符串xxlable(‘字符串y1yyy中元2x,yxy元素的值為縱坐標，一一MATLAB>>Y=[0.7,0.5,0.2,->>Y=[0.7,0.5,0.2,- %畫連線圖, %畫散點圖,圖1.1持續(xù)波形 圖1.2散點波形;..]holdon2holdholdon3plot 4plotyy命令，設有兩個縱坐標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ù)計算機語言同樣，MATLABMATLAB的次序構造就是由體現(xiàn)式構成的語句，語句之間用逗號或分號隔開，這些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為真時的語句組case值2體現(xiàn)式的值＝23forM文獻的運行時間,>>elapsed_time>>elapsed_time=00.38forfork=1:100在命令窗口運行,>>elapsed_time >>elapsed_time

2：MATLAB矩陣須用中括號[]>>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，：）表達提取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使用，這時使用*.mat>>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>load save途徑\文獻名 >>n=1:;x=cos(0.2*pi*n);>>savec:\temp.matx>>loadMATLABMATLAB命令窗口定義的變量、運算成果和有關數(shù)據(jù)，此內(nèi)存空間就是workspaceworkspace*.matFilesaveworkspaceas….項即可矩陣 矩陣 a-ABAB矩陣的除法從逆矩陣的概念引申，