版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

淺紫微藍(lán)

木蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1684.5
散金: 50
帖子: 6
在線: 8.4小時
蟲號: 1352694
注冊: 2011-07-21
性別: MM
專業(yè): 無機(jī)化工

[求助] MATLAB求解多元非線性方程組

我想求解一個非常復(fù)雜的方程組，原理很簡單，就是解三個方程求出三個未知數(shù)a、b、c，但每一個方程都包含exp、log、sqrt函數(shù)，不知道能否用MATLAB求出解來，我試著寫了程序，但是運行不出來，我第一次用MATLAB，不知道寫的對不對，請各位高手幫忙解出該方程組來，我有急用，謝謝了~
syms a b c
I=3*a+4*b+3*c;
f=-4*0.377/1.2*I*log(1+1.2*sqrt(I));
f1=-2*0.377*(sqrt(I)/(1+1.2*sqrt(I))+2/1.2*log(1+1.2*sqrt(I)));
B1=-0.03942+0.97225/2/I*(1-(1+2*sqrt(I))*exp(-2*sqrt(I)));
B11=0.97225/2/I/I*(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I)));
C1=0.0179/sqrt(8);
B2=0.23825+2.9605/0.98/I*(1-(1+1.4*sqrt(I))*exp(-1.4*sqrt(I)))-30.905/72/I*(1-(1+12*sqrt(I))*exp(-12*sqrt(I)));
B21=2.9605/0.98/I/I*(-1+(1+1.4*sqrt(I)+0.98*I)*exp(-1.4*sqrt(I)))-30.905/72/I/I*(-1+(1+12*sqrt(I)+72*I)*exp(-12*sqrt(I)));
C2=0.003;
B3=0.00488+0.4911/2/I*(1-(1+2*sqrt(I))*exp(-2*sqrt(I)));
B31=0.4911/2/I/I*(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I)));
C3=-0.00116/sqrt(8);
g=(1-(1+2*sqrt(I))*exp(-2*sqrt(I)))/2/I;
g1=(-1+(1+2*sqrt(I)+2*I)*exp(-2*sqrt(I)))/2/I/I;
w=1+(I*f1-f+4*a*(a+b+c)*(B1+I*B11+4*(a+b+c)*C1)+2*b*(a+b+c)*(B2+I*B21+4*(a+b+c)*C2)+4*c*(a+b+c)*(B3+I*B31+4*(a+b+c)*C3))/(3*a+2*b+3*c);

r1=exp(f1+(4/3*a+4/3*(a+b+c))*(B1+2*(a+b+c)*C1)+4/3*c*(B3+2*(a+b+c)*C3)+2/3*b*(B2+2*(a+b+c)*C2)+2*a*(a+b+c)*(2*B11+4/3*C1)+2*c*(a+b+c)*(2*B31+4/3*C3)+b*(a+b+c)*(2*B21+4/3*C2)+8/3*c*(0.125167+1.078184*g)+4/3*b*(2.521475-15.811274*g)+(4/3*c*(a+b+c)+4/6*a*c)*(-0.0080597)+(2/3*b*(a+b+c)+1/3*a*b)*(-1.016358)+1/3*c*b*0.0672253+4*a*c*1.078184*g1+2*a*b*(-15.811274)*g1+2*c*b*10.163009*g1);

r2=exp(2*f1+(b+(a+b+c))*(B2+2*(a+b+c)*C2)+2*a*(B1+2*(a+b+c)*C1)+2*c*(B3+2*(a+b+c)*C3)+2*a*(a+b+c)*(4*B11+2*C1)+b*(a+b+c)*(4*B21+2*C2)+2*c*(a+b+c)*(4*B31+2*C3)+2*a*(2.521475-15.811274*g)+2*c*(-0.0692694+10.163009*g)+(a*(a+b+c)+1/2*a*b)*(-1.016358)+(c*(a+b+c)+1/2*b*c)*0.0672253+a*c*(-0.0080597)+4*a*b*(-15.811274)*g1+4*c*b*10.163009*g1+8*a*c*1.078184*g1);

r3=exp(f1+(4/3*c+4/3*(a+b+c))*(B3+2*(a+b+c)*C3)+4/3*a*(B1+2*(a+b+c)*C1)+2/3*b*(B2+2*(a+b+c)*C2)+2*a*(a+b+c)*(2*B11+4/3*C1)+2*c*(a+b+c)*(2*B31+4/3*C3)+b*(a+b+c)*(2*B21+4/3*C2)+8/3*a*(0.125167+1.078184*g)+4/3*b*(-0.0692694+10.163009*g)+(4/3*a*(a+b+c)+4/6*a*c)*(-0.0080597)+(2/3*b*(a+b+c)+1/3*c*b)*0.0672253+1/3*a*b*(-1.016358)+4*a*c*1.078184*g1+2*a*b*(-15.811274)*g1+2*c*b*10.163009*g1);

eq1=b*r2*(exp(-(3*a+2*b+3*c)*w/55.51))^7-0.076849;
eq2=b*c*r2*r3*(exp(-(3*a+2*b+3*c)*w/55.51))^6-0.012198;
eq3=a*r1*(exp(-(3*a+2*b+3*c)*w/55.51))^10-0.092085;
[a,b,c]=solve(eq1,eq2,eq3,'a,b,c')

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)

有些人離開就是離開了，漸漸地，生活會變得沒有什么不同，仿佛那個人不是消失了，而是從未曾出現(xiàn)過。這是我們所希望的，也是必需承認(rèn)的，原來我們沒有那么重要

1樓 2013-03-22 10:39:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

matlab解非線性方程組解不開已經(jīng)有3人回復(fù)
matlab求解非線性方程組已經(jīng)有3人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
請教用matlab求解一個非線性偏微分方程組的數(shù)值解已經(jīng)有8人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
matlab非線性方程組求解，急�。�！已經(jīng)有5人回復(fù)
用不動點迭代求非線性方程組，求得的不是想要的解，怎么辦已經(jīng)有12人回復(fù)
初學(xué)者，請教大家關(guān)于matlab解多元二次方程組已經(jīng)有4人回復(fù)
matlab求解非線性方程組，求助！已經(jīng)有6人回復(fù)
如何使用matlab求解非線性方程組的所有整數(shù)解？已經(jīng)有9人回復(fù)
Matlab求解偏微分方程組已經(jīng)有22人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
請教Matlab高手解多元非線性方程組中的參數(shù) 已經(jīng)有6人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
MATLAB非線性方程組求解~~~反應(yīng)釜模擬已經(jīng)有9人回復(fù)
【求助】多元方程組求解及擬合已經(jīng)有9人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
非線性方程組的迭代法（數(shù)值計算高手請進(jìn)）已經(jīng)有7人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主淺紫微藍(lán) 的主題更新