| 查看: 1547 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助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 ] |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有240人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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人回复
pl408367247
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 36.9
- 红花: 1
- 帖子: 105
- 在线: 75.6小时
- 虫号: 1402271
- 注册: 2011-09-15
- 性别: GG
- 专业: 信号理论与信号处理
4楼2013-09-26 12:24:56
pl408367247
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 36.9
- 红花: 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
感谢参与，应助指数 +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
3楼2013-09-26 11:20:06
5楼2013-09-26 15:10:34
