1、、編程實(shí)現(xiàn)以下科學(xué)計(jì)算算法，并舉一例應(yīng)用之?！袄窭嗜粘俗臃s束最優(yōu)化”拉格朗日乘子法求約束最優(yōu)化問(wèn)題實(shí)例。采用拉格朗日乘子法如下最優(yōu)化問(wèn)題mini(x,九)二x+x+九(x2+x2)1211在MATLAB中編寫(xiě)函數(shù)ex1208.m來(lái)進(jìn)行求解，具體代碼如下所示。%ex1208.m拉格朗日乘子法求最優(yōu)化解x=zeros(1,2)symsxylamaf=x+y+lama*(xA2+yA2-2);dx=diff(f,x);dy=diff(f,y);dlama=diff(f,lama);xx=solve(dx,x);yy=solve(dy,y);ff=subs(dlama,x,y,xx,yy);lam

2、ao=solve(ff);xo=subs(xx,lama,lamao)yo=subs(yy,lama,lamao)fo=subs(f,x,y,lama,xo,yo,lamao)程序運(yùn)行結(jié)果為：xo=1-1yo=1-1fo=2-2%用syms表示出轉(zhuǎn)化后的無(wú)約束函數(shù)%分別求函數(shù)關(guān)于x、y、lama的偏導(dǎo)%令偏導(dǎo)為零，求解x、y%將x表示為lama函數(shù)%將y表示為lama函數(shù)%代入dlama得關(guān)于lama的一元函數(shù)%求解得lama0%求得取極值處的x0%取極值處的y0%取極值處的函數(shù)值流程圖：開(kāi)始丨/fflsyms表示出轉(zhuǎn)化/后6勺無(wú)約束函數(shù)/別求函數(shù)關(guān)于/x、y、lama的偏號(hào)/令偏導(dǎo)為零/求解

3、x、ydlama得關(guān)帀/Lama的一元函哪/金解得lamaO極值易/的xO、yO和函數(shù)粵/結(jié)束二、編程解決以下科學(xué)計(jì)算和工程實(shí)際問(wèn)題。、1、利用MATLAB提供的randn函數(shù)聲稱(chēng)符合正態(tài)分布的105隨機(jī)矩陣A,進(jìn)行如下操作:A各列元素的均值和標(biāo)準(zhǔn)方差。A的最大元素和最小元素。求A每行元素的和以及全部元素之和。分別對(duì)A的每列元素按升序、每行元素按降序排序。代碼：clearall;closeall;clc;A=randn(10,5);meanA=mean(A);stdA=std(A);maxA=max(max(A);minA=min(min(A);rowsumA=sum(A,2);sumA=su

4、m(rowsumA);sort1=sort(A);sort2=sort(A,2,descend);%(1)A各列元素的均值%(1)A各列元素的標(biāo)準(zhǔn)方差%(2)A的最大元素%(2)A的最小元素%(3)A每行元素的和%(3)A全部元素的和%(4)A的每列元素按升序排列%(4)A的每列元素按降序排列運(yùn)行結(jié)果：因生成矩陣隨機(jī)，故無(wú)固定結(jié)果流程圖：2、按要求對(duì)指定函數(shù)進(jìn)行插值和擬合。（1）按表6.4用3次樣條方法插值計(jì)算090度范圍內(nèi)整數(shù)點(diǎn)的正弦值和075度范圍內(nèi)整數(shù)點(diǎn)的正切值，然后用5次多項(xiàng)式擬合方法計(jì)算相同的函數(shù)值，并將兩種計(jì)算結(jié)果進(jìn)行比較。表6.4特殊角的正弦與正切值表a(度)0153045607

5、590asin00.25880.50000.70710.86600.96591.0000+atan00.26790.57741.00001.73203.73202）按表6.5用3次多項(xiàng)式方法插值計(jì)算1100之間的整數(shù)的平方根。表6.51100內(nèi)特殊值的平方根表N149162536496481100vW12345678910(1)代碼：clearall;closeall;clc;alpha1=0:15:90;sin_alpha1=sin(alpha1*pi/180);plot(alpha1,sin_alpha1,k:p);alpha2=0:90;sin_Y1=interp1(alpha1,sin

6、_alpha1,alpha2,spline);plot(alpha2,sin_Y1,r-*);holdon;P1=polyfit(alpha1,sin_alpha1,5);sin_Y2=polyval(P1,alpha2);plot(alpha2,sin_Y2,b-o);legendalpha3=0:15:75;tan_alpha3=tan(alpha3*pi/180);figure,plot(alpha3,tan_alpha3,k:p);holdon;alpha4=0:75;tan_Y1=interp1(alpha3,tan_alpha3,alpha4,spline);plot(alpha4

7、,tan_Y1,r-*);holdon;P2=polyfit(alpha3,tan_alpha3,5);tan_Y2=polyval(P2,alpha4);plot(alpha4,tan_Y2,b-o)legend運(yùn)行結(jié)果：正切值比較流程圖：5次多項(xiàng)式擬兮務(wù)次多項(xiàng)式擬合）結(jié)束求3次樣條正弦插值次多項(xiàng)式求夢(mèng)/4次多項(xiàng)式求粵/繪精確正切曲線繪3次樣條插值正切曲線繪精確正弦曲線繪3次樣條拯值正弦曲線繪5次多項(xiàng)式正切插值曲線繪5次多項(xiàng)式正弦插值曲線開(kāi)始(2)代碼：clearall;closeall;clc;X=1,4,9,16,25,36,49,64,81,100;Y=1:10;X1=1:100;Y1=interp1(X,Y,X1,cublc);plot(X,Y,r:o);holdon;plot(X1,Y11k-x);legend3、已知一組實(shí)驗(yàn)數(shù)據(jù)如表6.6所示。表6.6一組實(shí)驗(yàn)數(shù)據(jù)ii2345x165123150123141y.187126172125148求它的線性擬合曲線代碼：clearall;closeall;clc;x=165,123,150,123,141;y