Matlab练习——多项式和一元方程求解

1. 一元函数 f(x) = x^3 + (x - 0.98)^2 / (x + 1.25)^3 - 5 * (x + 1 / x），求f(x) = 0 的根

%在函数文件中定义函数

function y = fun(x)

y = x ^  + (x - 0.98) ^  / (x + 1.25) ^  -  * (x +  / x);

end

%其他文件中求解

clear;

clc;

a = fzero('fun',[-,])

2. 求上述函数的极值

clear;

clc;

a = fminbnd('fun',-,0.5)

3.  创建两个多项式，进行除法运算，然后交换位置，再进行除法运算。（多项式的乘法为卷积，除法为乘法的逆运算）

clear;

clc;

p1 = [0.69552, 0.436, 0.668, 1.35];

p2 = [, 0.876, 0.645, 1.23];

[a b] = deconv(p1, p2)

[c d] = deconv(p2, p1)

4. 多项式求导

clear;

clc;

p1 = [0.69552, 0.436, 0.668, 1.35];

p2 = [0.1,0.876,0.645,1.23]

d = polyder(p1)                     %求导

[a,b] = polyder(p1,p2)              %分式求导

5.  建立一个5×5矩阵。 分别用polyval函数和polyvalm函数将矩阵带入f(x) = 1.35 + 0.668*x + 0.436*x^2 +0.69552*x^3 计算结果，进行比较。

clear;

clc;

p = [0.69552, 0.436, 0.668, 1.35];

A = rand()

a = [,0.876,0.645,1.23];

pa1 = polyval(p,A)

pa2 = polyvalm(p,A)

以上两个函数的区别：

polyval是按照数组运算规则计算多项式的值
polyvalm是按照矩阵运算规则计算多项式的值
polyval(a,X)，是计算X每个元素的多项式值
polyvalm(a,X)实际上是计算：X^2+X+E的值