資源描述:
《matlab平時(shí)作業(yè)共享.doc》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、第一章習(xí)題13.請(qǐng)指出如下5個(gè)變量名中,哪些是合法的?abcd-2,xyz_3,ABCDefghabcd-2xyz_33chana變量ABCDefgh6.指令clear,clf,clc各有什么用處?Clear為清除工作空間變量指令,包括內(nèi)存中的變量和數(shù)據(jù);Clc為清命令窗口指令,而保留工作空間內(nèi)容;Clf為清除當(dāng)前圖形指令。7.以下兩種說法對(duì)嗎?(1)“MATLAB的數(shù)值表達(dá)精度與其指令窗中的數(shù)據(jù)顯示精度相同?!保?)MATLAB指令窗中顯示的數(shù)值有效位數(shù)不超過7位?!苯猓海?)錯(cuò)誤;(2)錯(cuò)誤,該顯示數(shù)值有效位可以在菜單中設(shè)置或通過vpa(X,n)進(jìn)行設(shè)置。8.
2、想要在MATLAB中產(chǎn)生二維數(shù)組,下面哪些指令能實(shí)現(xiàn)目的?S=[1,2,3;4,5,6;7,8;9]S=[123;456;789]S=[1,2,3;4,5,6;7,8,9]解:S=[123;456;789];S=[1,2,3;4,5,6;7,8,9]。第三章習(xí)題31.要求在閉區(qū)間上產(chǎn)生具有10個(gè)等距采樣點(diǎn)的一維數(shù)組。試用兩種不同的指令實(shí)現(xiàn)。解:a=linspace(0,2*pi,10)a=0:2*pi/9:2*pi2.由指令rand('state',0),A=rand(3,5)生成二維數(shù)組A,試求該數(shù)組中所有大于0.5的元素的位置,分別求出它們的“全下標(biāo)”和“單下標(biāo)”
3、。解:rand('state',0),A=rand(3,5);L(A>0.5)=NAN程序如下:%注釋:最終求出的該數(shù)組中所有大于0.5的元素的“全下標(biāo)”信息存儲(chǔ)于C中%%%%“單下標(biāo)”信息存儲(chǔ)于D中clear;ii=0;jj=0;x=1;C=zeros(15,2);D=zeros(1,15);rand('state',0),A=rand(3,5);forii=1:3forjj=1:5ifA(ii,jj)>0.5C(x,:)=[ii,jj];D(x)=(ii-1)*5+jj;x=x+1;endendendC,D結(jié)果如下:C=1115222425313233340000
4、00000000D=157910111213140000005.在時(shí)間區(qū)間[0,10]中,繪制曲線。要求分別采取“標(biāo)量循環(huán)運(yùn)算法”和“數(shù)組運(yùn)算法”編寫兩段程序繪圖。解:(1)標(biāo)量循環(huán)運(yùn)算法plot(t,y)fort=0:0.01:10y=1-exp(-0.5*t)*cos(2*t);plot(t,y);grid;holdon;end如圖:(2)數(shù)組運(yùn)算法t=0:0.01:10;y=1-exp(-0.5*t).*cos(2*t);plot(t,y)如圖:6.先運(yùn)行clear,formatlong,rand('state',1),A=rand(3,3),然后根據(jù)A寫出兩個(gè)
5、矩陣:一個(gè)對(duì)角陣B,其相應(yīng)元素由A的對(duì)角元素構(gòu)成;另一個(gè)矩陣C,其對(duì)角元素全為0,而其余元素與對(duì)應(yīng)的A陣元素相同。A={0.952782149656620.598158524172190.836819600676340.704062166775000.840743198113070.518703059724920.953877473592230.442818842235130.02220977857260}B={0.952782149656620000.840743198113070000.02220977857260}C={00.598158524172190.83
6、6819600676340.7040621667750000.518703059724920.953877473592230.442818842235130}第四章習(xí)題43.求函數(shù)的數(shù)值積分,并請(qǐng)采用符號(hào)計(jì)算嘗試復(fù)算。解:(1)數(shù)值積分:clearsymsxx=0:0.1:pi;y=exp(sin(x).^3);>>trapz(x,y)ans=5.0954(2)符號(hào)計(jì)算:symsxIsym=vpa(int(exp(sin(x)^3),x,0,pi))Warning:Explicitintegralcouldnotbefound.>Insym.intat58Isym=5
7、.13701355671761418934690897508764.用quad求取的數(shù)值積分,并保證積分的絕對(duì)精度為。解:Y=quad('exp(-abs(x)).*abs(sin(x))',-5*pi,1.7*pi,1e-9)Y=1.087849437547795.求函數(shù)在區(qū)間中的最小值點(diǎn)。解:clearsymsxy=((sin(5*x))^2)*(exp(0.06*x^2))-1.5*x*cos(2*x)+1.8*abs(x+0.5);yd=diff(y,x);a=solve(yd)yd_a=vpa(subs(yd,x,xs0),6)y_a=vp