MATLAB上機(jī)作業(yè)

專業(yè)班級(jí)：

學(xué)生姓名：

學(xué)生學(xué)號(hào)：

年 月

作業(yè)1

1.用MATLAB可以識(shí)別的格式輸入下面兩個(gè)矩陣

,

再求出它們的乘積矩陣,并將矩陣的右下角子矩陣賦給矩陣。賦值完成后,調(diào)用相應(yīng)的命令查看MATLAB工作空間的占有情況。A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]

B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;]

B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;]

C=A*B

D=C(4:5,4:6);

whos;

2.設(shè)矩陣,求,,,,,并求矩陣的特征值和特征向量。 A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 2;]

det(A)

inv(A)

A.^3

2*A+inv(A)

3*A-A'

[V,D]=eig(A)

abs (A)

3.解下列矩陣方程:

A=[0 1 0;1 0 0;0 0 1;];

B=[1 0 0;0 0 1;0 1 0;];

C=[1 -4 3; 2 0 -1;1 -2 0;];

X=inv(A)*C*inv(B)

4.一球從100米高度自由落下,每次落地后反跳回原高度的一半,再落下。求它在第10次落地時(shí),共經(jīng)過多少米?第10次反彈有多高?syms x y sum ;

x=100;

y=0.5*x;

sum=0;

for i=2:10,

sum=sum+x+y;

x=0.5*x;

y=0.5*x;

end

sum

y

5.用MATLAB語(yǔ)言實(shí)現(xiàn)下面的分段函數(shù)

function y=f(x)

if x>10

y=5;

else

if x>=-10&x<=10

y=1/2*x;

else

y=-5;

end

end

6.分別用for和while循環(huán)編寫程序,求出

并考慮一種避免循環(huán)的簡(jiǎn)潔方法來(lái)進(jìn)行求和,并比較各種算法的運(yùn)行時(shí)間。syms x K ;

x=2;

K=0;

for i=1:63,

K=K+x^i;

end

K

7.應(yīng)用MATLAB語(yǔ)言及二分法編寫求解一元方程在區(qū)間[3，6]的實(shí)數(shù)解的算法，要求絕對(duì)誤差不超過0.001。format long

A=[1 -14 59 -70];

a=3;

b=6;

c=0.001;

while c<0.5*(b-a), x1=(a+b)/2;

f1=polyval(A,x1);

fa=polyval(A,a);

fb=polyval(A,b);

if f1==0; x=x1;

else if f1*fa>0;

a=x1;

else f1*fb>0;

b=x1;

end

end

x=x1

8.二階系統(tǒng)的單位階躍響應(yīng)為，在同一平面繪制分別為0，0.3，0.5，0.707的單位階躍響應(yīng)曲線。要求：

(1)

四條曲線的顏色分別為藍(lán)、綠、紅、黃，線型分別為“——”、“……”、“oooooo”、“++++++”；

(2)添加橫坐標(biāo)軸和縱坐標(biāo)軸名分別為“時(shí)間t”和“響應(yīng)y”，并在平面圖上添加標(biāo)題“二階系統(tǒng)曲線”和網(wǎng)格；

(3)在右上角添加圖例(即用對(duì)應(yīng)的字符串區(qū)分圖形上的線)，并分別在對(duì)應(yīng)的響應(yīng)曲線的第一個(gè)峰值處標(biāo)示“zeta＝0”、“zeta＝0.3”、“zeta＝0.5”、“zeta＝0.707”。t=[0:0.1:10];

kos=0;

y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));

plot(t,y,'b-')

hold on;

kos=0.3;

y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));

plot(t,y,'g.')

hold on;

kos=0.5;

y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));

plot(t,y,'ro')

hold on;

kos=0.707;

y=1-1/sqrt(1-kos^2)*exp(-kos*t).*sin(sqrt(1-kos^2)*t+2*cos(kos));

plot(t,y,'y+')

hold on;

xlabel('時(shí)間t');

ylabel('響應(yīng)y');

grid on;

title(

總結(jié)

以上是生活随笔為你收集整理的表上作业法matlab解决,MATLAB上机作业.doc的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。