
matlab 多项式带入实际值求解
现在要求一个多项式的解,比如[10*(x+a)(x+b)(x+c)]/[(x+a1)(x+b1)(x+c1)]=1,其中abca1b1c1都是常数,不知道matlab程序...
现在要求一个多项式的解,比如[10*(x+a)(x+b)(x+c)]/[(x+a1)(x+b1)(x+c1)]=1,其中a b c a1 b1 c1都是常数,不知道matlab程序是怎么写的,或者我将x代入一特定值,算前面的结果也行。算一大概区间。期待您的回答。不胜感激!
要是我要自己把x设置为一个值,看等式前面是否为1呢,或者接近于1就行。因为多项式多了,matlab也许不认。就是[10*(x+a)(x+b)(x+c)]/[(x+a1)(x+b1)(x+c1)] 将x设置为我估计的值,看结果,能否实现! 展开
要是我要自己把x设置为一个值,看等式前面是否为1呢,或者接近于1就行。因为多项式多了,matlab也许不认。就是[10*(x+a)(x+b)(x+c)]/[(x+a1)(x+b1)(x+c1)] 将x设置为我估计的值,看结果,能否实现! 展开
3个回答
展开全部
先移个项,
解10*(x+a)*(x+b)*(x+c)-(x+a1)*(x+b1)*(x+c1)=0方程,要不然运算更长时间。
syms x a b c a1 b1 c1
solve('10*(x+a)*(x+b)*(x+c)-(x+a1)*(x+b1)*(x+c1)')
结果很长,你自己运行一下吧。
顺便说句,matlab符号计算能力并不强,如果是数字就先代入数字。
由于你补充了问题,那我也补充下:
你想检验等式左边是否为1;用上面说的subs是可以的。还有个更简单办法就是:
syms x a b c a1 b1 c1 y
y=10*(x+a)*(x+b)*(x+c)/((x+a1)*(x+b1)*(x+c1));
a=1;b=1;c=1;a1=1;b1=1;c1=1; %根据需要输入数据
x=1; %对x给值
y %得出对应的值
每次给个值运行下即行,本质上同subs一样,但或许更容易理解。
解10*(x+a)*(x+b)*(x+c)-(x+a1)*(x+b1)*(x+c1)=0方程,要不然运算更长时间。
syms x a b c a1 b1 c1
solve('10*(x+a)*(x+b)*(x+c)-(x+a1)*(x+b1)*(x+c1)')
结果很长,你自己运行一下吧。
顺便说句,matlab符号计算能力并不强,如果是数字就先代入数字。
由于你补充了问题,那我也补充下:
你想检验等式左边是否为1;用上面说的subs是可以的。还有个更简单办法就是:
syms x a b c a1 b1 c1 y
y=10*(x+a)*(x+b)*(x+c)/((x+a1)*(x+b1)*(x+c1));
a=1;b=1;c=1;a1=1;b1=1;c1=1; %根据需要输入数据
x=1; %对x给值
y %得出对应的值
每次给个值运行下即行,本质上同subs一样,但或许更容易理解。
展开全部
你可以编写一个f.m文件,内容如下:
function y=f(x)
a=1;
b=1;
c=1;
a1=1;
b1=1;
c1=1;
y=(10*(x+a)*(x+b)*(x+c))/((x+a1)*(x+b1)*(x+c1));
这里的x就是你设置的值,然后在主窗口中输入f(1),就可以得到多项式的值。
f.m中a b c a1 b1 c1是我测试时用的,你可以换成自己对应的值~~
function y=f(x)
a=1;
b=1;
c=1;
a1=1;
b1=1;
c1=1;
y=(10*(x+a)*(x+b)*(x+c))/((x+a1)*(x+b1)*(x+c1));
这里的x就是你设置的值,然后在主窗口中输入f(1),就可以得到多项式的值。
f.m中a b c a1 b1 c1是我测试时用的,你可以换成自己对应的值~~
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
展开全部
先给你个例子
>> syms x y;
>> f=x^2+y^2;
>> subs(f,findsym(f),2)
ans =
y^2 + 4
更多内容见参考文章(我写的哦)
【原创】MATLAB求解符号表达式数值的方法:subs函数
>> syms x y;
>> f=x^2+y^2;
>> subs(f,findsym(f),2)
ans =
y^2 + 4
更多内容见参考文章(我写的哦)
【原创】MATLAB求解符号表达式数值的方法:subs函数
参考资料: http://hi.baidu.com/2008zhenghui/blog/item/93b05c911c789687a877a45d.html
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询