内容发布更新时间 : 2024/12/23 17:36:14星期一 下面是文章的全部内容请认真阅读。

十二、已知y的表达式y?n=100; i=1:n;

f= ；y= 十三、x是数据点的横坐标数组，li(x)是拟合数据点x?[x1,x2,?1/ii?1n2，当n?100时，求y的值，请补充完整。

,xi,x,n?1，]y?[0,0,,1,,0]的多项式，p是一个矩阵，它的第i行即为li(x)的幂系数，请补充完整。

y = zeros(1,np); y(j) = 1; p(j,:)= ;

end

function p = shape_pw(x) np = length(x); for j=1:np

十四、下面的程序是解线性方程组的一个函数文件，请补充完整。 function [x,y]=line_solution(A,b) [m,n]=size(A);y=[];

if norm (b)>0 if rank(A)==rank([A,b]) if rank(A)== x=A\\b;

else

disp('原方程组有有无穷个解，其齐次方程组的基础解系为y，特解为x');

x');

x=null(A,'r'); end ');

x=A\\b; end

else if rank(A)==n x= ;

else disp('方程组有无穷个解，基础解系为

y=null(A,'r'); x= ; end end

else

disp('方程组的最小二乘法解是：

9

return

十五、下面是求解非线性方程的牛顿法的函数文件，公式为xn?xn?1?f(xn?1)，

请补充完整。

function x = Newt_n(f_name, x0) x = x0; xb=x-999;n=0; h = 0.01; while abs(x-xb)>0.0001 n= ; xb=x; if n>300 break; end y=feval(f_name, x);

yd=(feval(f_name, x+h) - y)/ ; x = ; end fprintf('\\n Final answer = .6e\\n', x);

f?(xn?1)

10

十六、编写函数文件求半径为r的圆的周长和面积。

十七、求解下列微分方程

设著名的Lorenz方程可表示为

dx(t)???x(t)?y(t)z(t)dtdy(t) ???y(t)??z(t)dtdz(t)??x(t)y(t)??y(t)?z(t)dt 其中，??8,??10,??28。若令其初值为x(0)?y(0)?0,z(0)?1?10?10 3（1）写出用desolve 命令求解该方程的语句，并指出能否用desolve求出解析解

（2）写出用ode45求解该微分方程组的解并分别绘出状态变量的时间响应曲线和相空间三维图的程序

11

12