怎么用c语言解四元方程,求四元一次方程组的解

满意答案

哎ya088

2013.05.14

采纳率：53%    等级：12

已帮助：6014人

用matlab计算就简单了，指令：

y=solve('x^3*a+x^2*b+x*c+d=y1','(x+t)^3*a+(x+t)^2*b+(x+t)*c+d=y2','(x+2*t)^3*a+(x+2*t)^2*b+(x+2*t)*c+d=y3','(x+3*t)^3*a+(x+3*t)^2*b+(x+3*t)*c+d=y4','a,b,c,d')

输出：

y =

a: [1x1 sym]

b: [1x1 sym]

c: [1x1 sym]

d: [1x1 sym]

再指令：

[y.a;y.b;y.c;y.d]

ans =

[ -1/6*(y1-3*y2+3*y3-y4)/t^3]

[ 1/2*(x*y1-3*x*y2+3*x*y3-x*y4+2*y1*t-5*y2*t+4*y3*t-y4*t)/t^3]

[ -1/6*(-30*x*t*y2+24*x*t*y3-6*x*t*y4+12*x*t*y1+11*y1*t^2-18*y2*t^2+9*y3*t^2-2*y4*t^2+3*x^2*y1-9*x^2*y2+9*x^2*y3-3*x^2*y4)/t^3]

[ 1/6*(12*x^2*y3*t-3*x^2*y4*t+9*x*y3*t^2-18*x*y2*t^2+11*x*y1*t^2-15*x^2*y2*t+6*x^2*y1*t-2*x*y4*t^2+6*y1*t^3+x^3*y1-3*x^3*y2+3*x^3*y3-x^3*y4)/t^3]追问： 非常感谢！本人也稍微会点matlab，上面也大致上看懂了，可否对您的两个指令和两个结果稍作介绍？先增加50个悬赏表示对您回答的感谢！

追答：y=solve('方程1','方程2',...,'变量1,变量2,...') %这是求解代数方程组的指令

但输出结果是如下形式：

y =

a: [1x1 sym]

b: [1x1 sym]

c: [1x1 sym]

d: [1x1 sym]

需要再追加指令使其显示出来：

y.a %这是显示a的。即y 后跟小数点再跟a，其余类推。

00分享举报