太原理工大学数值计算某实验报告材料

太原理工大学数值计算某实验报告材料

ID:29447239

大小:266.50 KB

页数:16页

时间:2018-12-19

太原理工大学数值计算某实验报告材料_第1页
太原理工大学数值计算某实验报告材料_第2页
太原理工大学数值计算某实验报告材料_第3页
太原理工大学数值计算某实验报告材料_第4页
太原理工大学数值计算某实验报告材料_第5页
资源描述:

《太原理工大学数值计算某实验报告材料》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用标准文案太原理工大学数值计算实验报告学院:专业班级:姓名:学号:课程名称:数值计算上机实验指导老师:精彩文档实用标准文案求非线性方程的根。一、实验目的了解数值计算方法的必要性。学会使用数值方法求解非线性方程的根。学会使用牛顿切线法求方程的解。学会使用简单迭代法求解。二、实验题目1、求方程在附近的是根,要求精度满足.(牛顿切线法)2、求方

2、程在附近的是根,求出具有思维有效数字的根近似值..(简单迭代法)三、实验内容1、程序代码clearclcphi=inline('(x^2+0.8)^(1/3)');%迭代函数x0=input('x0=');del=input('del=');N=input('N=');n=1;fprintf('%2d%f',0,x0);whilen

3、n',N);2、实验结果与分析第一题结果NewtonIterationx0=1.5del=0.001N=20kx(k)01.50000010.784472结果:0.739519即方程在x=1.5附近的根为0.739519第二题结果精彩文档实用标准文案x0=1del=0.00001N=2001.00000011.21644021.31611631.36300441.38518051.39568861.40067171.40303481.40415591.404687101.404939111.405059121.405116131.405143141.405155近似解=1.405162精

4、彩文档实用标准文案求解线性方程组(直接法或迭代法)一、实验目的1、学会求解线性方程组2、体会不同方法求解的优劣3、使用高斯-赛德尔求解下列方程组二、实验题目1、精彩文档实用标准文案三、实验代码及结果使用高斯-赛德尔迭代法求解1、程序代码clearclcn=input('n=');%矩阵的阶数A=input('A=');%系数矩阵b=input('b=');x=input('x=');%自变量epsilon=input('精度=');N=input('最大迭代次数N=');fprintf('%d:',0);fori=1:nfprintf('%f',x(i));end%以下是迭代

5、过程fork=1:N%这是第k步迭代,迭代前的向量在x0[]中,迭代后的向量在x[]中;normal=0;fori=1:nt=x(i);x(i)=b(i);forj=1:nifj~=ix(i)=x(i)-A(i,j)*x(j);endendx(i)=x(i)/A(i,i);temp=abs(x(i)-t);%求范数于迭代在同一个循环中;iftemp>normalnormal=temp;%这里用的是无穷范数endend%第i不迭代结束;fprintf('%d:',k);fori=1:nfprintf('%f',x(i));%输出迭代过程精彩文档实用标准文案endifnormal

6、ilonreturnendendfprintf('迭代%d次后仍未求得满足精度的解',N);实验结果第一题n=4A=[2,2,1,-3;-2,1,-1,-3;8,-1,3,8;10,4,4,3]b=[8,1,-1,8]x=[1,-1,2,-2]精度=0.001最大迭代次数N=100:1.000000-1.0000002.000000-2.0000001:1.000000-1.0000002.000000-2.000000>>故原方程的解为第二题n=6A=[17.031,-0.615,-2.991,1.007,-1.006,0;-1,34.211,-1,-2.1,6.3,-1.

7、7;0,0.5,13,-0.5,1,-1.5;4.501,3.11,-3.907,-61.705,12.17,8.99;0.101,-0.812,-0.017,-0.91,4.918,0.1;1,2,3,4.5,5,21.803]b=[0.23,-22.322,54,240.236,29.304,-117.818]x=[0,0,0,0,0,0]精度=0.000001最大迭代次数N=300:0.0000000.0000000.0000

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。