| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1319 | 回復(fù): 4 | ||
wangcy1989鐵蟲(chóng) (初入文壇)
|
[求助]
求助matlab求解三元超越方程
|
|
本人因?yàn)閯傋?cè)不久,金幣很少,懸賞金幣也較少,請(qǐng)您諒解,但我可以現(xiàn)金酬謝,我的郵箱wangcy0929@163.com,比較著急,謝謝了! 1、下面的圖片是我要求的三元超越方程。 2、這是我自己寫(xiě)的代碼,但程序已經(jīng)運(yùn)行一整天了,既沒(méi)有報(bào)錯(cuò),也沒(méi)有運(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 ] |

鐵蟲(chóng) (小有名氣)
|
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 |

鐵蟲(chóng) (初入文壇)

鐵蟲(chóng) (小有名氣)

鐵蟲(chóng) (初入文壇)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +9 | 懷瑾握瑜l 2026-03-20 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +4 | 13831862839 2026-03-24 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿武理085500機(jī)械專(zhuān)業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 08工學(xué)調(diào)劑 +12 | 用戶(hù)573181 2026-03-20 | 17/850 |
|
|
[考研] 一志愿南航材料專(zhuān)317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[考博] 申博26年 +4 | 八6八68 2026-03-19 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
[考研] 材料專(zhuān)碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專(zhuān)材355 2026-03-19 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|