資源描述:
《實驗一 ATLAB環(huán)境與數(shù)值計算.doc》由會員上傳分享,免費在線閱讀,更多相關內容在行業(yè)資料-天天文庫。
1、實驗一MATLAB環(huán)境與數(shù)值計算一:實驗目的和要求(1)熟練掌握MATLAB的啟動和退出、MATLAB的命令窗口;常用選單和工具欄,通過實例初步認識對MATLAB的命令和程序運行等。(2)熟練掌握MATLAB變量的使用;(3)熟練掌握矩陣的創(chuàng)建、掌握MATLAB的矩陣和數(shù)組的運算;(4)熟悉MATLAB多項式的運用,掌握MATLAB的擬合和插值方法、積分運算和統(tǒng)計運算二:實驗內容:1:用MALTAB命令和M程序兩種方式求解自建的三元一次方程組。(1):MALTAB命令A=[2,-1,3;3,1,-5;4,-1,1]A=2-1331-54
2、-11>>B=[5;5;9]B=559>>X1=ABX1=2-10(2):M程序.創(chuàng)建M程序后運行結果如下:A=[2,-1,3;3,1,-5;4,-1,1]B=[5;5;9]X1=ABA=2-1331-54-11B=559X1=2-102:用MALTAB命令和M程序兩種方式求解自建的3×3復數(shù)矩陣的模和相角。(1):用MALTAB命令>>A=[1+2*i,3-2*i,4-5*i;5+6*i,7+8*i,9-1*i;9+8*i,9+7*i,9+6*i]A=1.0000+2.0000i3.0000-2.0000i4.0000-5.000
3、0i5.0000+6.0000i7.0000+8.0000i9.0000-1.0000i9.0000+8.0000i9.0000+7.0000i9.0000+6.0000i>>norm(A)ans=25.5949>>angle(A)*180/pians=63.4349-33.6901-51.340250.194448.8141-6.340241.633537.875033.6901M文件運行結果:>>typehc2.mA=[1+2*i,3-2*i,4-5*i;5+6*i,7+8*i,9-1*i;9+8*i,9+7*i,9+6*i]nor
4、m(A)angle(A)*180/pi>>hc2.mA=1.0000+2.0000i3.0000-2.0000i4.0000-5.0000i5.0000+6.0000i7.0000+8.0000i9.0000-1.0000i9.0000+8.0000i9.0000+7.0000i9.0000+6.0000ians=25.5949ans=63.4349-33.6901-51.340250.194448.8141-6.340241.633537.875033.69013:自建一個4╳4矩陣,求逆矩陣、產生對角矩陣、求矩陣的秩等。>>A=[3
5、765;4218;6197;9799]A=3765421861979799>>inv(A)ans=-0.2839-0.0766-0.12600.32380.1077-0.0215-0.12920.05980.0973-0.08610.1499-0.09410.10290.17940.0766-0.1651>>diag(A)ans=3299>>rank(A)ans=44:自建兩個4╳4矩陣進行加、減、乘法和左除運算。>>A=magic(4)A=16231351110897612414151>>B=[2-135;31-55;4-119;97
6、99]B=2-13531-554-1199799>>C=A+BC=18161881251313672113212410>>E=A-BE=143082101535853-576-8>>D=A*BD=167741582341555242242171761062421192-34234>>BAans=-474.3333-35.6667-92.3333-304.3333283.000023.000057.0000181.000012.83330.16671.83337.8333241.833319.166747.8333155.83335:自
7、建一個一元四次多項式,進行多項式求值,求根計算。>>A=[19799];>>X1=1:10X1=12345678910>>polyval(A,X1)ans=Columns1through63514342398919793555Columns7through10590392331377919799>>R=roots(A)R=-8.26920.1055+1.0698i0.1055-1.0698i-0.94196:自建一個一元三次多項式,取十個不同自變量值,計算多項式的值,并分別進行一階、二階擬合,圖示原始數(shù)據(jù)和擬合結果.>>typehc1.
8、mx1=1:10A1=[9689];Y1=polyval(A1,x1)A2=polyfit(x1,Y1,1)Y2=polyval(A2,x1)A3=polyfit(x1,Y1,2)Y3=polyval(A3