怎样含有字符的sym格式的转化数组乘积

lihaitao123

首先是我编的切比雪夫的公式：

1. function T=Chebyshev2(n)
   
2.    syms x
   
3.    T(1:n)=sym(zeros(1,n));
   
4.    T(1)=1;
   
5.    T(2)=x;
   
6.    for i=3:n
   
7.        T(i)=2*x*T(i-1)-T(i-2);
   
8.    end

复制代码

然后是用最小二乘法拟合函数

1. function f=li2_23(c,x)
   
2.     T=Chebyshev2(4);
   
3.      f=c(1)*T(1)+c(2)*T(2)+c(3)*T(3)+c(4)*T(4);

复制代码

最后是主程序：

1. clear all
   
2. f=@(x)sin(x)
   
3. x=-6:0.2:6;
   
4. y=f(x);
   
5. plot(x,y,'*')
   
6. grid on;
   
7. hold on;
   
8. c0=ones(1,4);  %迭代初值
   
9. c=lsqcurvefit(@li2_23,c0,x,y);  %调用函数
   
10. g=@(x)c(1)*T(1)+c(2)*T(2)+c(3)*T(3)+c(4)*T(4)  %构建拟合函数
    
11. x_g=-6:0.1:6;
    
12. y_g=g(x_g);
    
13. plot(x_g,y_g)
    
14. hold off
    
15. legend('原始数据','拟合数据');
    
16. c

复制代码

??? Error using ==> lsqncommon at 191
Function value and YDATA sizes are incommensurate
但是把T=Chebyshev2(4)单独拿出运行，并且在多项式乘方加点。T=[ 1, x, 2*x.^2 - 1, 2.*x.*(2.*x.^2 - 1) - x];然后把T(1),T(2),T(3),T(4)代入就可以运行。所以我的问题是怎么Chebyshev2(n)生成 的sym矩阵方便的转化为T=[ 1, x, 2*x.^2 - 1, 2.*x.*(2.*x.^2 - 1) - x]；（乘积带点）。毕竟当阶数较高时候，挨个展开带进去不太方便。谢谢