西安科技大學(xué)

MATLAB程序設(shè)計

專業(yè)：信息與計算科學(xué)

班級：1001班

學(xué)號：1008060129

姓名：劉仲能

2012年6月27日

實驗一

2.已

知：CommandWindQW

12344131

?A=[12,34,-4;34,7*87,3,6*7]

A34787，B203

3657327C=AF6*B

求下列表達(dá)式的值：

(1)人+6*8和人6+1(其中I為單位矩陣)

H=A,:3

(2)A*B和A.*BIM

J=&A

(3)AS和A.*3K二[XB]

LF[l([Lf9]f

⑷A/B及B\A

(5)[A,B]和:);BA2]

8=

1234713-11852-10

34F8720346r105

36575349

3-2721

12□2-26844C2

121024

33S85309-725S5

680261

1581

154-52419-130<19

H=

37226233824486041?2839304-G-4

24'37014918860076639300343658503

7868845414211882027貯5343

J=

16*4000-13.0000F.SO00109.400D-13K20DU322.8000

35.8000-76.200050.2000-53.000085.0000-17L0000

67.0000-134.0000S3.-61.60003S.S000-136.2000

1234-4

3657

12-41?-1451

34873

12011019

373-27

20-540

123453016

6789101769

3.設(shè)有矩陣A和BA1112131415，B0234

1617181920970

212223242541311

(1)求它們的乘積C

(2)將矩陣C的右下角3X2子矩陣賦給D

⑶查看MATLAB工作空間的使用情況

CommandWindaw

JiFtl,2,3\5,5,20:21,22,23,24,25]

15.11,12,\3(14,ISJSj17,15,

B二[覓①IS：r-6,3;0,23,-4;9?7.0;4,411]

G二

D二C(end-2:end,end-1:end)

A二B

]2J453016

&7891017-S9

1]1213H15023-4

16\?181920970

212223242b41311

(1)⑵

D-

150U

25fi335237

42339770555?

訊0717

5gg70555?

753S9D717?

12345

二

1I235

767Sg10

115

0II121114

416171也20

5212222425

gArrayEditor-B亡ArrayEditor-G

BI,t.st

=E(XIfelffi]3]*切Stw

123123

1[30161的15077

217弋922535237

3

3D23■442352a397

44

9705670567

54131157533S717

ArrayEditor-A

虢^電馨詹匹"托st堿:日曲eVI

ArrayBJitor-□

a

?EIJtIS匡■

12

iirE2H337

2I705557

3|890717

4.完成下列操作

(1)求［100,999之間能被21整除的數(shù)的個數(shù)。

(2)建立一個字符串向量，刪除其中的大寫字母

(1)(2)

QommandWindowCommandWindow

?cei_irLt=0;?ch^sldjADAkdALjrGETdk1.

for1=100;D90Eubch=ch(firLd(ch>?a"7)0

if10]"(l,21)==0

couivt=count:+l;subch

end

end

sldjkdijfdk

coimt

coun,t=

43

實驗二

3建立一個5X5矩陣，求它的行列式值、跡、秩和范數(shù)

運行截圖：

*CmimndWindow?A=randC6)ardet(A)r=rank(A)t-trace(A)il=iicrni(A"1■)a2=luorni(A)ainf-norm

(A,iiLf)

0.95010.76210.C1540.4057

0.23110.45550.79190.93550-3529

0.60680.01850.92130.8132

0.4S600.82140.73820.41030-0D93

0.89130.444?0.17630.89300.

t=

-0.00712.

范數(shù)如下：

al=

a2二ainf

3.5620

2.8868L2772

29618

4.已知A20512

885

求A的特征值及特征向量，并分析其數(shù)學(xué)意義

運行截圖：

A矩陣的行列式值、跡、秩分別如下：

1/21/31/4為0.95

5.下面是一個線性方程組：1/31/41/5Xa0.67

1/41/51/6X30-52

CommAHdWindowCoIriiimandWlnicfow

1A&

?dF[-29,6A18;20,IZi-SjBjB]-EH巨IS

zu日iz

[VjD](瑞nobalaruce')-eas

ComnnaHdWihddu/

D二

1.0000636?40.3132

-25.310?00

-0.8533-1.00001.0000

0-10.51820

0.4390Q.69920.4642

0016.B351

(1)求方程的解;

(2)將方程右邊向量元素：改為0.53,在求解，并比較：的變化和解的相對變化;

(3)計算系數(shù)矩陣A的條件數(shù)并分析結(jié)論。

CommandWindow

CommandWindow

?lftl/4;：怎I?1/5:血血1/fi];

土口?悶。?記。?52f；

血El呱Q■肝；

血購山4XI二

張ira(iW

L2000

cl=t13.0000

o.

-6.6D0D

c2=3orjd(A)flora

6.6DD0

:1=:orM(再加:

(2):變大，其解中，相對未變化前的：的解：xl變大,x2變小，x3變大。

c3

cl=匚2二

2.0150e+003

2,0150&+003L3533e+Q03

⑶

由于A矩陣的條件數(shù)很大，故當(dāng)線性方程組中的b變大時，x也將發(fā)生很大的變

化，即數(shù)值穩(wěn)定性較差。

實驗三

3.硅谷公司員工的工資計算方法如下：

(1)工作時數(shù)超過120小時者，超過部分加發(fā)15%;

(2)工作時數(shù)低于60小時者，扣發(fā)700元；

(3)其余按每小時84元計發(fā)

Commithdwindew

I1CommandWindow

l>J"=11岳箱元L號(輸入0可赫翰JE

K-input(PWJA「寫入OFlSSSRIliA”，)iifx-0沽輸人工時數(shù)：1

In斑

end

t=nipu+C「％勺皴).

if+>120

y=12D*S4+(T-120)(l+fl.15)；

alseii?t*60

y=84*t-700.

P].SP

y-M't

y-

試編程按輸入的工號和該員工的工時數(shù)，計算

應(yīng)發(fā)工資-cie

and

詩輸入工號(輸入u可結(jié)求輸;2清輸入請輸入工號(輸入o可結(jié)束輸入！“3

工時數(shù)；同請輸入工時數(shù)；78

2

x

y

ond

40045552

實驗四

1根據(jù)1,求的近似值。當(dāng)n分別取1001000

2

n

VqEEondWindowCommandWllow

?y

11OU.21□□□；

1K1zn■for1=1:n

4/a.

?■nrl

10000時，結(jié)果是多少？要求：分

一

一別用循環(huán)結(jié)構(gòu)和向量運算來實現(xiàn)

3?13213.1400

CcnimjiHiiciWilidow

>>y=o.

1n=too：1-n=1000|;1'B=10000|;

2.Eg2-1=1：FL；2-i=1:n;

3-'2;3-f=L/l-*2;3-"2;

4-y=siun(f);4-y=suui(f)P4-y=suni(f).

匚5-sqrt(6*y)5-sqrt(67)

t-sqrt(6*y)

ans=

ans二arLS=

313213.3,1406考慮3,1415以下迭代公

式:Ya°其中a、b為正的常數(shù)

X

J00D0”

年H*i=1:n.

向量運算:

（1）編寫程序求迭代的結(jié)果，迭代的終止條件為，迭代初值xo=l.o,迭

代次數(shù)不超過500次

(2)如果迭代過程收斂于r,那么r的準(zhǔn)確值是b//24a,當(dāng)(a,

T~

b)的值取(1，1)'

(8,3)、(10,0.1)時，分別對迭代結(jié)果和準(zhǔn)確值進(jìn)行

比較

Conmandyvlndowxrl=V：jt-jd.

>>K=1,a=1,b_1.ifabs<Kl*X)<=10"<-5)brook

end

■fqxi=l;SQQend

■1

3cl二

CorrirrianclWihdc?w

?3E=1::b=3：

fori=1=600

KI=a/Hl-K);

i-f-sbs(x1—x)<=LCPC—5)hrII

end

1;

end

xl

CommandWindow

?K=L；aFao7b=0LL:

■feni=1:500

s:1----a/a

i-fabsCxi—3i.1<=IV'5)

bxcak

grid

x=xl;

end

3Kl

X1■

CommandWindowCommandWindowCommandWindow

?a=S;b=3;?'=if3\b=0.］；

r1=<-b+sqi;t(b"2+<*a))f2rI1-(—b^oq-ft(b~2+4*a))/2r1=(-b+sqrtCb"2+4*a)V2

r2=<-b-叩工t(b,2+4*a>)/2r^=(-b-sq.r-t(b2+4*aj)/Zr2=(-b-sqrtCb*2+4*a))/2

rl=i1-r1=

9180i-701S3.1127

r2=r2=r2=

-I.G130一勺?-3.2127

5.若兩個連續(xù)自然數(shù)的乘積減1是素數(shù)，則稱這兩個連續(xù)自然數(shù)是親密數(shù)對，

該素數(shù)是親密素數(shù)。例如，2X3-1=5是素數(shù)，所以2和3是親密數(shù)對，5是

親密素數(shù)。求［2,50區(qū)間內(nèi)：

(1)親密數(shù)對的對數(shù)。

(2)與上述親密數(shù)對對應(yīng)的所有親密素數(shù)之和。

1-

2-sum=0.

3-P=2:2500;

4-for1=2:00

5-n=findtrem(p,i)-=04A=i);

€-pW=［1i

7-end

電-P

S-fari=2:49

Ak=

11-if(find(m.==p)>Dfefind(m==p)<=leng+h(p))

28

12-k=k+l;

13-stnh=?ujif47i(find(m==p));

14-end

sun=

15-end

1C-k

21066

17-sum

實驗五

、實驗內(nèi)容

4設(shè)f(X)1-----------------L，編寫一個MATLAB函數(shù)

(X2)20.1(x3)'O.01

文件fx.m,使得調(diào)用f(x)時,x可用矩陣代入，得出的f(x)為同階矩陣。

1￡unction￡=fx匕)

2-fori=l:length(KC1::))

3-orj=I:length(K1:end))

4-f=j)-2)_*2+0.!)+!/(j)-3).4HL01);

5-end

6-endl

1Aunctianf=Af2(n)

2-f=nA(ix+l);

?A二口：閔21;3c3~56:23,34,8];

尸弦⑷

0.97160.31010.0028

0,00120.0011a.ooo3

0.0023o.aoio0.0293

5.已知f知0)

yf(30)f(20)

⑴當(dāng)f(n)n101n：5時，求y的值。

當(dāng)時f(n)122334nn1，求y的值。

⑴

1ftinctionf=f1(r)、、11(30)+f1(2Q))

2f=n+10i.(n."2代);

y=

o.6390

⑵

Comm-andWindow

1Jfunction

!>y「3(40X(f3(30)+f3(20y)

2i=f2(n-l)+f2(a),

y=

1.230&

實驗六

1.設(shè)3sinx,在x=0~2區(qū)間取101點，繪制

函數(shù)

y0.5------廠cosx

1x34

的曲線。

?x~linsi)ace(0,2*51,101);

plotCi,(0.5+3Hfsintz)/(1-ht.*2))*cos(x))

3?Aa-T;b=l,n=lA

t=llnspac&(0?2Api,360):p=a*sin(b-hrt*t);polar(tjp)

4?繪制極坐標(biāo)曲線aSinbn

并分析參數(shù)a、b、n對曲線形狀的影響

CommandWindow

為偶數(shù)時畫出的圖有該偶數(shù)的兩倍個環(huán)。參數(shù)a控制極坐標(biāo)的半徑，參數(shù)b

3、4、7時的情況，不難發(fā)現(xiàn)，當(dāng)n為奇數(shù)時畫出的圖有奇數(shù)個環(huán)，而當(dāng)n

可對圖進(jìn)行角度旋轉(zhuǎn)。

6.繪制曲面圖形，并進(jìn)行插值著色處理

xcosScostycos

Ssint0s7*7101

zsinS

2-t=linsj)ace((;

(SJl-ieshEnitsj);

4-K-COE.i：>COSIT)

5?AC05(3,*sin,T)

G-Fain(5);

r*

8-曲dingmtf:[:

ilabelfi-nisAylfibslC『曲)zAelC丁祖iF);

JU-

實驗七

2.利用曲面對象繪制曲面v(x,t)10e(，01xsin(2000t0.2x)

先利用默認(rèn)屬性繪制曲線，然后通過圖形句柄操作來改變曲線的顏色、線

型和線寬，并利用文字對象給曲線添加文字標(biāo)注

]1(0,10,1000);

2-t-1inspace(□,10,100);

T"I=iiQh-Pr4-d+)?

4-v=10*&xp(—0.01*E)-*sin(20D0*^pi*T-0?2*U+pi):

5-axesCviev\[-37.5,30]);

tis=surface(x,tjvi'FaceColor1'：'I：L<i-I)；

gridon;

Dause:

setChs,'FaceColor,'flat,);

實驗八

1.利用MATLAB提供的rand函數(shù)生成30000個符合均勻分布的隨機(jī)數(shù),

然后檢驗隨機(jī)數(shù)的性質(zhì)：

(1)均值和標(biāo)準(zhǔn)方差。

(2)最大元素和最小元素。

(3)大于0.5的隨機(jī)數(shù)個數(shù)占總數(shù)的百分比。

1|i=rand(3000c/n；

2-M=m$arL(A)

3-d(A)

4-E匕max(A)

5-F=jnin(A)

6-5))/30005

2.某氣象觀測站測得某日6:00~18:00之間每隔2h的室內(nèi)外溫度「C）

如實驗表1所示

實驗表1室內(nèi)外溫度觀測結(jié)果（C）

時間h681012141618

室內(nèi)溫度tl18.020.022.025.030.028.024.0

室外溫度t215.019.024.028.034.032.030.0

試用三次樣條插值分別求出該日室內(nèi)外6:30^17:30之間每隔2h各

點的近似溫度（C）

t-

2-hl=6.5:2:17.5;

3-11=[18,0,20.0,22,0,25,0,30.0,28.Q,24.0];

4t2=[15*C,19-0,2'.0,28.0,34.C,32.0,30.0].

5-七3=Eplin已⑴】七1,hl）

2

5有3個多項式，p（x）x2x4x5'p2?x2,時進(jìn)

6-t4=spline（hj

CommandWirdow

t3=tl二

Coluiuis]ihrcmghi4Colujins[thrAUfh◎

二50202山解鋌2占弭毋2-3”7515.655320.3366常第[購2工時彩

Lolmne5ikiTongh6Colimns5through6

30.205126134,255B30,3594

行下列操作：(D求p(x)p'X)P2(X)P3(X)o

(2)求P(x)的根。

(3)當(dāng)x取矩陣A的每一元素時，求P(x)的值。其中:

11.21.4

A0.7523.5

052.5

(4)當(dāng)以矩陣A為自變量時，求P(x)的值。其中A的值與第(3)

題相同。

1-pl=Ll.2.4,0,SJ：

2-p2;EL2J:

3-p3=[l,2.,3]:

4-￡=[%conv(p2jp3)];

5-p=pl+￡

6-s(p)

7L2-1.4.0.‘5,2二5:0,5,2.5].

8-xIL一polyval(pjA)(2)

9-x?=polyvalm(p,A)

(1)

Window-L3840+L83171

-1.3840-L83171

-0.1160+1.44001

-0.1160-1.44001

(3)

(2)

L0ft+003*

xi=

0.0076-0.1281-0.07?E

t.0&+003*

0.1J281.3900i.isaa

0.用羽1-73641,5136

Q.01000.03820.0125

0.02230.09700,4122

0.01101,.24600,LS44

實驗九

1.求函數(shù)在指定點的數(shù)值導(dǎo)數(shù)。

XXx

2x1,2,3

f(x)12xX

Q6x

2.用數(shù)值為法求定積分

2I22

()Iiocost4sin(2t)川的近似值。

1-0.Cl：2*pi:

2-l=5^+(cosf-t.?'2)+4*sin(.”2)+1):

3-Itrapzd,!)6*79S3+3*148fii

3.分別用3