24小时热门版块排行榜

返回列表

wangcy1989

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 69
帖子: 5
在线: 7.6小时
虫号: 2173681
注册: 2012-12-08
性别: GG
专业: 信息经济学

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

本人因为刚注册不久，金币很少，悬赏金币也较少，请您谅解，但我可以现金酬谢，我的邮箱wangcy0929@163.com，比较着急，谢谢了！
1、下面的图片是我要求的三元超越方程。

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人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复
物理学I论文润色/翻译怎么收费? 已经有182人回复
大豆异黄酮分离已经有0人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有4人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用matlab如何解方程？已经有6人回复
求Matlab高手解决线性方程组的迭代求解问题已经有17人回复
如何求解下列方程：xlnx=a(x-b) 已经有14人回复
各位朋友，如何用matlab编写这个方程，并求解？已经有15人回复
求助用matlab拟合三元线性回归方程及分析已经有9人回复
matlab求解方程中的参数已经有21人回复
求助matlab---fsolve解非线性方程组已经有6人回复
MATLAB用于三元回归参数已经有3人回复
Matlab求解偏微分方程组已经有22人回复
matlab求解非线性方程组已经有16人回复
matlab解微分方程组已经有15人回复
求matlab求解一个方程组的问题已经有6人回复
MATLAB解方程出现??? Subscripted assignment dimension mismatch. 已经有3人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】matlab 解矩阵方程重金悬赏解决方案已经有10人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】三元二次方程求最大值已经有3人回复
求助一个差分方程的MATLAB实现已经有3人回复
matlab求解超越方程已经有3人回复
【求助】求教matlab解非线性方程组已经有9人回复

止于至善，严谨求实；诚朴雄伟，厚德载物；博学明辨，慎思笃行。

1楼 2013-09-25 10:10:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

pl408367247

铁虫 (小有名气)

应助: 3 (幼儿园)
金币: 35.4
红花: 1
帖子: 105
在线: 75.6小时
虫号: 1402271
注册: 2011-09-15
性别: GG
专业: 信号理论与信号处理

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 69
帖子: 5
在线: 7.6小时
虫号: 2173681
注册: 2012-12-08
性别: GG
专业: 信息经济学

引用回帖:

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的回帖