《数值计算方法》课程设计报告-曲线拟合-迭代法解方程组.docVIP

《数值计算方法》课程设计报告 曲线拟合 迭代法解方程组 班 级 1314132 专 业 数学与应用数学 学号姓名 指导教师 时 间 2015-01-06 前言： 本此课题我们设计了两个程序，一个是解决拟合二元函数曲线的程序，但此程序需自己先把函数化成y。=a。x。+b。的形式。另一个是雅克比迭代解方程组的程序，但此程序只能解决三元一次和四元一次方程组。这两个程序不能运算数据中含π的问题，若含π则结果不准确。本课题我们以以下两道题为例： 参数说明： [] ............................ 矩阵 (:,i) ............................ 提取矩阵中第i行的数据 ./ ............................ 点除 ~= ............................ 不等于 Sum(sum（）)......................... 矩阵中所有元素的和 .* ............................ 点乘 \ ............................ 右除 pi ............................ π for end ............................ 循环语句 while end............................ 循环语句 if else end.......................... 判断语句 abs（） ............................ 调用取绝对值函数 ............................ 与 ‘’ ............................ 字符 exp（n） ............................ e的n次方 matlab程序流程图： 6. D=[1 60;2 30;3 20;4 15]; %x y的值 x0=D(:,1); %调用矩阵第一列的数据 y0=D(:,2); %调用矩阵第二列的数据 y=log(y0); %对矩阵中的元素进行运算 x=x0; A=x~=pi; %取不等于π的元素并表示为1 i0=sum(sum(A)); %求矩阵中元素的和 i1=sum(sum(x)); B=x.*x; %矩阵中每个元素平方后生成新矩阵 i2=sum(sum(B)); C=x.*y; %x矩阵与y矩阵对应位置元素相乘 e0=sum(sum(y)); e1=sum(sum(C)); E=[i0 i1;i1 i2]; F=[e0;e1]; H=E\F; %矩阵左除运算 a0=H(1); %取单行矩阵第一个元素 b0=H(2); %取单行矩阵第二个元素 a=exp(a0) %调用函数e的a0次方 b=b0 22. function f=f3(A); %新建一个函数 Format long; %增加小数点后位数 B=A(:,1); %调用矩阵第一列数据 C=B~=pi; %取不等于π的元素并表示为1 e=sum(sum(C)); %矩阵元素和 x=0;y=0; for(i=1:e) %循环语句 x(i)=0;y(i)=1; %给x，y赋初始值 end; if e==3; %判断语句，下为循环语句 while(abs((y(1)-x(1))/y(1))10^(-3)abs((y(2)-x(2))/y(2))10^(-3)