clear;clc;
syms x a;
m=5；%自己改
y=(11/6-3*x 3/2*x^2-1/3*x^3)^a
f=taylor(y,m 1,x);
w=sym(zeros(m 1,1));
w(1)=subs(f,x,0);
f=f-w(1);
for n=m:-1:2
w(n 1)=subs(f-subs(f,x^n,0),x^n,1);
f=f-w(n 1)*x^n;
end
w(2)=subs(f,x,1)
注意，因為matlab數組下標從1開始，因此我這裡w(1)是常數項，w(2)是一次項，以此類推，即
y=w(1) w(2)*x w(3)*x^2 .... w(m 1)*x^m
【1】把函數變形
>>f=sym('2*x^3 3*x^2 21*x 4-(3*a*x^3 b*x^2 c*x d)=0')
f =
2*x^3 3*x^2 21*x 4-(3*a*x^3 b*x^2 c*x d)=0
【2】用collect合併同類項
>>ff=collect(f):
#(2-3*a)*x^3 (3-b)*x^2 (21-c)*x 4-d = 0
【3】用maple擷取多項式係數，如果很多，可以用迴圈語句。
>>c3=maple('coeff',ff,x,3)
c3 =2-3*a
>>c1=maple('coeff',ff,x,1)
c1 =21-c
>>c2=maple('coeff',ff,x,2)
c2 =3-b
>>c0=maple('coeff',ff,x,0)
c0 =4-d
補充：
這次變成這樣，程序倒是通了，我不是很滿意，咱們共同把它整理好，如何？
syms a b c d x
%【1】把函數變形
f=sym('2*x^3 3*x^2 21*x 4-(3*a*x^3 b*x^2 c*x d)')
N=3;
for i=0:N
temp=maple('coeff',f,x,N-i);
cp(1,i 1)={temp};
end
celldisp(cp);
再補充：這次總算是解決了，就是看起來很笨，不是很理想，湊合著用，當然我相信可以修改的很漂亮。
syms a b c d x
f=sym('2*x^3 3*x^2 21*x 4-(3*a*x^3 b*x^2 c*x d)')
N=3;
for i=0:N
temp=maple('coeff',f,x,N-i);
temp1(i 1)=temp;
#end
cp=temp1
a=solve(cp(1)), b=solve(cp(2)), c=solve(cp(3)), d=solve(cp(4))
運行結果：
a =2/3
b =3
c =21
d =4
首先，多項式是動態的，所以這必須是matlab的輸入項；
其次，多項式的matlab表達要清楚，是將多項式降冪排列後提取其係數來表示該多項式的-n次多項式用n 1維向量表示；例如多項式3*x^2 5 在matlab中的表示為[3 0 5];
最後，多項式函數值的matlab法要明白，就是指令polyval。
綜合上述，M文件如下：
function val = fpolyval(p,x)
% 函數 fpolyval 函數：多項式 p 在 x 處的函數值 val.
% 輸入項 p 是多項式依降冪排列後的係數；
val = polyval(p,x);
例如：3*x^2 5在x=1,2處的值
>>p=[3 0 5];
>>x=[1 2];
>>val=fpolyval(p,x)
val =
8 17
以上是使用MATLAB計算多項式的泰勒級數展開係數的詳細內容。更多資訊請關注PHP中文網其他相關文章！
