应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助matlab求解三元超越方程
51/1返回列表
查看: 1280  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wangcy1989

铁虫 (初入文坛)

[求助] 求助matlab求解三元超越方程

本人因为刚注册不久，金币很少，悬赏金币也较少，请您谅解，但我可以现金酬谢，我的邮箱wangcy0929@163.com，比较着急，谢谢了！
1、下面的图片是我要求的三元超越方程。
求助matlab求解三元超越方程
2、这是我自己写的代码，但程序已经运行一整天了，既没有报错，也没有运行出结果，请各位大牛指教，不甚感激！
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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
止于至善，严谨求实；诚朴雄伟，厚德载物；博学明辨，慎思笃行。
1楼 2013-09-25 10:10:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangcy1989

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by pl408367247 at 2013-09-26 12:24:56
由于计算机性能问题，我只运算了a=的 ，运算有没有结果你看matlab中的workspace中的变量c

我找到同学帮忙刚画出来了，还是非常感谢您的热心帮助！
赞一下 回复此楼
止于至善，严谨求实；诚朴雄伟，厚德载物；博学明辨，慎思笃行。
5楼2013-09-26 15:10:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

pl408367247

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dingd: 金币+1, 谢谢交流。 2013-09-25 22:15:30
wangcy1989: 金币+5, ★有帮助, 非常感谢您的热心帮助，金币有限，全部给您了~ 2013-09-26 15:11:31
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
赞一下 回复此楼
学习中！
2楼2013-09-25 21:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangcy1989

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by pl408367247 at 2013-09-25 21:35:28
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* ...

我运行了一下，没有结果。
另外，我想要的是画出h与a,b的三维图形
赞一下 回复此楼
止于至善，严谨求实；诚朴雄伟，厚德载物；博学明辨，慎思笃行。
3楼2013-09-26 11:20:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pl408367247

铁虫 (小有名气)

【答案】应助回帖

★
csgt0: 金币+1, 谢谢 2013-09-26 13:19:59
由于计算机性能问题，我只运算了a=[0，0.1]的 ，运算有没有结果你看matlab中的workspace中的变量c
赞一下 回复此楼
学习中！
4楼2013-09-26 12:24:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭