24小时热门版块排行榜    

CyRhmU.jpeg
查看: 948  |  回复: 4

mmffdd26

金虫 (著名写手)

[求助] 求matlab高手帮忙求解一个含xyz的方程,求出y=f(x,z)~~~~~~~~~~~~~~~~

求matlab高手帮忙求解一个含xyz的方程,求出y=f(x,z)~~~~~~~~~~~~~~~~
回复此楼

» 猜你喜欢

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

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

mmffdd26

金虫 (著名写手)

方程为
1./2.*(820000.*z.*9./0.05415).*[y.^2.*(x-y./3)-4./3.*((2000000-1000000)./(1820000-1000000)).*x.^3-(((2.*x-y).*y).^0.5).^2.*((-2.*((((2.*x-y).*y).^0.5)./x).*((((2.*x-y).*y).^0.5)./x.*cos(5.*pi./6)+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).*sin(5.*pi./6)).*0.05415./(1820000.*x.^2.*z.*9)-(4./3).*2000000./1820000+(2./3-((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)+(1./3).*((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).^3)+1000000./1820000.*(2./3+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)-(1./3).*(((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)).^3)).*x./[(1000000./1820000-1).*((((2.*x-y).*y).^0.5)./x).^2])]-((((2.*x-y).*y).^0.5).*sin(acos(y./x-1)-(5.*pi./6)))=0



方程里面含x,y,z三个变量,我想算出y=f(x,z) 这种形式

但是用solve说方程复杂了,算不出结果,方程如下



syms x y z

solve(1./2.*(820000.*z.*9./0.05415).*[y.^2.*(x-y./3)-4./3.*((2000000-1000000)./(1820000-1000000)).*x.^3-(((2.*x-y).*y).^0.5).^2.*((-2.*((((2.*x-y).*y).^0.5)./x).*((((2.*x-y).*y).^0.5)./x.*cos(5.*pi./6)+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).*sin(5.*pi./6)).*0.05415./(1820000.*x.^2.*z.*9)-(4./3).*2000000./1820000+(2./3-((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)+(1./3).*((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).^3)+1000000./1820000.*(2./3+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)-(1./3).*(((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)).^3)).*x./[(1000000./1820000-1).*((((2.*x-y).*y).^0.5)./x).^2])]-((((2.*x-y).*y).^0.5).*sin(acos(y./x-1)-(5.*pi./6)))==0,y,x>0,z>0)



好像可以算出z=f(x,y),但我想算出y=f(x,z) ,求高手指点怎么列程序算出y=f(x,z)

x y z


可以把限制在一定区间内,看能得到结果不
x区间 0,0.000001

y区间 0,0.000001

z区间 0,50000
2楼2013-06-21 16:48:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mmffdd26

金虫 (著名写手)

刚打错了,我算的是下面这个,没 x>0,z>0
syms x y z
solve(1./2.*(820000.*z.*9./0.05415).*[y.^2.*(x-y./3)-4./3.*((2000000-1000000)./(1820000-1000000)).*x.^3-(((2.*x-y).*y).^0.5).^2.*((-2.*((((2.*x-y).*y).^0.5)./x).*((((2.*x-y).*y).^0.5)./x.*cos(5.*pi./6)+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).*sin(5.*pi./6)).*0.05415./(1820000.*x.^2.*z.*9)-(4./3).*2000000./1820000+(2./3-((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)+(1./3).*((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x).^3)+1000000./1820000.*(2./3+((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)-(1./3).*(((x.^2-(((2.*x-y).*y).^0.5).^2).^0.5./x)).^3)).*x./[(1000000./1820000-1).*((((2.*x-y).*y).^0.5)./x).^2])]-((((2.*x-y).*y).^0.5).*sin(acos(y./x-1)-(5.*pi./6)))==0,y)
3楼2013-06-21 16:51:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mmffdd26

金虫 (著名写手)

高手在哪啊
4楼2013-06-21 21:30:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyandaisha

金虫 (正式写手)

看起来好复杂的样子
一直在努力
5楼2013-06-23 18:28:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mmffdd26 的主题更新
信息提示
请填处理意见