查看单个帖子
旧 2009-06-09, 09:55   #2
slgu
高级会员
 
注册日期: 2007-06-24
年龄: 70
帖子: 188
声望力: 21
slgu 正向着好的方向发展
默认 回复: 左除怎么求不出正确的线性代数方程组的解!重点部分标示

试试下列小程序:
x=[0.0032, 0.0074, 0.0107, 0.0136, 0.0166, 0.0224, 0.0282, 0.0341, 0.0399, 0.0457, 0.0545, 0.0632, 0.072, 0.0807, 0.0903, 0.1102, 0.131, 0.1557, 0.1807, 0.2107, 0.2449, 0.2852, 0.3391, 0.4204, 0.5282, 0.591, 0.671, 0.7899, 0.9693, 1.086, 1.2282, 1.4004, 1.6016, 1.8382, 2.1216, 2.466, 2.8771, 3.3882, 3.9971, 4.7438, 5.6816, 6.8193, 7.9749]';
y=[0.04, 0.07, 0.1, 0.13, 0.16, 0.22, 0.28, 0.34, 0.41, 0.47, 0.56, 0.65, 0.72, 0.82, 0.91, 1.09, 1.27, 1.46, 1.63, 1.82, 2.0, 2.17, 2.35, 2.53, 2.67, 2.73, 2.79, 2.85, 2.91, 2.94, 2.97, 3.0, 3.03, 3.06, 3.09, 3.12, 3.15, 3.18, 3.21, 3.21, 3.27, 3.3, 3.33]';
fx=@(b,x)(b(1)*x+b(2)*x.^2+b(3)*x.^3)./(1+b(4)*x+b(5)*x.^2+b(6)*x.^3);
b=[9.7865318, 32.775392, 27.227736, 1.967689, 13.118039, 7.8103828];
x1=linspace(0,8)
plot(x,y,'o','markerfacecolor','k','linewidth',3);
hold on
y1=fx(b,x1);
plot(x1,y1,'linewidth',2);
slgu 当前离线   回复时引用此帖