資源描述：
《matlab控制與系統(tǒng)仿真 實驗四.doc》由會員上傳分享，免費在線閱讀，更多相關內容在應用文檔-天天文庫。
1、黃淮學院電子科學與工程系matlab控制與系統(tǒng)仿真課程驗證性實驗報告實驗名稱MATLAB的數(shù)值計算實驗時間2013年05月21日學生姓名王茂勝實驗地點同組人員無專業(yè)班級電技1001B1、實驗目的1．掌握積分的數(shù)值計算方法；2．掌握微分的數(shù)值計算方法。2、實驗主要儀器設備和材料：計算機PC一臺、MATLAB軟件7.03、實驗內容和原理：實驗內容：（1）求定積分：1）2）（2）求二重定積分1)2)3)設X由區(qū)間內均勻分布的10個點組成，求向量的1～3階差分。4）設用3種不同的方法求函數(shù)的數(shù)值導數(shù)，并在同一個坐標系中作出的圖像。思考與練
2、習：（1）簡述數(shù)值微分與積分的計算過程。（2）試用函數(shù)quad、quadl和trapz、求積分，比較3種算法的精度（）（3）求定積分1）2）（4）求3重積分：（5）設，用不同的方法求函數(shù)的數(shù)值導數(shù)，并在同一個坐標系中作出的圖像。1）用一個5次多項式擬合函數(shù)，并對求一般意義下的導數(shù)，求出在假設點的值。2）用diff函數(shù)直接求在假設點的數(shù)值導數(shù)。3）先求出導函數(shù)，然后直接求在假設點的導數(shù)。4、實驗方法、步驟：（1）、求定積分：1）、新建M文件m4_1_1.m，輸入內容：symsrx;r=int('sin(x)/x',x,0,2)運行，
3、觀察結果。2）、新建M文件m4_1_2.m，輸入內容：symsx;int('((x-0.3)^2+0.01)^(-1)-((x-0.09)^2+0.04)^(-1)',0,1)運行，觀察結果。（2）、求二重定積分：1）、新建M文件m4_2_1.m，輸入內容：f=@(x,y)exp(-(x.^2+y.^2));I=dblquad(f,0,1,0,1)運行，觀察結果。2）、新建M文件m4_2_2.m，輸入內容：f=@(x,y)abs(cos(x+y));I=dblquad(f,0,pi,0,pi)運行，觀察結果。（3）、新建M文件m4
4、_3.m，輸入內容：x=linspace(0,2*pi,10);y=sin(x);dy1=diff(y)dy2=diff(y,2)dy3=diff(y,3)運行，觀察結果。（4）、新建M文件m4_4.m，輸入內容：x=[-pi:0.1*pi:pi]f=sin(x)./(x+cos(2*x))plot(x,f,'k-')holdonf1=diff(f)./diff(x)plot(x(1:end-1),f1,'r--')運行，觀察結果。5、實驗現(xiàn)象、實驗數(shù)據(jù)記錄：（1）、1）、運行結果：2）、運行結果：（2）、1）、運行結果：2）、運
5、行結果：（3）、運行結果：（3）、運行結果：6、實驗現(xiàn)象、實驗數(shù)據(jù)的分析：分析：1、symsrx;r=int('sin(x)/x',x,0,2)。使用syms定義式中的符號變量，然后直接調用int(f(x),a,b)函數(shù)求解表達式的積分。其中f(x)為所求積分的函數(shù)表達式，a，b分別為積分的上限和下限。2、f=@(x,y)exp(-(x.^2+y.^2));I=dblquad(f,0,1,0,1)。使用兩次一重積分求二重積分，f=@(x,y)exp(-(x.^2+y.^2))，y從0積分到1，x從0積分到1。3、x=linspac
6、e(0,2*pi,10)為定義x的取值范圍，調用差分函數(shù)diff()求差分。其中dy1=diff(y)為計算一階差分，dy2=diff(y,2)為計算二階差分，dy3=diff(y,3)為計算三階差分。4、，利用syms定義符號，利用diff函數(shù)求符號導數(shù)：Y=diff(X,n)，求函數(shù)X的n階導數(shù)，然后再用plot（）函數(shù)繪制函數(shù)的圖形。7、實驗結論：通過本次練習，掌握積分的數(shù)值計算方法和微分的數(shù)值計算方法，達到了預期的實驗目的。指導教師評語和成績評定：實驗報告成績：指導教師簽字：年月日
