作业帮 > 数学 > 作业

matlab 求解带参数方程组

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/30 07:00:48
matlab 求解带参数方程组
在matlab里求解下面的方程,得不到解.
f1(a,b)=-1/(2*b)*(a^3/((a^2+b^2)^1.5)-1/((1+b^2)^1.5);
f2(a,b)=-1/(24*b^3)*(a^3*(20b^4+7*b^2*a^2+2*a^4)/((a^2+b^2)^3.5)-(20*b^4+7*b^2+2)/((1+b^2)^3.5));
d用a,b表示
方程:a*f1(a,b)-c*f1(c,d)=0;
a^3*f2(a,b)-c^3*f2(c,d)=0
我自己用下式求解,说第一个方程is not a valid expression or equation.
[c,d]=solve('-a/(2*b)*(a^3/((a^2+b^2)^1.5)-1/((1+b^2)^1.5)-(c/(2*d)*(c^3/((c^2+d^2)^1.5)-1/((1+d^2)^1.5)))=0','-a^3/(24*b^3)*(a^3*(20*b^4+7*b^2*a^2+2*a^4)/((a^2+b^2)^3.5)-(20*b^4+7*b^2+2)/((1+b^2)^3.5))-(c^3/(24*d^3)*(c^3*(20*d^4+7*d^2*c^2+2*c^4)/((c^2+d^2)^3.5)-(20*d^4+7*d^2+2)/((1+d^2)^3.5))=0')
syms a b c d
f1ab=-1/(2*b)*(a^3/((a^2+b^2)^1.5)-1/((1+b^2)^1.5));
f2ab=-1/(24*b^3)*(a^3*(20*b^4+7*b^2*a^2+2*a^4)/((a^2+b^2)^3.5)-(20*b^4+7*b^2+2)/((1+b^2)^3.5));
f1cd=subs(subs(f1ab,a,c),b,d);
f2cd=subs(subs(f2ab,a,c),b,d);
p=a*f1ab-c*f1cd;
q=a^3*f2ab-c^3*f2cd;
[cc,dd]=solve(p,q,'c,d')
首先很遗憾没有解析解.你的表达式太复杂了.
但是可以告诉你的是方法是可行的.
不信,你可以将f1ab,f2ab改成下面的试一下:
f1ab=a+b;
f2ab=a-b;
我记起来了,好像我对你的公式里面有修改了一下,你再仔细地看你的公式,和我上面的吧.我记得我加了一两个括号