2021年matlab數(shù)學(xué)實(shí)驗(yàn)報(bào)告_2

1、精編word文檔 下載可編輯matlab數(shù)學(xué)實(shí)驗(yàn)報(bào)告班級(jí)09b姓名陳榮菲學(xué)號(hào)201*1611223實(shí)驗(yàn)名齊次與非齊次線(xiàn)性方程組的解稱(chēng)問(wèn)題背景描述在學(xué)習(xí)高等代數(shù)的過(guò)程中，我們常常需要求解一些線(xiàn)性方程組，而這一過(guò)程又必須經(jīng)常涉及矩陣的有關(guān)計(jì)算，很多同學(xué)都抱怨這些計(jì)算太過(guò)繁瑣，費(fèi)時(shí)且結(jié)果也不一定正確，但是計(jì)算機(jī)恰恰彌補(bǔ)了這個(gè)不足，使你可以很方便的求到方程組的解。實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)、掌握matlab軟件的有關(guān)命令；掌握線(xiàn)性方程組的求解。實(shí)驗(yàn)原理與數(shù)學(xué)模型對(duì)于齊次線(xiàn)性方程組而言，當(dāng)系數(shù)行列式的秩等于未知數(shù)的個(gè)數(shù)時(shí)，方程組只有零解，否則有無(wú)窮解，只要找到其中一組解，在其基礎(chǔ)上乘以一常數(shù)k就得到方程組的解。對(duì)于非

2、齊次線(xiàn)性方程組而言，當(dāng)系數(shù)矩陣的秩=增廣矩陣的秩=未知數(shù)的個(gè)數(shù)時(shí)，方程組有唯一解；當(dāng)系數(shù)矩陣的秩=增廣矩陣的秩2x13x2x30（2）4x12x2x30x13x30解(1)a=1-111;1-11-2;1-1-21;formatratn=4;ra=rank(a)ra=3if(ra=n)elseb=null(a,r)endb=1100symskx=k*bx=kk00(2)a=23-1;4-21;103;formatratn=3;ra=rank(a)ra=3if(ra=n)x=000elseb=null(a,r)endx=000求解下列非齊次線(xiàn)性方程組4x12x2x323x1x22x31011x1

3、3x28x38解a=42-1;3-12;1138;b=2108;b=ab;n=3;ra=rank(a)ra=3rb=rank(b)rb=3if(ra=rb&ra=n)x=abelseif(ra=rb&ra002的解為x=97/40-169/40-3/4實(shí)驗(yàn)總結(jié)這是我第一次用matlab進(jìn)行編程來(lái)求解實(shí)際問(wèn)題，雖然過(guò)程有點(diǎn)艱辛，但每一步都親力親為，這讓我收獲很多，通過(guò)做次實(shí)驗(yàn)，讓我對(duì)matlab有了進(jìn)一步的了解，了解了它的強(qiáng)大的功能和他如何求解實(shí)際問(wèn)題，激發(fā)了我學(xué)好matlab的決心。思考與深入；求解線(xiàn)性方程組還有沒(méi)有別的方法嗎？答案是有的，對(duì)于第2題這里給出另外一種方法a=42-1;3-12;

4、1138;b=2108;b=ab;n=3;ra=rank(a)ra=3formatratrb=rank(b)rb=3l,u=lu(a)l=4/11-1/213/1110100u=11380-20/11-2/1100-4x=u(lb)x=97/40-169/40-3/4教師評(píng)語(yǔ)擴(kuò)展閱讀matlab數(shù)學(xué)實(shí)驗(yàn)報(bào)告-1-matlab作業(yè)matlab與科學(xué)計(jì)算實(shí)驗(yàn)報(bào)告指導(dǎo)老師易昆南實(shí)驗(yàn)時(shí)間大二年級(jí)上學(xué)期學(xué)院數(shù)學(xué)科學(xué)與計(jì)算技術(shù)學(xué)院專(zhuān)業(yè)班級(jí)統(tǒng)計(jì)0902姓名丁克明學(xué)號(hào)1304090112-2-matlab作業(yè)目錄學(xué)生成績(jī)管理3函數(shù)極限問(wèn)題5蛛網(wǎng)模型7水塔模型10混沌問(wèn)題12koch曲線(xiàn)14拉格朗日與三次樣條插

5、值以及曲線(xiàn)擬合17解線(xiàn)性方程組20矩陣對(duì)角化2310.矩陣標(biāo)準(zhǔn)化.241摸球?qū)嶒?yàn)261釘板問(wèn)題261火車(chē)問(wèn)題28-3-matlab作業(yè)一、學(xué)生成績(jī)管理學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南學(xué)生成績(jī)管理程序1、設(shè)計(jì)（實(shí)習(xí)）目的1結(jié)合實(shí)際問(wèn)題展現(xiàn)matlab在生活和學(xué)習(xí)方面的廣泛應(yīng)用2學(xué)會(huì)利用matlab編程并求解實(shí)際問(wèn)題3學(xué)會(huì)并運(yùn)用for循環(huán)和switch結(jié)構(gòu)，以及matlab中已有函數(shù)如sum4了解單元數(shù)組cell和結(jié)構(gòu)數(shù)組struct的作用，學(xué)會(huì)創(chuàng)建這些數(shù)組.2、實(shí)驗(yàn)內(nèi)容編寫(xiě)一個(gè)學(xué)生成績(jī)管理的函數(shù)程序?qū)W(xué)生成績(jī)按五個(gè)等級(jí)劃分優(yōu)秀（90以上），良好（80-

6、90），中等（70-79），及格（60-69），不及格（60以下）對(duì)輸入的學(xué)生成績(jī)按五個(gè)等級(jí)劃分后，打印學(xué)生姓名，得分，等級(jí)，并計(jì)算全班最高分、最低分、平均分。3、詳細(xì)設(shè)計(jì)clearfori=1:10ai=89+i;bi=79+i;ci=69+i;di=59+i;endc=d,c;name=input(pleaseinputname:);score=input(pleaseinputscore:);n=length(score);rank=cell(1,n);s=struct(name,name,score,score,rank,rank);fori=1:nswitchs(i).scoreca

7、se100s(i).rank=滿(mǎn)分;caseas(i).rank=優(yōu)秀;casebs(i).rank=良好;casecs(i).rank=及格;caseds(i).rank=不及格;endenddisp(學(xué)生姓名,得分,等級(jí))；-4-matlab作業(yè)fori=1:ndisp(s(i).name,blanks(6),num2str(s(i).score),blanks(6),s(i).rank);ends=0;fori=1:ns=s(i).score+s;endaverscore=s/n;t=s(1).score;fori=1:(n-1)if(s(i).scores(i+1).score)m=s

8、(i+1).score;endenddisp(平均成績(jī));disp(averscore);disp(最高分);disp(t);disp(最低分);disp(m);4、實(shí)驗(yàn)結(jié)果-5-matlab作業(yè)5、實(shí)驗(yàn)總結(jié)這是我第一次接觸matlab并運(yùn)用matlab解決實(shí)際問(wèn)題，在以前學(xué)習(xí)c+的時(shí)候，解決一個(gè)問(wèn)題的代碼會(huì)很繁瑣。但相同的問(wèn)題，在matlab中卻變得很簡(jiǎn)單。這就是matlab的強(qiáng)大之處。運(yùn)用它可以方便的解決許多實(shí)際問(wèn)題。知道了這一點(diǎn)，我決心以后會(huì)認(rèn)真的學(xué)習(xí)這樣一門(mén)課程。同時(shí)，它也給我們枯燥的數(shù)學(xué)公式的學(xué)習(xí)帶來(lái)了動(dòng)手解決實(shí)際問(wèn)題的機(jī)會(huì)。是將理論與實(shí)際相結(jié)合的方法。讓我們更加體會(huì)到知識(shí)的力量是強(qiáng)

9、大的，我們應(yīng)該更好的掌握科學(xué)技術(shù)和相關(guān)理論知識(shí)，并能夠?qū)⑵溥\(yùn)用于實(shí)際生活當(dāng)中。在解決這個(gè)問(wèn)題的時(shí)候也確實(shí)遇到了一些問(wèn)題。比如如何輸入符號(hào)數(shù)組和數(shù)值數(shù)組的問(wèn)題，當(dāng)時(shí)一直不理解。當(dāng)通過(guò)查詢(xún)資料，和同學(xué)討論，最后明白了。輸入符號(hào)時(shí)，符號(hào)要用單引號(hào)，而數(shù)據(jù)就不需要了。通過(guò)這樣一次實(shí)踐，我更加明白了動(dòng)手的必要性。只有理論知識(shí)是遠(yuǎn)遠(yuǎn)不夠的。所以以后一定要加強(qiáng)自己的動(dòng)手能力，勤動(dòng)手。二、函數(shù)極限問(wèn)題學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南函數(shù)極限問(wèn)題1、設(shè)計(jì)（實(shí)習(xí)）目的1結(jié)合實(shí)際問(wèn)題展現(xiàn)matlab在生活和學(xué)習(xí)方面的廣泛應(yīng)用2學(xué)會(huì)利用matlab編程并求解實(shí)際問(wèn)題3學(xué)會(huì)并

10、運(yùn)用movie和moviein函數(shù)，以及plot結(jié)構(gòu)。4理解getframe以及學(xué)會(huì)運(yùn)用axis調(diào)節(jié)坐標(biāo)抽2、實(shí)驗(yàn)內(nèi)容自選函數(shù)，運(yùn)用plot進(jìn)行繪圖。并使用movie以及moviein，制作函數(shù)極限動(dòng)畫(huà)。3、詳細(xì)設(shè)計(jì)clearm=moviein(16);holdonfori=1:50forj=1:it1=(j-1)*5;t2=j*5;t=t1:0.01:t2;axis(0,40,-1,1)%限制動(dòng)畫(huà)的坐標(biāo)顯示大小plot(t,0,.r)plot(t,1.)endm(:,i)=getframe;endmovie(m,50)clearm=moviein(16)-6-matlab作業(yè)holdonfo

11、ri=1:10forj=1:it1=(j-1)*pi/10;t2=j*pi/10;t=t1:0.01:t2;axis(0,5*pi/2,-1,1)%限制動(dòng)畫(huà)的坐標(biāo)顯示大小plot(pi-0.1,t,.k)plot(pi-0.1,-t,.k)plot(pi+0.1,-t,.k)plot(pi+0.1,t,.k)plot(t,0,.r)plot(t,sin(t),.)t1=(20-j)*pi/10;t2=(21-j)*pi/10;t=t1:0.01:t2;axis(0,5*pi/2,-1,1)%限制動(dòng)畫(huà)的坐標(biāo)顯示大小plot(t,0,.r)plot(t,sin(t),.)endm(:,i)=get

12、frame;endmovie(m,50)4、實(shí)驗(yàn)結(jié)果-7-matlab作業(yè)5、實(shí)驗(yàn)總結(jié)使用matlab進(jìn)行繪圖覺(jué)得很有意思，簡(jiǎn)單的代碼，卻能繪出各種顏色的圖形動(dòng)畫(huà)。在學(xué)習(xí)的過(guò)程中還是有很多的問(wèn)題。還有很多方面的東西不夠了解，只能邊嘗試邊查詢(xún)資料，讓自己能夠更了解，更明白。以便繪出更好的圖形。三、蛛網(wǎng)模型學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南蛛網(wǎng)模型1、設(shè)計(jì)（實(shí)習(xí)）目的學(xué)會(huì)運(yùn)用matlab解決實(shí)際問(wèn)題進(jìn)一步掌握matlab的一些基本操作通過(guò)范例體會(huì)matlab的初步建模過(guò)程4學(xué)會(huì)用做動(dòng)畫(huà)2、實(shí)驗(yàn)內(nèi)容已知前兩年的豬肉的產(chǎn)量和豬肉的價(jià)格分別為39噸，28噸，12

13、元/公斤，17元/公斤，根據(jù)前一年的豬肉價(jià)格影響后一年豬肉產(chǎn)量當(dāng)年豬肉產(chǎn)量影響當(dāng)年豬肉價(jià)格的線(xiàn)性關(guān)系，編寫(xiě)程序，利用動(dòng)畫(huà)原理畫(huà)出前十年豬肉的產(chǎn)量?jī)r(jià)格的動(dòng)態(tài)圖形。（參數(shù)設(shè)置為，-8-matlab作業(yè)c1=39,c2=28,c3=36,r1=12,r2=17,k=16)。3、詳細(xì)設(shè)計(jì)clear%c1為產(chǎn)量1,c2為產(chǎn)量2,c3為產(chǎn)量3,r1為%肉價(jià)1,r2為肉價(jià)2,k為k年后產(chǎn)量與肉價(jià)%是否穩(wěn)定holdoffc1=39;c2=28;c3=36;r1=12;r2=17;k=16;a1=c11;c21;%系數(shù)矩陣b1=r1,r2;%列向量a2=r11;r21;b2=c2,c3;a=a1b1;b=a2b

14、2;%x0(1)=c1;forn=1:30y0(n)=a(1)*x0(n)+a(2);%a(1)a(2)為矩陣元素的引用，下行類(lèi)似x0(n+1)=b(1)*y0(n)+b(2);x(n)=x0(n);y(n)=x0(n+1);endplot(x,y0,-g,y,y0,-b)holdonm=moviein(100);forn=1:kforj=1:30t1=x0(n)+(j-1)*(x0(n+1)-x0(n)/30;t2=x0(n)+j*(x0(n+1)-x0(n)/30;ift2-9-matlab作業(yè)t2=t;elseendt=t1:0.01:t2;plot(x(n+1),t,.r)%劃豎線(xiàn)en

15、dm(:,n)=getframe;endmovie(m,20)4、實(shí)驗(yàn)結(jié)果ans=columns1through630000200003000030.18183413233358columns7through12357393946231300326913895234398columns13through18377103530237053357793670636032columns19through24365223616636425362363637436274columns25through30363463629436332363043632436310column31363205、實(shí)驗(yàn)圖像8

16、864278642635353635373538353935-10-matlab作業(yè)6、實(shí)驗(yàn)總結(jié)通過(guò)做此實(shí)驗(yàn)，讓我對(duì)matlab有更進(jìn)一步的了解，學(xué)會(huì)怎樣才能正確運(yùn)用matlab求解實(shí)際問(wèn)題，了解如何利用數(shù)學(xué)模型去解釋和分析社會(huì)經(jīng)濟(jì)問(wèn)題，特別是這個(gè)典型經(jīng)濟(jì)問(wèn)題的求解。我對(duì)matlab還不是特別熟悉，過(guò)程中遇到了很多問(wèn)題，經(jīng)過(guò)與同學(xué)討論得到了良好的解決，希望以后會(huì)有進(jìn)步。四、水塔模型學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南水塔問(wèn)題1、設(shè)計(jì)（實(shí)習(xí)）目的1掌握用matlab來(lái)求最值2學(xué)會(huì)在一個(gè)閉區(qū)間求最值3學(xué)會(huì)解決實(shí)際問(wèn)題。2、實(shí)驗(yàn)內(nèi)容水塔問(wèn)題在地面上建有一座圓柱

17、形水塔,水塔內(nèi)部的直徑為d,并且在地面處開(kāi)了一個(gè)高為h的小門(mén).現(xiàn)在要對(duì)水塔內(nèi)部進(jìn)行維修施工,施工方案要求把一根長(zhǎng)為l(ld)的水管運(yùn)到水塔內(nèi)部.請(qǐng)問(wèn)水塔的門(mén)高h(yuǎn)多高時(shí),才有可能成功地把水管搬進(jìn)水塔內(nèi)3、詳細(xì)設(shè)計(jì)一clearfunctionx=lt523(l,d,h)k1=0;a=0l=7;d=3ymax=0;k2=0;b=pi/2;h=3;while(b-a)10(-2)k2=k2+1;m=0;a=0;ifk1=0n=ceil(b/h)-1elsen=ceil(b/h);endfori=1:n+1;-11-matlab作業(yè)x(i)=a+(i-1)*(b-a)/n;y(i)=l*sin(x(i)

18、-d*tan(x(i);endfori=1:nify(i)ymaxymax=y(i);a=x(i-1);elsey(i)-12-matlab作業(yè)plot(x,8*sin(x)-2*tan(x),-b)4、實(shí)驗(yàn)圖像5、實(shí)驗(yàn)總結(jié)通過(guò)做此實(shí)驗(yàn)，讓我對(duì)matlab有更進(jìn)一步的了解，學(xué)會(huì)怎樣才能正確運(yùn)用matlab求解實(shí)際問(wèn)題，了解如何利用數(shù)學(xué)模型去解釋和分析社會(huì)經(jīng)濟(jì)問(wèn)題，特別是這個(gè)典型經(jīng)濟(jì)問(wèn)題的求解。我對(duì)matlab還不是特別熟悉，過(guò)程中遇到了很多問(wèn)題，經(jīng)過(guò)與同學(xué)討論得到了良好的解決，希望以后會(huì)有進(jìn)步。五、混沌問(wèn)題學(xué)號(hào)1304090112實(shí)驗(yàn)題目1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)

19、會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容利用迭代方法繪制混沌圖形3、詳細(xì)設(shè)計(jì)x0=0.3;holdonfork=1:100班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南混沌問(wèn)題-13-matlab作業(yè)a=3+(k-1)/100;fori=1:300x0=a*x0*(1-x0);ifi100plot(a,x0,r.)endendend4、實(shí)驗(yàn)圖像5、實(shí)驗(yàn)總結(jié)實(shí)際的繪圖進(jìn)一步讓我體會(huì)到了matlab功能的強(qiáng)大。他能在實(shí)際生活中解決很多人工無(wú)法解決的復(fù)雜問(wèn)題。通過(guò)做此實(shí)驗(yàn)，讓我對(duì)matlab有更進(jìn)一步的了解，學(xué)會(huì)怎樣才能正確運(yùn)用matlab求解實(shí)際問(wèn)題，了解如何利用數(shù)學(xué)模型去解釋和分析社會(huì)經(jīng)濟(jì)

20、問(wèn)題，特別是這個(gè)典型經(jīng)濟(jì)問(wèn)題的求解。我對(duì)matlab還不是特別熟悉，過(guò)程中遇到了很多問(wèn)題，經(jīng)過(guò)與同學(xué)討論得到了良好的解決，希望以后會(huì)有進(jìn)步。-14-matlab作業(yè)六、koch曲線(xiàn)學(xué)號(hào)1304090112實(shí)驗(yàn)題目1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容繪制雪花曲線(xiàn)和其他一個(gè)美麗圖形3、詳細(xì)設(shè)計(jì)1、雪花曲線(xiàn)p=00;100;n=2;a=cos(pi/3)-sin(pi/3);sin(pi/3)cos(pi/3);fork=1:4d=diff(p)/3;m=4*n-3;q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);

21、p(2:4:m,:)=q+d;p(3:4:m,:)=q+d+d*a;p(4:4:m,:)=q+2*d;n=m;endplot(p(:,1),p(:,2)axis(010010)班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南koch曲線(xiàn)2、p=010;100;0-10;-100;010;n=5;a=0-1;10;fork=1:5d=diff(p)/3;m=4*n-3;q=p(1:n-1,:);p(5:4:m,:)=p(2:n,:);p(2:4:m,:)=q+d;p(3:4:m,:)=q+2*d+d*a;p(4:4:m,:)=q+2*d;n=m;endplot(p(:,1),p(:,2)axis(-

22、1010-1010)3、花草樹(shù)木-15-matlab作業(yè)p=50;510;n=2;a=cos(pi/3)-sin(pi/3);sin(pi/3)cos(pi/3);b=cos(-pi/3)-sin(-pi/3);sin(-pi/3)cos(-pi/3);fork=1:4d=diff(p)/3;d1=d(1:2:n,:);m=5*n;q1=p(1:2:n-1,:);p(10:10:m,:)=p(2:2:n,:);p(1:10:m,:)=p(1:2:n,:);p(2:10:m,:)=q1+d1;p(3:10:m,:)=p(2:10:m,:);p(4:10:m,:)=q1+d1+d1*a;p(5:1

23、0:m,:)=p(2:10:m,:);p(6:10:m,:)=q1+2*d1;p(7:10:m,:)=p(6:10:m,:);p(8:10:m,:)=q1+2*d1+d1*b;p(9:10:m,:)=p(6:10:m,:);n=m;endplot(p(:,1),p(:,2)axis(010010)4、實(shí)驗(yàn)圖像1、雪花曲線(xiàn)-16-matlab作業(yè)2、3、-17-matlab作業(yè)5、實(shí)驗(yàn)總結(jié)通過(guò)做此實(shí)驗(yàn)，讓我對(duì)matlab有更進(jìn)一步的了解，學(xué)會(huì)怎樣才能正確運(yùn)用matlab求解實(shí)際問(wèn)題，了解如何利用數(shù)學(xué)模型去解釋和分析社會(huì)經(jīng)濟(jì)問(wèn)題，特別是這個(gè)典型經(jīng)濟(jì)問(wèn)題的求解。我對(duì)matlab還不是特別熟悉，過(guò)程中

24、遇到了很多問(wèn)題，經(jīng)過(guò)與同學(xué)討論得到了良好的解決，希望以后會(huì)有進(jìn)步。七、拉格朗日插值與三次樣條插值以及曲線(xiàn)擬合學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師易昆南評(píng)分拉格朗日插值與三次樣條插值以及曲線(xiàn)擬合1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容拉格朗日插值和曲線(xiàn)擬合（擬合）用下面一組數(shù)據(jù)擬合c(t)abe0.0.2kt中的參數(shù)a，b，ktj100201*004005006007008009001000ci549935659010263950593、詳細(xì)設(shè)計(jì)拉格朗日插值x1=1:0.1:5;x0=1,3,5;y0=

25、1,9,25;x=2,5,3,5,4,5;fori=1:6s=0;fork=1:3p=1;forj=1:3ifj=kp=p*(x(i)-x0(j)/(x0(k)-x0(j);endend-18-matlab作業(yè)s=p*y0(k)+s;endy(i)=s;endplot(x1,x2,b,x,y,*r)三次樣條插值x0=0:0.1:5;y0=x0.2;x=0:0.5:5;y=interp1(x0,y0,x,spline);plot(x0,y0,b,x,y,*r)曲線(xiàn)擬合f=inline(x(1)+x(2)*exp(0.02*x(3)*tdata),x,tdata);tdata=100:100:10

26、00;cdata=54,99,35,65,90,10,26,39,50,59;x,jm=lsqcurvefit(f,111,tdata,cdata);xjmplot(tdata,cdata)4、實(shí)驗(yàn)圖像拉格朗日插值-19-matlab作業(yè)三次樣條插值擬合-20-matlab作業(yè)5、實(shí)驗(yàn)總結(jié)這次的試驗(yàn)個(gè)人覺(jué)得就拉格郎日的振蕩比較難做，雖然寫(xiě)出來(lái)了一眼看的就很清楚，不過(guò)其思想還是找了挺多資料后在同學(xué)的幫助下才得以完成；至于其他的試驗(yàn)讓我學(xué)習(xí)到了更多的知識(shí)，開(kāi)闊了眼界，明白了多種方法，對(duì)于以后的解題思路有了更多的思考余地，這次試驗(yàn)中用到了拉格郎日插值法，分段線(xiàn)性插值，三次樣條插值，擬合的基本原理八、

27、解線(xiàn)性方程組學(xué)號(hào)1304090112實(shí)驗(yàn)題目班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南解線(xiàn)性方程組1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容1/比較用逆矩陣法、除法、克拉默法則解方程ax=b的用時(shí)和誤差2/求解方程組最多零解3、詳細(xì)設(shè)計(jì)1、a=rand(100)*e2;x=ones(100,1);b=a*x;ticy=inv(a)*b;tocerr=norm(y-x)res=norm(a*y-b)ticy=ab;tocerr=norm(y-x)res=norm(a*y-b)tica=det(a)fori=1:1-21-matlab

28、作業(yè)b=a;b(1:100,i)=b;y(i)=det(b)/a;endtocerr=norm(y-x)res=norm(a*y-b)2、functionc=solution(a)l,u=lu(a);m1=;m2=;e=;f=;c=;r=rank(a);rf=size(a);n=rf(1,2);k=1;fori=1:rwhilek-22-matlab作業(yè)elapsedtimeis0.004000seconds.err=9797e-012res=1628e-010elapsedtimeis0.003000seconds.err=0572e-012res=1790e-011a=-3898e+225

29、elapsedtimeis0.083000seconds.err=0860e-012res=4895e-0092、-23-matlab作業(yè)5、實(shí)驗(yàn)總結(jié)通過(guò)本次實(shí)驗(yàn)，我基本上掌握了matlab求矩陣的秩、行列式和逆矩陣的命令，并通過(guò)編程比較可知用逆矩陣法、除法、克拉默法則求解方程ax=b時(shí)除法最省時(shí)，最精確；而且利用了以上方法求解給定的恰定，超定和欠定線(xiàn)性方程組，可以很快的求出結(jié)果；可見(jiàn)，利用matlab求解線(xiàn)性方程組是非常方便的，很有利于我們各方面涉及到這類(lèi)問(wèn)題的學(xué)習(xí)和工作。而由于本來(lái)對(duì)線(xiàn)性代數(shù)的知識(shí)就掌握得不是很好，所以做起這樣的題目來(lái)確實(shí)很吃力，不過(guò)使用matlab軟件解決了計(jì)算煩瑣的問(wèn)題

30、，所以相對(duì)來(lái)說(shuō)，題目也基本上能得到解決。在做實(shí)驗(yàn)過(guò)程中，我對(duì)線(xiàn)性代數(shù)的知識(shí)又有了進(jìn)一步的認(rèn)識(shí)，而且對(duì)于它結(jié)合matlab在現(xiàn)實(shí)問(wèn)題中的應(yīng)用也有了初步的了解，這也是一個(gè)不小的收獲吧。九、矩陣對(duì)角化學(xué)號(hào)1304090112實(shí)驗(yàn)題目1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容矩陣對(duì)角化3、詳細(xì)設(shè)計(jì)functiony=trigle(a)y=1;c=size(a);ifc(1)=c(2)y=0;return;ende=eig(a);n=length(a);while1ifisempty(e)return;endd=e(1);f=sum(abs(

31、e-d)-24-matlab作業(yè)y=0;return;ende(find(abs(e-d)-25-matlab作業(yè)1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容對(duì)矩陣進(jìn)行標(biāo)準(zhǔn)化3、詳細(xì)設(shè)計(jì)a=1-11;-1-3-3;1-30p,d=eig(a)symsy1y2y3;y=y1;y2;y3f=y*d*y4、實(shí)驗(yàn)結(jié)果a=1-11-1-3-31-30p=-0.0580-0.8018-0.5948-0.8554-0.26730.4437-0.51470.5345-0.6703d=-87300000.00000008730y=y1y2y3f=-34

32、2905718471287/70368744177664*conj(y1)*y1+3196254119/19342813113834066795298816*conj(y2)*y2+3234691681855341/1125899906842624*conj(y3)*y35、實(shí)驗(yàn)總結(jié)這個(gè)實(shí)驗(yàn)進(jìn)一步讓我感受到了matlab的實(shí)用性，以及其在數(shù)學(xué)無(wú)論是代數(shù)還是幾何整個(gè)領(lǐng)域里的極大用處。個(gè)人覺(jué)得，學(xué)習(xí)數(shù)學(xué)專(zhuān)業(yè)的我們大家應(yīng)該很好的學(xué)學(xué)這樣一個(gè)多功能的軟件。以便在以后的生活工作中游刃有余。-26-matlab作業(yè)1摸球?qū)嶒?yàn)學(xué)號(hào)1304090112實(shí)驗(yàn)題目1、設(shè)計(jì)（實(shí)習(xí)）目的了解matlab在實(shí)際問(wèn)題中的

33、應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容不同情況下，計(jì)算條件摸球的概率。一共有是個(gè)球三個(gè)黑球，七個(gè)白球。不放回摸球。1、第三次摸到黑球2、第三次才摸到黑球3、三次都摸到黑球3、詳細(xì)設(shè)計(jì)a=rand(1000000,3);a(:,1)=round(a(:,1)-0.2);a(:,2)=round(a(:,2)*0.9-0.2-0.1*(a(:,1)-1);a(:,3)=round(a(:,3)*0.8-0.2-0.1*(a(:,1)-1)-0.1*(a(:,2)-1);fori=1:6b(i)=sum(a(1:10i,3)/(10i);m=a(1:10i,1)&a(1:10i,2)

34、&a(1:10i,3);c(i)=sum(m)/10i;d=a(1:10i,1)&a(1:10i,2)&a(1:10i,3);e(i)=sum(d)/(10i);endbce4、實(shí)驗(yàn)結(jié)果b=0.40000.2201*.26400.30070.30130.3004c=0.10000.14000.14500.17560.17580.1752e=000.00300.00800.00790.00835、實(shí)驗(yàn)總結(jié):第一次將軟件和正在學(xué)習(xí)的概率論理論知識(shí)聯(lián)系在一起。將理論在計(jì)算機(jī)中模擬實(shí)現(xiàn)。覺(jué)得非常有意思。更加加強(qiáng)了我學(xué)習(xí)matlab的信心。雖然這個(gè)學(xué)期快要結(jié)束，雖然matlab課程快要結(jié)束，但我相信自己會(huì)繼續(xù)學(xué)習(xí)這個(gè)軟件。班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南摸球?qū)嶒?yàn)1釘板問(wèn)題學(xué)號(hào)1304090112實(shí)驗(yàn)題目1、設(shè)計(jì)（實(shí)習(xí)）目的班級(jí)統(tǒng)計(jì)0902姓名丁克明指導(dǎo)教師評(píng)分易昆南釘板問(wèn)題-27-matlab作業(yè)了解matlab在實(shí)際問(wèn)題中的應(yīng)用學(xué)會(huì)利用matlab做圖并求解實(shí)際問(wèn)題2、實(shí)驗(yàn)內(nèi)容釘板問(wèn)題3、詳細(xì)設(shè)計(jì)c