| 5 | 1/1 | 返回列表 |
| 查看: 1320 | 回復(fù): 4 | |||
wangcy1989鐵蟲 (初入文壇)
|
[求助]
求助matlab求解三元超越方程
|
|
本人因?yàn)閯傋?cè)不久,金幣很少,懸賞金幣也較少,請(qǐng)您諒解,但我可以現(xiàn)金酬謝,我的郵箱wangcy0929@163.com,比較著急,謝謝了! 1、下面的圖片是我要求的三元超越方程。 2、這是我自己寫的代碼,但程序已經(jīng)運(yùn)行一整天了,既沒有報(bào)錯(cuò),也沒有運(yùn)行出結(jié)果,請(qǐng)各位大牛指教,不甚感激! alpha=0.1;bata=0.1;r=1; M=0.5;TA=0.2; A=10*M; k=0.5;l=0.9; syms a b h X; X=1.5; g1=r*(alpha+bata*(r*h+A-M)-TA).^a/(2*A.*(a+2)); g2=(1+k*l).*(a+1)-k*(alpha+bata*(r*h+A-M)-TA)/(A*bata); g3=k*X*r*((TA-alpha).^(b+1))/(2*A^2*bata.*(b+1)); f=@(a,b,h)g1*g2-g3-h; [a,b,h]=meshgrid(linspace(0,1)); val=f(a,b,h); [p,v]=isosurface(a,b,h,val,0); patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat'); view(3);grid on; axis equal [ Last edited by wangcy1989 on 2013-9-25 at 10:22 ] |

鐵蟲 (小有名氣)
|
clc; clear; alpha=0.1;bata=0.1;r=1; M=0.5;TA=0.2; A=10*M; k=0.5;l=0.9; X=1.5; syms h; a=0:0.01:1-0.01; b=a; N=length(a); c=zeros(N,N); for i=1:10 g1=r*(alpha+bata*(r*h+A-M)-TA)^a(i)/(2*A*(a(i)+2)); g2=(1+k*l)*(a(i)+1)-k*(alpha+bata*(r*h+A-M)-TA)/(A*bata); for j=i+1:11 g3=k*X*r*((TA-alpha)^(b(j)+1))/(2*A^2*bata*(b(j)+1))-h; d=solve(g1*g2-g3,'h'); N_d=length(d); for mm=1:N_d %從眾多解中選合適解 if(d(mm)==real(d(mm))) gg=g1*g2-g3; v_gg=subs(gg,d(mm)); if(double(abs(v_gg))<0.0000001) break; end end end c(i,j)=d(mm);%c就是你要求的h end end |

鐵蟲 (初入文壇)

鐵蟲 (小有名氣)

鐵蟲 (初入文壇)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +6 | 路癡小琪 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +4 | 新時(shí)代材料 2026-03-27 | 4/200 |
|
|
[考研] 311求調(diào)劑 +7 | lin0039 2026-03-26 | 7/350 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 | 4/200 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研]
|
y7czhao 2026-03-26 | 8/400 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 一志愿南大,0703化學(xué),分?jǐn)?shù)336,求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|