数学建模社区-数学中国

标题: 新手征程，希望大家来解答~关于多项式和常微分方程部分的问题 [打印本页]

作者: 爱的侍者 时间: 2012-8-26 13:51
标题: 新手征程，希望大家来解答~关于多项式和常微分方程部分的问题
第一个问题：p=poly(r) 和 p=poly(A) 前者由跟来计算多项式系数，后者书上写到计算A的特征多项式P，这里的A是什么？A和r不一样嘛？

第二个问题：pp=spline（x,y）;
pp
pp=
form:'pp'
breaks:[0 1 2 3 4 5 6 7 8 9 10]
coefs:[10*4 double]
pieces:10
order:4
dim:1
这里的coefs为什么是10*4，而不是10*5，order显示多项式阶数是4阶，那就是应该是需要求出5个系数啊，还有一个常数项啊？

function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=y(2)*y(3);
dy(2)=-y(1)*y(3);
dy(3)=-0.51*y(1)*y(2);

options=odeset('RelTol',le-4,'AbsTol',[le-4 le-4 le-5]);
[T,Y]=ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')

为什么要创立一个dy的空数组呢，还是3*1的，它用来干嘛的，dy=rigid(t,y)，直接创立dy=[]不可？
('RelTol',le-4,'AbsTol',[le-4 le-4 le-5]);什么意思？控制精度？我想不通
[T,y]中0 1 1]是初始条件中的常数项，[0 12]是用来干嘛的？

希望大家多多指教感激不尽。

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)