Matlab編程和模擬模型1?編程訓(xùn)練題目一?有一個(gè)向量V二加電％，編程求解下列式子:n2iJS二送3ii72i2 | , 2S2=ai?83M../au其中U為不超過(guò)n的最大奇數(shù)。自己給出測(cè)試?yán)樱⒔o出運(yùn)算結(jié)果二?模擬一次撲克牌洗牌過(guò)程，輸出牌的順序%4*13+2paihao-[11:23,31:43,51:63,71:83,98,99];iflength(paihao)^=54,error(*errorinnumbenendPAIXU二randperm(54);fapai=paihao(PAIXU)三?模擬1000次投幣過(guò)程，輸出出現(xiàn)正面的頻率。程序1:num_zheng=0;fori=1:1000,cur=rand;ifcur<0.5,num_zheng二num_zheng+1;endendf=num_zheng/1000四?問(wèn)題：現(xiàn)有一個(gè)文本文件(mydata?txt),包含了火車(chē)的出發(fā)和到達(dá)時(shí)刻，基本格式如下(實(shí)際中該文件超過(guò)1000行)19:3507:4719:1407:1219:2107:1918:5308:2319:0707:0519:0006:5821:3708:0420:4007:4021:4407:2221:3007:2519:2806:4822:0007:1019:3507:4719:1407:1219:2107:1918:5308:2319:0707:0519:0006:5821:3708:0420:4007:4021:4407:2221:3007:2519:2806:4822:0007:1019:2107:1918:5308:2319:0707:05要求把時(shí)間讀取并存為矩陣數(shù)據(jù)，并且需要將其以始0:00時(shí)刻記為初時(shí)刻，把奇數(shù)行數(shù)據(jù)轉(zhuǎn)為分鐘。準(zhǔn)備知識(shí)：文件打開(kāi)(fopen)，尖閉(fclose);讀取各行數(shù)據(jù)(fgetl)(3)解析各行數(shù)據(jù)，并轉(zhuǎn)換數(shù)據(jù)注：處理方法不唯一。文本文件：fopen，fclose，fgetl。%先初始化矩陣mat=[];row=0;fid=fopen(,mydata.txt,);while1tline=fgetl(fid);iffeof(fid),break,enddisp(tline)cur=str2num(tline);row=row+1;%將當(dāng)前行的(行)向量保存到mat中

ifisempty(mat),mat-cur;elselen=length(cur);mat(row,1:len)二cur;endendfclose(fid);19357478319421914712000019217190000185382300001977500001906580000213784000020407400000214472200002130725000019286480000五.(適中)有一個(gè)文件mytest?txt,內(nèi)容如下:%請(qǐng)將此文件中的數(shù)值讀取為1個(gè)向量1 4 7 2 5 83 6 9特點(diǎn)：第2行開(kāi)始為數(shù)值數(shù)據(jù)。請(qǐng)編程將這個(gè)文件中從第2行開(kāi)始的所有數(shù)保存為1個(gè)列向量(按行先后順序排列這些數(shù))程序%先初始化矩陣vec=[];tline=fgetl(fid);%跳過(guò)第1行字符串，不處理fid=fopen(,mytest.txt,);while1tline=fgetl(fid);iffeof(fid),break,enddisp(tline)cur=str2num(tline);%將當(dāng)前行的（行）向量保存到vec中vec二[vec,cur];endfclose(fid);vec二vec'%請(qǐng)將此文件中的數(shù)值讀取為1個(gè)向量1 4 7 2 5 83 6 9vec=147258369L' JX1X2+X3+X：y=a十~2六?（較難）有一個(gè)函數(shù)人如下。X5X3 ,另有一個(gè)函數(shù)論）定義f(X)丄X2」XsjX4JXs「 L、L、LL、CX1 GX2CX3CX4CX5O編寫(xiě)函數(shù)完成計(jì)算函數(shù)f（x）在任意點(diǎn)的函數(shù)值的功能1? 創(chuàng)建符號(hào)變量X,編寫(xiě)符號(hào)表達(dá)式y(tǒng)；symsx1x2x3x4X5e1-sqrt(x1*x2+x3+x4八2)/(x5+x3八2);y=sqrt(1+e1);利用diff計(jì)算偏導(dǎo)并創(chuàng)建f(x)的表達(dá)式funx=diff(y,x1)*x1+ diff(y,x2)*x2+diff(y,x3)*x3+diff(y,x4)*x4+diff(y,x5)*x5;復(fù)制funx表達(dá)式文本、并粘貼到一個(gè)程序文件在命令行輸入myfun(1:5)檢查語(yǔ)法是否正確%diff(exp,x)%exp為符號(hào)表達(dá)式2模擬模型案例：修理廠的模擬某修理廠設(shè)有3個(gè)停車(chē)位置，其中一個(gè)位置供正在修理的汽車(chē)停放?現(xiàn)以一天為一個(gè)時(shí)段，每天到達(dá)修理站的汽車(chē)數(shù)有如下概率分布：到達(dá)數(shù)X012概率0.30.50.2假定在一個(gè)時(shí)段內(nèi)能修好汽車(chē)數(shù)有如下概率分布:修好數(shù)丫012概率0.250.70.05本時(shí)段內(nèi)未能完成修理的汽車(chē)與正在等待修理的汽車(chē)一起進(jìn)入下一時(shí)段 ?試問(wèn):該停車(chē)廠有無(wú)必要增加停車(chē)位置？并說(shuō)明理由.3?已知非線(xiàn)性整數(shù)規(guī)劃為：22222Maxz=%x23x34x42x5-8Xi?2x??3X3?X4?2Xs0<\<99 5)為+X2十X3十X4+X5蘭400s.t.*捲+2X2+2X3+X4+6X5蘭8002Xj+X2+6X3蘭200X3X45X5乞200對(duì)該題，目前尚無(wú)有效方法求出準(zhǔn)確解。 如果用顯枚舉法試探，共需計(jì)算（100）5=1010個(gè)點(diǎn)，其計(jì)算量非常之大。然而應(yīng)用蒙特卡洛去隨機(jī)計(jì)算 1個(gè)點(diǎn)，便可找到滿(mǎn)意解，那么這種方法的可信度究竟怎樣呢？下面就分析隨機(jī)取樣采集106個(gè)點(diǎn)計(jì)算時(shí)，應(yīng)用概率理論來(lái)估計(jì)一下可信度。不失一般性,假定一個(gè)整數(shù)規(guī)劃的最優(yōu)點(diǎn)不是孤立的奇點(diǎn)。假設(shè)目標(biāo)函數(shù)落在高值區(qū)的概率分別為 0.01,0.00001，則當(dāng)計(jì)算106個(gè)點(diǎn)后,有任一個(gè)點(diǎn)能落在高值區(qū)的概率分別為1_o.991000000,0.99 99（100多位），1-O.999991000000:0.999954602。

解(i)首先編寫(xiě)M文件mente.m定義目標(biāo)函數(shù)f和約束向量函數(shù)g，程序如下:[f,g]=mengte(x);f=x(1)a2+x⑵人2+3妝⑶a2+4*x(4)A2+2*x⑸?8*x(1)-2*x(2)-3*x(3) ????x(4)-2*x(5);g(1)=sum(x)-400;g(2)=x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-800;g(3)=2*x(1)+x(2)+6*x(3)-200;g(4)=x(3)+x(4)+5*x(5)-200;(ii)編寫(xiě)如下程序求問(wèn)題的解：rand('state',sum(clock));p0=0;ticfori=1:10A5x=99*rand(5,1);x1=floor(x)