全文預(yù)覽已結(jié)束
下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
沙石形骨料混凝土的數(shù)學(xué)模型與計(jì)算機(jī)實(shí)現(xiàn) 宋來(lái)忠設(shè)有沙石形骨料混凝土塊,長(zhǎng)、寬、高分別為a、b、c(單位: cm), 其中有n塊骨料,記為Pi, i=1,n,每塊骨料的最大直徑分別為2ri(單位: cm),i=1,n,。rie, f,不失一般性,建立如圖1所示的坐標(biāo)系,可假定骨料Pi剛好位于球心在(xi, yi, zi)半徑為ri的球Si中。(優(yōu)化) 假如隨機(jī)的用N塊骨料與水泥細(xì)沙攪拌成混凝土泥漿,制作長(zhǎng)、寬、高分別為a、b、c混凝土塊,用料中含骨料n,(n塊都放進(jìn),所以n要適當(dāng),可以由程序1實(shí)驗(yàn)得到),此時(shí),上述半徑為ri的球Si的球心坐標(biāo)(xi, yi, zi)滿足(1)、(2),有優(yōu)化模型:約束函數(shù):s. t 當(dāng)n不太大時(shí),比如小于100;求解較為可行;程序2是取n=60運(yùn)行的情形:function c,ceg=shuiliY(x)global n r n=60;load r;b=15;d=15;f=15;N=1;for i=1:n-1 for j=i+1:n A(N)=(x(i)-x(j)2+(x(n+i)-x(n+j)2+(x(2*n+i)-x(2*n+j)2-(r(i)+r(j)2;N=N+1; endendA1=A(:);%Nc=-A1;-x+r;r;r;(x-b)+r;r;r; ceg=;目標(biāo)函數(shù):function f=shuiliM(x)global n rn=60;s=0;%r=unifrnd(0.02,0.2,n,1);b=15;d=15;f=15;x=x;for i=1:n for j=1:n %i+1:n s=s+(x(i)-x(j)2+(x(n+i)-x(n+j)2+(x(2*n+i)-x(2*n+j)2;ends=s+x(i)2+(x(i)-b)2+x(n+i)2+(x(n+i)-d)2+x(2*n+i)2+(x(2*n+i)-f)2;endf=s;主程序:clear;close;global n r n=60;%r=unifrnd(0.2,2,n,1);x0=5*ones(3*n,1); %b=15;d=15;f=15;%VLB=zeros(3*n,1)+0.2;VUB=15*ones(3*n,1)-0.2;load r;VLB=zeros(3*n,1);VUB=15*ones(3*n,1);Aeg=;beg=;OPT=optimset; OPT.LargeScale=off; x,fval=fmincon(shuiliM,x0,Aeg,beg,VLB,VUB,shuiliY,OPT);X=x(1:n);Y=x(n+1:2*n);Z=x(2*n+1:3*n);X(:),Y(:),Z(:),r(:)size(X(:),Y(:),Z(:),r(:)for i=1:n k = 5;m = 2k-1;theta = pi*(-m:2:m+2)/m;phi = (pi/2)*(-m:2:m)/m;s1=unifrnd(0.3,0.9);s2=unifrnd(0.3,0.9);s3=unifrnd(0.3,0.9);X1 =s1*r(i)*cos(phi)*cos(theta)+X(i);Y1 =s2*r(i)*cos(phi)*sin(theta)+Y(i);Z1 =s3*r(i)*sin(phi)*ones(size(theta)+Z(i);colormap(pink); C = hadamard(2k); surf(X1,Y1,Z1,C); axis(0 15 0 15 0 15)hold on;end當(dāng)n較大時(shí),比如大于100;求解較為復(fù)雜,一般用時(shí)(計(jì)算機(jī)運(yùn)行)較長(zhǎng),而在工程實(shí)際中n往往是很大的。因此,有必要用一些非標(biāo)準(zhǔn)的方法求解,比如遺傳算法、模
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2024物業(yè)管理服務(wù)協(xié)議:智能住宅小區(qū)安全監(jiān)控合同3篇
- 2024年標(biāo)準(zhǔn)社區(qū)安保協(xié)議終止書(shū)版B版
- 2024年知識(shí)產(chǎn)權(quán)許可合同:專利權(quán)人與被許可人之間的專利使用
- 2024年隱名股東權(quán)益分配合同版B版
- 2024食用油包裝設(shè)計(jì)及印刷服務(wù)合同3篇
- 2024建筑工程施工管理與質(zhì)量安全保障合同
- 2025年度電影劇本創(chuàng)作編劇助理及現(xiàn)場(chǎng)工作合同3篇
- 2024年鐘點(diǎn)工雇傭合同3篇
- 2024年重點(diǎn)交通樞紐土方運(yùn)輸工程承包合同書(shū)范本3篇
- 2024年智能消防系統(tǒng)研發(fā)與實(shí)施合同3篇
- 《皮膚病中成藥導(dǎo)引》課件
- Chinese Tea Culture.中國(guó)茶文化(課堂PPT)
- 排洪溝工程設(shè)計(jì)說(shuō)明
- (整理版高中英語(yǔ))杭高第一學(xué)期期末考試高一英語(yǔ)試卷
- 學(xué)校食堂生鮮肉(豬肉、牛肉、羊肉、雞鴨鵝肉)配送服務(wù)方案
- SMT工資方案(原創(chuàng))
- 關(guān)于礦棉裝飾吸聲板檢驗(yàn)報(bào)告加圖標(biāo)版
- 大紅色節(jié)word感謝信信紙背景模板
- 安全檢查匯報(bào)材料
- 2005年海南高考理科綜合真題及答案
- 機(jī)房巡檢記錄表.doc
評(píng)論
0/150
提交評(píng)論