作业帮 > 综合 > 作业

求matlab高手帮忙解一个二元方程组,

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/04/26 01:06:06
求matlab高手帮忙解一个二元方程组,
有两个方程
0.015=d+p
0.050=0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5
求高手帮忙用matlab求解出d和p
方程无解.可以很简单地画出两个方程的曲线来:ezplot('-0.050+0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5',[-0.02 0.02])
hold on
ezplot('-0.015+d+p',[-0.02 0.02])由图可见,两条线无交点,也就意味着方程组无解. 或者换另一种方式来分析:把方程1代入方程2:syms d p
eq1=-0.015+d+p;
eq2=-0.050+0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5;
D=solve(eq1,d)
eq=subs(eq2,d,D)得到方程:>> eq=subs(eq2,d,D) 
eq = 
7/200 - 2*(129/400000 - (43*p)/2000)^(1/2) - 2*((63*p)/1250)^(1/2)方程中两个包含p的项都要进行开平方,应该大于0,即:(63*p)/1250 >= 0129/400000 - (43*p)/2000 >= 0联立解得 0<=p<=0.015,画出该方程的曲线如下:不可能与0轴相交,也即方程无解.
再问: 不好意思,有可能是单位的原因,我现在把数值乘以1000,您看看能求出p和d吗,p和d都大于0 15=d+p 50=15+70-2.*(d.*21.5).^0.5-2.*(p.*50.4).^0.5 又加了100分,呵呵
再答: 用上面相同的方法,可以画出曲线来:ezplot('-15+d+p',[-1 1]*20)
hold on
ezplot('-50+15+70-2.*(d.*21.5).^0.5-2.*(p.*50.4).^0.5',[-1 1]*20)由图可见,两条曲线仍然没有交点。 同样,把方程1代入方程2:syms d p
eq1=-15+d+p;
eq2=-50+15+70-2.*(d.*21.5).^0.5-2.*(p.*50.4).^0.5;
D=solve(eq1,d)
eq=subs(eq2,d,D)
ezplot(eq,[-1 16]) 在p的有效取值范围内,仍然不可能和0轴相交(也即无解)。 请你再仔细看一下方程有没有其他问题。