计算方法实验报告一谜底 许晶
计算方法实验报告
计算方法实验报告计算方法实验报告概述：计算方法是一门研究如何用计算机解决数学问题的学科。
在本次实验中，我们将学习和应用几种常见的计算方法，包括数值逼近、插值、数值积分和常微分方程求解。
通过实验，我们将深入了解这些方法的原理、应用场景以及其在计算机科学和工程领域的重要性。
数值逼近：数值逼近是一种通过使用近似值来计算复杂函数的方法。
在实验中，我们通过使用泰勒级数展开和牛顿迭代法等数值逼近技术，来计算函数的近似值。
这些方法在科学计算和工程领域中广泛应用，例如在信号处理、图像处理和优化问题中。
插值：插值是一种通过已知数据点来估算未知数据点的方法。
在实验中，我们将学习和应用拉格朗日插值和牛顿插值等方法，以及使用这些方法来构造函数的近似曲线。
插值技术在数据分析、图像处理和计算机图形学等领域中具有重要的应用价值。
数值积分：数值积分是一种通过将函数曲线划分为小矩形或梯形来估算函数的积分值的方法。
在实验中，我们将学习和应用矩形法和梯形法等数值积分技术，以及使用这些方法来计算函数的近似积分值。
数值积分在物理学、金融学和统计学等领域中被广泛使用。
常微分方程求解：常微分方程求解是一种通过数值方法来求解微分方程的方法。
在实验中，我们将学习和应用欧拉法和龙格-库塔法等常微分方程求解技术，以及使用这些方法来求解一些常见的微分方程。
常微分方程求解在物理学、生物学和工程学等领域中具有广泛的应用。
实验结果：通过实验，我们成功地应用了数值逼近、插值、数值积分和常微分方程求解等计算方法。
我们得到了准确的结果，并且在不同的应用场景中验证了这些方法的有效性和可靠性。
这些实验结果将对我们进一步理解和应用计算方法提供重要的指导和支持。
结论：计算方法是计算机科学和工程领域中的重要学科，它提供了解决复杂数学问题的有效工具和方法。
通过本次实验，我们深入了解了数值逼近、插值、数值积分和常微分方程求解等计算方法的原理和应用。
这些方法在科学研究、工程设计和数据分析等领域中具有广泛的应用价值。
(完整word版)计算方法A上机实验报告
计算方法A上机实验报告姓名：苏福班级：硕4020 学号：3114161019一、上机练习目的1）复习和巩固数值计算方法的基本数学模型，全面掌握运用计算机进行数值计算的具体过程及相关问题。
2）利用计算机语言独立编写、调试数值计算方法程序，培养学生利用计算机和所学理论知识分析解决实际问题的能力。
二、上机练习任务1）利用计算机语言编写并调试一系列数值方法计算通用程序，并能正确计算给定题目，掌握调试技能。
2）掌握文件使用编程技能，如文件的各类操作，数据格式设计、通用程序运行过程中文件输入输出运行方式设计等。
3）写出上机练习报告。
三、上机题目1. 共轭梯度法求解线性方程组。
（第三章）2. 三次样条插值（第四章）3. 龙贝格积分（第六章）4. 四阶龙格－库塔法求解常微分方程的初值问题四、上机报告题目1：共轭梯度法求解线性方程组1．算法原理共轭梯度法是把求解线性方程组的问题转化为求解一个与之等价的二次函数极小值的问题。
从任意给定的初始点出发，沿一组关于矩阵A共轭的方向进行线性搜索，在无舍入误差的假定下，最多迭代n 次（其中n 为矩阵A 的阶数），就可求得二次函数的极小值，也就求得了线性方程组Ax b =的解。
定理：设A 是n 阶对称正定矩阵，则x *是方程组Ax b =的解得充分必要条件是x *是二次函数1()2TT f x x Ax b x =-的极小点，即 ()()min nx R Ax b f x f x **∈=⇔=共轭梯度法的计算公式：(0)(0)(0)()()()()(1)()()(1)(1)(1)()()()(1)(1)()k T k k k T k k k k k k k k T k k k T k k k k k d r b Ax r d d Ad xx d r b Ax r Ad d Ad d r d ααββ++++++⎧==-⎪⎪=⎪⎪=+⎪⎨=-⎪⎪⎪=-⎪⎪=+⎩2. 程序框图（1）编写共轭梯度法求解对称正定矩阵的线性方程组见附录（myge.m）：function x=myge(A,b)输入对称正定矩阵及对应的列向量，初始向量设为0，精度取为810 。
计算方法实验报告.doc
计算方法实验报告实验报告一、求方程fxx3-sinx-12x1的全部根, ε1e-6 1、用一般迭代法; 2、用牛顿迭代法; 并比较两种迭代的收敛速度。
一、首先，由题可求得. 其次，分析得到其根所在的区间。
①令，可得到. ②用一阶导数分析得到和两个函数的增减区间；再用二阶导数分析得到两个函数的拐点以及凹凸区间. ③在直角坐标轴上描摹出和的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。
经过估计，得到根所在的区间为，和. 1、一般迭代法（1）算法步骤设为给定的允许精度，迭代法的计算步骤为①选定初值.由确定函数，得等价形式. ②计算.由迭代公式得.③如果，则迭代结束，取为解的近似值；否则，用代替，重复步骤②和步骤③. （2）程序代码①在区间内，代码clc x0-3.5; 初值iter_max100; 迭代的最大次数ep1e-6; 允许精度k0; while kiter_max k从0开始到iter_max循环x1sinx012*x0-1.1/3; 代入，算出的值if absx1-x0ep 与允许精度作比较break; 条件成立，跳出循环end x0x1; 条件不成立，用代替kk1; k加1 end x_starx1, iterk 为解的近似值，iter为迭代次数运行结果x_star -3.4101 ；iter 14 ②在区间内，代码clc x00.5; 初值iter_max100; 迭代的最大次数ep1e-6; 允许精度k0; while kiter_max k从0开始到iter_max循环x11/12*x0.3-sinx01; 代入，算出的值if absx1-x0ep 与允许精度作比较break; 条件成立，跳出循环end x0x1; 条件不成立，用代替kk1; k加1 end x_starx1, iterk 为解的近似值，iter 为迭代次数结果x_star 0.07696；iter 6 ③在区间内，代码clc x03.5; 初值iter_max100; 迭代的最大次数ep1e-6; 允许精度k0; while kiter_max k从0开始到iter_max循环x1sinx012*x0-1.1/3; 代入，算出的值if absx1-x0ep 与允许精度作比较break; 条件成立，跳出循环end x0x1; 条件不成立，用代替kk1; k加1 end x_starx1, iterk 为解的近似值，iter 为迭代次数运行结果x_star 3.4101 ；iter 10 2、牛顿迭代法（1）算法步骤①选定初值，计算，. ②按公式迭代，得新的近似值，并计算，. ③对于给定的允许精度，如果，则终止迭代，取；否则，，在转回步骤②计算. （2）程序代码①在区间内，clc x1-3.5; 初值k0; while k100 k从0开始到100循环x0x1; 将初值赋给f0x0.3-sinx0-12*x01; 计算f13*x0.2-cosx0-12; 计算x1x0-f0/f1; 计算得到新的近似值if absx1-x0 1.0e-6 与允许精度作比较break; 条件成立，跳出循环end kk1; 条件不成立，k加1 end x_starx1, iterk 为解的近似值，iter为迭代次数运行结果x_star -3.4911；iter 2 ②在区间内，clc x10.5; 初值k0; while k100 k从0开始到100循环x0x1; 将初值赋给f0x0.3-sinx0-12*x01; 计算f13*x0.2-cosx0-12; 计算x1x0-f0/f1; 计算得到新的近似值if absx1-x0 1.0e-6 与允许精度作比较break; 条件成立，跳出循环end kk1; 条件不成立，k加1 end x_starx1, iterk 为解的近似值，iter为迭代次数运行结果x_star 0.07696 ；iter 3 ③在区间内，clc x13.5; 初值k0; while k100 k从0开始到100循环x0x1; 将初值赋给f0x0.3-sinx0-12*x01; 计算f13*x0.2-cosx0-12; 计算x1x0-f0/f1; 计算得到新的近似值if absx1-x0 1.0e-6 与允许精度作比较break; 条件成立，跳出循环end kk1; 条件不成立，k加1 end x_starx1, iterk 为解的近似值，iter为迭代次数运行结果x_star 3.4101；iter 3 3、运行结果x_star iter x_star iter x_star iter 一般迭代3.4101 14 0.7696 6 3.4101 10 牛顿法3.4911 2 0.70696 3 3.4101 3 4、结果分析从这题的结果可以看出，牛顿迭代法的迭代速度比一般迭代法的速度要快，牛顿法的迭代次数比较少。
大学计算方法实验报告
《计算方法》实验报告实验题目实验报告1：非线性方程组的求解···················P1~2实验报告2：线性方程组解法·······················P3~4 实验报告3：Lagrange 插值多项式··················P5~7姓名：学号：班级：指导老师：时间：专业 序号 日期实验报告1：非线性方程组的求解【实验目的】1．用MATLAB 来实践进行牛顿法的变形，即对牛顿法进行了修正，使其应用更为方便，掌握用MATLAB 运用割线法求解非线性方程组。
2．运用MATLAB 进行隐函数作图。
【实验内容】[方法] 设a,b 为迭代初值,求两点(a,f(a)) 与 (b,f(b)) 的连线(割线)与 x 轴的交点记为 c ，再把迭代初值换成 b,c,重复计算.[要求] 把下面程序复制为新的 M-文件,去掉开头的 %再把 '?' 部分改写正确就是一个完整的程序,找前面一个例子试算【解】在牛顿迭代公式中用差商代替导数。
带入初值（a,f(a)）,（b,f(b)）,两点的连线与x 轴的交点作为c ，再把迭代初值换为b ，c ，重复计算。
【计算机求解】以y= x-exp(-x)为例初值a=0，b=1，误差不超过1.0*10^(-5)进行计算。
计算方法实验报告1
目录实验一牛顿下山法求解非线性方程的根 (3)一、实验目的 (3)二、实验内容 (3)三、基本原理 (3)四、算法设计与实现 (3)五、输入与输出 (4)六、源代码 (4)实验二高斯——赛德尔法求线性方程组 (6)一、实验目的 (6)二、实验内容 (6)三、算法基本原理 (6)四、算法设计与实现 (6)五、计算用例的参考输出 (8)六、源代码 (8)实验三高斯消去法 (10)一、实验目的 (10)二、实验内容 (10)三、算法基本原理 (10)四、算法设计与实现 (10)五、计算用例的参考输出 (12)六、源代码 (12)实验四ROMBERG算法 (14)一、实验目的 (14)二、实验内容 (14)三、算法基本原理 (14)四、算法设计与实现 (14)五、计算用例的参考输出 (16)六、源代码 (16)实验五RUNGEKUTTA算法 (18)一、实验目的 (18)二、实验内容 (18)三、算法基本原理 (18)四、算法设计与实现 (19)五、计算用例的参考输出 (19)六、源代码 (20)感悟与体会 (21)实验一 牛顿下山法求解非线性方程的根一、实验目的(1) 熟悉非线性方程求根简单迭代法，牛顿迭代及牛顿下山法 (2) 能编程实现简单迭代法，牛顿迭代及牛顿下山法 (3) 认识选择迭代格式的重要性(4)对迭代速度建立感性的认识；分析实验结果体会初值对迭代的影响二、实验内容用牛顿下山法解方程 013=--x x （初值为0.6）输入：初值，误差限，迭代最大次数，下山最大次数 输出：近似根各步下山因子三、基本原理牛顿下山公式：)()(1k k k k x f x f x x '-=+λ 下山因子 ，，，，322121211=λ 下山条件|)(||)(|1k k x f x f <+四、算法设计与实现流程图关键点：while(fabs(f(x1))>fabs(f(x0))) { cout<<"λ="<<r<<" "<<endl; x1=x0-r*f(x0)/df(x0); r=r/2; k++;if(k>n){cout<<"迭代失败"<<endl;break;} }//求解下山因子 while (fabs(x0-x1)>=eps) { if(df(x0)==0){cout<<"***无法迭代***"<<endl;break;} x0=x1; x1=x0-r*f(x0)/df(x0); c++; if(c>m){cout<<"迭代失败"<<endl;break;}}//迭代过程五、输入与输出x0=0.6;e=0.0000001;m=100;n=100输出六、源代码#include<iostream> #include<math.h> using namespace std;图3.2牛顿下山算法流程⇐⇐⇐⇐double f(double x){return (x*x-1)*x-1;}double df(double x){return 3*x*x-1;}double newton(double x0,double eps,int n,int m){if(df(x0)==0)cout<<"***无法迭代***"<<endl;else{int k,c;//下山次数,迭代次数double x1=x0-f(x0)/df(x0),r=1;while(fabs(f(x1))>fabs(f(x0))){cout<<"λ="<<r<<" "<<endl;x1=x0-r*f(x0)/df(x0);r=r/2;k++;if(k>n){cout<<"迭代失败"<<endl;break;}}//求解下山因子while (fabs(x0-x1)>=eps){if(df(x0)==0){cout<<"***无法迭代***"<<endl;break;}x0=x1;x1=x0-r*f(x0)/df(x0);c++;if(c>m){cout<<"迭代失败"<<endl;break;}}//迭代过程return x1;}}void main(){double x0,x1,e;x0=0.6;e=0.0000001;x1=newton(x0,e,100,100);cout<<"x="<<x1<<endl;}实验二 高斯——赛德尔法求线性方程组一、实验目的(1) 熟悉求解线性方程组的有关理论和方法； (2) 能编程实现雅可比及高斯-塞德尔迭代法；(4)根据不同类型的方程组，选择合适的数值方法。
计算方法与实习上机实验报告
计算方法与实习上机实验报告一、引言本文旨在介绍和展示我们在“计算方法”课程中的实习上机实验环节所完成的一些关键任务和所取得的成果。
该实验课程的目标是让我们更深入地理解和应用各种计算方法，并在实际操作中提高我们的编程和问题解决能力。
二、实验内容与目标实验的主要内容是利用各种计算方法解决实际数学问题。
我们被要求使用编程语言（如Python或Java）来实现和解决这些问题。
这些问题包括使用牛顿法求解平方根，使用蒙特卡洛方法计算圆周率，以及使用最优化方法求解函数的最小值等。
实验的目标不仅是让我们掌握计算方法的基本理论，更是要让我们能够在实际操作中运用这些方法。
我们需要在实习过程中，通过与同伴们合作，共同解决问题，提高我们的团队合作能力和问题解决能力。
三、实验过程与问题解决策略在实验过程中，我们遇到了许多问题，如编程错误、理解困难和时间压力等。
我们通过相互讨论、查阅资料和寻求教师帮助等方式，成功地解决了这些问题。
例如，在实现牛顿法求解平方根时，我们一开始对导数的计算和理解出现了一些错误。
但我们通过查阅相关资料和讨论，最终理解了导数的正确计算方法，并成功地实现了牛顿法。
四、实验结果与结论通过这次实习上机实验，我们不仅深入理解了计算方法的基本理论，还在实际操作中提高了我们的编程和问题解决能力。
我们的成果包括编写出了能有效求解平方根、计算圆周率和求解函数最小值的程序。
这次实习上机实验非常成功。
我们的团队不仅在理论学习和实践操作上取得了显著的进步，还在团队合作和问题解决方面积累了宝贵的经验。
这次实验使我们对计算方法有了更深的理解和认识，也提高了我们的编程技能和解决问题的能力。
五、反思与展望回顾这次实验，我们意识到在实验过程中，我们需要更好地管理我们的时间和压力。
在解决问题时，我们需要更有效地利用我们的知识和资源。
在未来，我们希望能够更加熟练地运用计算方法，并能够更有效地解决问题。
我们也希望能够将所学的计算方法应用到更广泛的领域中，如数据分析、科学研究和工业生产等。
计算课实验报告总结(3篇)
第1篇一、实验背景随着信息技术的飞速发展，计算课已成为现代教育中不可或缺的一部分。
通过计算课的学习，学生可以掌握计算机基本操作、编程语言以及算法设计等知识，为今后从事相关工作奠定基础。
本次实验旨在通过实际操作，加深对所学知识的理解，提高动手能力和团队协作能力。
二、实验目的1. 熟悉计算机基本操作，掌握常用软件的使用方法；2. 学习一种编程语言，理解编程思想，实现基本算法；3. 培养团队协作精神，提高动手实践能力；4. 提高对计算课重要性的认识，激发学习兴趣。
三、实验内容本次实验主要包括以下内容：1. 计算机基本操作：熟练使用计算机操作系统，掌握文件管理、系统设置等基本操作；2. 编程语言学习：选择一种编程语言（如Python、Java等），学习基本语法、数据结构、算法等知识；3. 算法实现：设计并实现一个简单算法，如排序、查找等；4. 项目实践：分组完成一个小型项目，如制作一个简单的网页、编写一个计算器程序等。
四、实验过程1. 实验准备：了解实验内容，预习相关理论知识，准备好实验所需的计算机和软件；2. 实验操作：按照实验指导书进行操作，记录实验步骤和结果；3. 团队协作：分组讨论，分工合作，共同完成实验任务；4. 结果分析：对实验结果进行分析，总结经验教训。
五、实验结果与分析1. 计算机基本操作：通过实验，掌握了计算机基本操作，如文件管理、系统设置等，提高了计算机应用能力；2. 编程语言学习：学习了所选编程语言的基本语法、数据结构、算法等知识，为今后深入学习打下了基础；3. 算法实现：实现了排序、查找等基本算法，加深了对算法原理的理解；4. 项目实践：分组完成了一个小型项目，如制作了一个简单的网页、编写了一个计算器程序等，提高了团队协作能力和动手实践能力。
六、实验总结1. 计算课实验对提高学生计算机应用能力具有重要意义，有助于培养学生动手实践能力和团队协作精神；2. 实验过程中，要注重理论与实践相结合，不断总结经验教训，提高实验效果；3. 在今后的学习中，要继续努力，深入学习计算课相关知识，为将来从事相关工作打下坚实基础。
计算方法实验报告(一)答案 许晶
计算方法实验报告年级班级：学号：姓名：日期：计算方法实验（一）实验目的：1、掌握MATLAB控制语句2、熟悉数组运算3、MATLAB图形处理功能4、MATLAB程序初步设计实验内容：一、数组运算1、已知向量a=1:12, 用reshape命令将其编程2*6矩阵。
>> a=1:12;>> reshape(a,2,6)ans =1 3 5 7 9 112 4 6 8 10 122、已知a=[1 2 3;2 3 4;3 4 5]; b=[1 1 1;2 2 2;3 3 3]; 在计算机上运行下面的命令，并写出计算结果。
a.*b a*b a.\b a./b3、已知b=[1 1 1;2 2 2;3 3 3]; 在计算机上运行 b^3 b.^3 b*b*b并写出计算结果。
二、程序设计（写出程序，并给出计算结果）1、利用循环语句编写程序，求1^2+2^2+3^2+…+100^2。
2、利用选择语句编写函数文件demo3实现sgn函数功能，并计算demo3(0) demo3(90) demo3(-12)function f=demo3(x)if x>0f=1;elseif x==0f=0;else x<0f=-1;enddemo3(0)ans =demo3(90)ans =1demo3(90)ans =-1三、MATLAB图形处理（写出程序，给出图形）1、在0≤x≤2π区间内，绘制曲线 y=2e^(0.5x)cos(4πx)2、在[0, 20π]之间，绘制三维曲线 x=sin(t); y=cos(t); z=tsin(t)cos(t);图形要求：标题为“'Line in 3-D Space”，x轴，y轴, z轴分别添加标注“X”“Y”“Z”。
计算方法实验指导与实验报告
6
实验题目 2 龙贝格(Romberg)积分法
方法概要：利用复化梯形求积公式、复化辛普生求积公式、复化柯特斯求积公式 的误差估计式计算积分 计算公式：
b
a
f ( x)dx 。记 h
ba ， xk a k h ， k 0,1, n
, n ，其
Tn
1 n h k 1[ f ( xk 1 ) f ( xk )] 2
2
实验题目 1 拉格朗日(Lagrange)插值
方法概要： 给定平面上 n 1 个不同的数据点 ( xk , f ( xk )) ， k 0,1, 则满足条件
, n ， xi x j ， i j ；
Pn ( xk ) f ( xk ) ， k 0,1,
的 n 次拉格朗日插值多项式
, n ，构造 Pn ( x ) ，利用拉格朗日插值多项式 Pn ( x ) 作
为 f ( x) 的近似值。分别取 n 5 ， n 10 ， n 20 ，同时计算 Pn ( x ) 在 x 0.95 ，
x 0.05 ， x 0.05 ， x 0.95 处的函数值。
（2）设 f ( x) e ， x [1,1] ，考虑非等距节点的拉格朗日插值多项式 Pn ( x ) ，
问题 2 插值区间越小越好吗？ 考虑下面两个拉格朗日插值问题：
1 ， x [1,1] ，考虑等距节点的拉格朗日插值多项式 Pn ( x ) ， 1 x2 2.0 xk 1.0 k h ， 即将区间 [1,1] 进行 n 等分， 记h ， 构造 Pn ( x ) ， k 0,1, , n ， n
问题 4 考虑拉格朗日插值问题，内插比外推更可靠吗？ 考虑下面两个拉格朗日插值问题： （1）设 f ( x)
计算方法上机实验报告
《计算方法》上机实验报告班级:XXXXXX小组成员：XXXXXXXXXXXXXXXXXXXXXXXXXXXX任课教师：XXX二〇一八年五月二十五日前言通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写.以下为本次上机实验报告，按照实验内容共分为六部分.实验一：一、实验名称及题目： Newton 迭代法例2。
7（P38):应用Newton 迭代法求在附近的数值解，并使其满足.二、解题思路：设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交点的横坐标)(')(0001x f x f x x -=,称1x 为'x 的一次近似值，过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标)(')(1112x f x f x x -=称2x 为'x 的二次近似值，重复以上过程，得'x 的近似值序列{}n x ,把)(')(1n n n n x f x f x x -=+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。
三、Matlab 程序代码：function newton_iteration （x0,tol) syms z ％定义自变量 format long ％定义精度 f=z*z*z —z-1;f1=diff(f ）；%求导 y=subs(f ，z ，x0);y1=subs （f1,z,x0)；%向函数中代值 x1=x0-y/y1; k=1;while abs(x1-x0）>=tol x0=x1；y=subs(f,z,x0)； y1=subs （f1,z,x0）; x1=x0—y/y1;k=k+1; endx=double(x1) K四、运行结果：实验二：一、实验名称及题目：Jacobi 迭代法例3。
计算方法实验报告
计算方法实验报告实验目的：1.掌握计算方法的基本概念和算法；2.熟悉常见计算方法的实现步骤和注意事项；3.学会使用计算方法解决实际问题。
实验内容：1.实现二分法；2.实现牛顿迭代法；3.实现高斯消去法。
实验步骤：1.实现二分法：1.1定义函数f(x)；1.2 确定初始区间[a, b]和精度tol；1.3计算区间中点c；1.4判断f(a)和f(c)的符号关系并更新区间；1.5重复步骤1.3和1.4直到满足精度要求。
2.实现牛顿迭代法：2.1定义函数f(x)；2.2定义函数f的导数；2.3给定初始点x0；2.4计算f(x0)和f'(x0)；2.5计算下一个点的近似值x1=x0-f(x0)/f'(x0)；2.6重复步骤2.4和2.5直到满足收敛条件。
3.实现高斯消去法：3.1输入线性方程组的系数矩阵A和右端向量b；3.2构造增广矩阵[A，b]；3.3进行主元素消去，得到梯形矩阵U和新的右端向量b；3.4回代求解，得到解向量x。
实验结果分析：1.二分法的主要优点是收敛稳定，但需要事先给定初始区间；2.牛顿迭代法的主要优点是收敛速度快，但需要事先给定初始点和收敛条件；3.高斯消去法的主要优点是适用于任何线性方程组，但需要事先进行主元素消去和回代的操作。
实验总结：通过本次实验，我深入理解了计算方法的基本概念和算法，并掌握了二分法、牛顿迭代法和高斯消去法的实现步骤和注意事项。
这些方法在解决实际问题中具有重要的应用价值。
实验过程中，我也遇到了一些困难和挑战，例如初始值的选择和收敛条件的判断。
通过不断的调试和优化，最终成功解决了这些问题。
本次实验不仅提高了我的编程能力，也增加了我的数学建模能力。
希望今后能够继续深入学习计算方法，并将其应用于更加复杂的实际问题中。
计算方法与实习的实验报告范文
计算方法与实习的实验报告范文1舍入误差与数值稳定性1.1目的与要求（1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；（2）通过上机计算，了解舍入误差所引起的数值不稳定性。
1.2舍入误差和数值稳定性1.2.1概要舍入误差在计算方法中是一个很重要的概念。
在实际计算中如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。
因此，选取稳定的算法在实际计算中是十分重要的。
1.2.2程序和实例对n=0,1,2，…，40计算定积分1某n某5n0d某。
算法利用递推公式yn=ln6-ln50.182322。
程序如下：#include#includevoidmain(){doubley_0=log(6.0/5.0),y_1;intn=1;printf(\while(1){y_1=1.0/n-5某y_0;printf(\if(n>=40)break;y_0=y_1;n++;if(n%2==0)printf(\}}1n5y（n=1,2,…，40）取y0=11某50d某=2方程求根2.1实验目的（1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点；（2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法作比较。
2.2二分法2.2.1算法给定区间[a,b],并设f(a)与f(b)符号相反，取ε为根的容许误差，δ为|f(某)|的容许误差。
1令c=(a+b)/2；2如果(c-a)0,则令a=c;否则令b=c,重复1,2,3。
2.2.2程序与实例求方程f(某)=某4某100在1.5附近的根。
程序如下：#include#include#defineep5e-6#definedelta1e-6floatBiection(floata,floatb,float(某f)(float)){floatc,fc,fa=(某f)(a),fb=(某f)(b);intn=1;printf(\二分次数\\t\\tc\\t\\tf(c)\\n\while(1){if(fa某fb>0){printf(\不能用二分法求解\c=(a+b)/2,fc=(某f)(c);printf(\if(fab(fc)returnc;}floatf(float某){return某某某某某+4某某某某-10;32}intmain(){floata=1,b=2;float某;某=Biection(a,b,f);printf(\方程的根为%f\return0;}3线性方程组数值解法3.1目的与要求（1）熟悉求解线性方程组的有关理论和方法；（2）会编制列主元消去法、LU分解法、雅克比及高斯赛德尔迭代法的程序；（3）通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。
计算方法实验报告
1. 熟悉并掌握常用的计算方法，包括数值积分、数值微分、线性方程组求解等。
2. 培养运用计算机进行数值计算的能力。
3. 增强对数值计算误差的分析和判断能力。
二、实验环境1. 操作系统：Windows 102. 编程语言：Python3.83. 库：NumPy、SciPy、Matplotlib三、实验内容1. 数值积分（1）函数：f(x) = x^2（2）区间：[0, 1]（3）方法：梯形法、辛普森法、复合梯形法2. 数值微分（1）函数：f(x) = e^x（2）点：x = 1（3）方法：有限差分法、中点法、牛顿法3. 线性方程组求解（1）方程组：2x + 3y - z = 8-x + 2y + 2z = -3x - y + 3z = 5（2）方法：高斯消元法、LU分解法1. 数值积分（1）编写函数f(x) = x^2（2）定义积分区间[0, 1]（3）实现梯形法、辛普森法、复合梯形法（4）计算积分结果2. 数值微分（1）编写函数f(x) = e^x（2）定义点x = 1（3）实现有限差分法、中点法、牛顿法（4）计算导数结果3. 线性方程组求解（1）定义方程组系数矩阵A和常数向量b（2）实现高斯消元法、LU分解法（3）求解方程组（4）输出解向量x五、实验结果与分析1. 数值积分（1）梯形法：积分结果约为1.6667（2）辛普森法：积分结果约为1.6447（3）复合梯形法：积分结果约为1.6458分析：三种方法计算结果接近，但辛普森法误差最小。
2. 数值微分（1）有限差分法：导数结果约为2.7183（2）中点法：导数结果约为2.7183（3）牛顿法：导数结果约为2.7183分析：三种方法计算结果一致，误差较小。
3. 线性方程组求解（1）高斯消元法：解向量x = [2, 1, 1]（2）LU分解法：解向量x = [2, 1, 1]分析：两种方法求解结果一致，且解向量正确。
六、实验总结本次实验通过Python编程，实现了数值积分、数值微分和线性方程组求解。
计算方法实习报告
一．给出一个有效的算法和无效的算法计算积分y(n)=∫(x^n)/(4x+1)dx,n=0,1,2,…,10,积分限为（0，1）1．有效算法利用递推公式y(n)=-y(n-1)/4+1/(4n)，取y0=(log5)/4程序为：#include<iostream.h>#include<math.h>void main(){ double y0,y1;y0=1/4.0*log(5.0);cout<<"y0="<<y0<<" ";for(int n=1;n<=10;n++){ y1=-1.0/4.0*y0+1.0/(4.0*n);cout<<"y"<<n<<"="<<y1<<" ";y0=y1;if(n%3==0) cout<<endl;}}其结果为： y0=0.402359 y1=0.14941 y2=0.0876475 y3=0.0614215y4=0.0471446 y5=0.0382138 y6=0.0321132y7=0.027686 y8=0.0243285 y9=0.0216957y10=0.0195761 Press any key to continue2．无效算法利用递推公式y(n-1)=-4y(n)+1/n, 又由广义积分中值定理可得y(n)=1/((4n+1)(4ζ+1)), ζ∈(0,1),则1/(5(n+1))<y(n)<1/(n+1), 所以可取y(n)≈[1/(5(n+1))+ 1/(n+1)]/2=3/(5(n+1)程序为：#include<iostream.h>#include<math.h>void main(){ float y9,y10;y10=3.0 /55.0;cout<<"y10="<<y10<<" ";for(int n=9;n>=0;n--){ y9=-4.0*y10+1.0/(n+1);cout<<"y"<<n<<"="<<y9<<" ";y10=y9;if(n%3==0) cout<<endl;}}其结果为：y10=0.0545455 y9=-0.118182y8=0.583838 y7=-2.21035 y6=8.98427y5=-35.7704 y4=143.282 y3=-572.877y2=2291.84 y1=-9166.86 y0=36668.43．心得体会由有效算法与无效算法的结果可以知道：有效算法的误差的传递是逐步缩小的，而无效算法的误差的传递是逐步扩大的。
计算方法上机实习报告5 - 华中科技大学
计算方法上机实习报告 5一． 提出问题（1） 给出概率积分 的数据表：试用二次插值计算(0.472)f 。
（2）已知sin()y x =的函数表试构造出差商表，利用二次Newton 插值公式计算sin1.609（保留五位有效数字）。
（3） 求不高于4次的多项式()H x ，使它满足(1)2,'(1)4,H H =-=(2)'(2)0,(3)2,H H H ===并写出余项表达式。
（4） 用最小二乘法求一个形如2y a bx =+的经验公式，使与下列数据相拟合二． 分析问题（1）题目给出概率积分的四个插值节点，要求用二次插值计算，而我们知道二次插值只需要三个插值节点即可，在该题中我们尝试取前三个点作为插值节点。
该题的重点在求插值基函数11(),(),(),k k k l x l x l x -+而它们由公式()nj n j i jj ix x l x x x =≠-=-∏(1.1)给出。
而0()()nn i i i L x y l x ==∑给出插值公式。
（2）Newton 基本插值公式00100101()()(,)()(,,,)()()...().n n n N x f x f x x x x f x x x x x x x x x =+-++--- (1.2)它的各项系数就是函数的各阶差商，每增加一个插值节点，只需要在原来的基础上多计算一项，这一性质被称作承袭性。
（3）该题需要确定一个4次插值多项式()H x ，就是要获得其各项系数的解。
4次多项式()H x 有5个系数，而题目正好给予了5个条件，这样我们会获得关于系数的非奇异五元一次方程组，在运用第三章的选列主元消元法求出各系数。
（4）本题要求用最小二乘法求经验公式。
实际上也就是确定拟合曲线的各项系数，关于系数(0,1,...,)i a i m =的线性方程组如下： 三.解决问题(1)C 语言代码如下：#include"stdio.h" #include"iostream.h" #include"math.h" #define N 3void Get_l(int k,double x,double xk[N],double L[N]); void main() { int i;double l[N],L[N],x,Lx=0;double xk[N]={0.46,0.47,0.48};double f[N]={0.4846555,0.4937542,0.5027498}; printf("输入需插值节点x:\n"); scanf("%lf",&x); for(i=0;i<N;i++) { Get_l(i,x,xk,l); L[i]=l[i]*f[i]; Lx+=L[i]; }cout.precision(7);cout<<"f(x)="<<Lx<<endl; }/*获得基函数lk(x)*/void Get_l(int i,double x,double xk[N],double l[N]) { int k;double y1=1,y2=1; for(k=0;k<N;k++)if(k!=i){ y1*=(x-xk[k]);y2*=(xk[i]-xk[k]);}l[i]=y1/y2;}结果如下：（2）C语言代码如下:#include"stdio.h"#include"math.h"#define N 3#define M 4void Put_out(int i,double a[N][M]);double Get_value(double x,double a[N][M]);void main(){ int i,j;double a[N][M],sinx,x1;double x[N]={1.5,1.6,1.7};double f[N]={0.99749,0.99957,0.99166};for(i=0;i<N;i++){ a[i][0]=x[i];a[i][1]=f[i];}for(i=0;i<N;i++)Put_out(i,a);printf("k\tf(xk)\t\tf(x0,xk)\tf(x0,x1,xk)\tf(x0,x1,x2,xk)\n");for(i=0;i<N;i++){ printf("%d\t",i);for(j=0;j<i+2;j++){ printf("%lf\t",a[i][j]);}printf("\n");}printf("Putin x:\n");scanf("%lf",&x1);sinx=Get_value(x1,a);printf("sin(%lf)=%lf\n",x1,sinx);}/*输出差商表*/void Put_out(int i,double a[N][M]){ int j;for(j=2;j<M;j++){ if(j<i)a[i][j]=(a[i][j-1]-a[0][j-1])/(a[i][0]-a[j-1][0]);if(j<i+2)a[i][j]=(a[i][j-1]-a[j-2][j-1])/(a[i][0]-a[j-2][0]);elsea[i][j]=0;}}以上代码仅是主函数和输出差商部分，求值用的是选列主元消元法，前面有所涉及，在这里不再赘述。
石大 计算方法实验报告3
二级学院： 专 业：
计算机学院
指导教师： 班级学号： 姓 名：
ห้องสมุดไป่ตู้
实验三
1、 实验目的： 在科学研究与工程技术中，常常需要从一组测量数据出发，寻找变 量的函数关系的近似表达式，使得逼近函数从总体上与已知函数的偏差 按某种方法度量能达到最小而又不一定过全部的点。这是工程中引入最 小二曲线拟合法的出发点。 充分掌握：1．最小二乘法的基本原理；2．用多项式作最小二乘曲 线拟合原理的基础上，通过编程实现一组实验数据的最小二乘拟合曲 线。 2、 实验要求: 1） 认真分析题目的条件和要求，复习相关的理论知识，选择适 当的解决方案和算法； 2） 编写上机实验程序，作好上机前的准备工作； 3）上机调试程序，并试算各种方案，记录计算的结果（包括必要的中 间结果）； 4）分析和解释计算结果；
决一些复杂烦琐的问题。 最小二乘法在编程中也有着广泛的应用，它可以很好的拟合一组比较多 的数据，而误差可以做到最小，利用它可以方便的求出拟合曲线的函数 表达式，此外，在以后的上机实验中，应该掌握一些必要的算法和程序 （比如，如何利用程序解线性方程组），对C语言进行进一步的巩固， 为以后的编程打下良好的基础。
for(k=j;k<=m;k++) X[i][k]+=X[j][k]*tmp; } } a[m]=Y[m]/X[m][m]; for(i=m-1;i>=0;i--) { a[i]=Y[i]; for(j=i+1;j<=m;j++) a[i]-=X[i][j]*a[j]; a[i]/=X[i][i]; } printf("\n 所求的幂函数为:\n"); printf("P(x)=%f*x^%f\n",exp(a[0]),a[1]); } 运行结果：
《计算方法》课内实验报告(实验2)2014
《计算方法》课内实验报告学生姓名:张靖2012309010111及学号：学院:理学院班级:信计121课程名称：计算方法实验题目:插值法与函数逼近指导教师周硕教授姓名及职称：朱振菊实验师2014年11月03日目录一、实验题目 (1)二、实验目的 (1)三、实验内容 (1)四、实验结果 (2)五、实验体会或遇到问题 (8)一、实验题目1．熟悉matlab 编写及运行数值计算程序的方法。
2．进一步理解数值积分的基础理论。
3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。
二、实验目的1．熟悉matlab 编写及运行数值计算程序的方法。
2．进一步理解插值法及函数逼近方法的理论基础。
3．进一步掌握给定数据后应用插值法及函数逼近方法进行数据处理并给出图示结果的实际操作过程。
三、实验内容1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 1⎰,要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ⎰+3021,使误差不超过510-.3．用3=n 的高斯-勒让德公式计算积分⎰31sin x e x ,给出计算结果.4．用辛普森公式 (取2==M N ) 计算二重积分.5.005.00dydx e x y ⎰⎰-四、 实验结果1．问题1：计算结果如下表表1问题1求解表复合梯形求积公式：取1210-，n=，n为迭代次数，当迭代12次后，精度达到4 n-；节点数为21=4095复合辛普森求积公式：取1000010-，节点数为n=，n为区间数，取精度为4n+=。
1100012．问题2：计算结果如下表表2问题2求解表龙贝格数值积分：给定被积函数0，被积上限3，精度为510-，龙贝格积分表中行的最大数目13，计算出龙贝格数值积分近似解为10.20759362。
3．问题3：计算结果如下表表3问题3求解表高斯-勒让德积分公式：取3n = ，节点横坐标k x 取,n k A 取585999，，，2n 阶导数e sin x x -,求得高斯-勒让德积分近似解为10.94840256。
计算方法实验报告二答案许晶
计算方法实验报告二答案许晶计算方法实验报告年级班级：11月24日计算方法实验（二）实验目的：1、利用MATLAB插值命令进行计算2、利用MATLAB 数值积分命令进行计算实验内容：1、从1点12点的11小时内，每隔1小时测量一次温度，测得的温度的数值依次为：5，8，9，15，25，29，31，30，22，25，27，24．试估计每隔1/20小时的温度值．并画出图形。
答:hours=1:12;temps=[589152529313022252724];h=1:0.05:12;t=interp1 (hours,temps,h,spline);plot(hours,temps,+,h,t,hours,temps,r:)2、图（1）是欧洲一个国家的地图，为计算出它的国土面积，首先对该国地图做出如下测量：以由西向东方向为x轴，由南向北方向为y轴，选择方便的原点，并从最西边界点到最东边界点在x轴上的区间适当的划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2,得到表格1（单位mm）根据地图的比例尺我们知道18mm 相当于40km,试由测量数据统计该国国土的近似面积，与它的精确值41288km 作比较。
图（1）答：表（1）x=[7.010.513.017.534.040.544.548.056.061.068.576.580.591.096.0101.0 104.0106.5111.5118.0123.5136.5142.0146.0150.0157.0158.0];y1=[44454 7505038303034363441454643373328326555545250666668];y2=[44597 07293100110110110117118116118118121124121121121122116838182868568];plot(x,y1,o)holdonplot(x,y2,rp)holdonxi=10:10:100;yi=interp1(x, y1,xi,spline);yi=interp1(x,y2,xi,spline);xx=min(x):0.1:max(x);yy1=interp 1(x,y1,xx,spline);plot(xx,yy1)holdonyy2=interp1(x,y2,xx,spline);plot(xx, yy2,r)holdonsymsx;s=int(-0.0107*x +1.8693*x+41.1473-0.0035*x +0. 4648*x-51.3154,7,158);s=vpa(s);s=s/18 *40 s=43819.9887242798387 649472*********、下表给出的x，y数据位于机翼剖面的轮廓线上，y1和y2分别对应轮廓的上下线。
计算方法课程实验报告
double lagrange(int n,vector<double>&x,vector<double>&y,double x){ double result=0; for(int i=0;i<n;i++){ double temp=y[i]; for(int j=0;j<n;j++){ if(i!=j){ temp = temp*(x-x[j]); temp = temp/(x[i]-x[j]); } } result += temp;
课程实验报告
课程名称：
专业班级： 学 号： 姓 名： 指导教师： 报告日期：
计算方法
信安 1302 班 U201315130
王超 付才 2015.5.2
计算机科学与技术学院
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，通力根1保过据护管生高线产中敷工资设艺料技高试术中卷0资不配料仅置试可技卷以术要解是求决指，吊机对顶组电层在气配进设置行备不继进规电行范保空高护载中高与资中带料资负试料荷卷试下问卷高题总中2体2资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况1卷中下安，与全要过，加度并强工且看作尽护下可1都关能可于地以管缩正路小常高故工中障作资高；料中对试资于卷料继连试电接卷保管破护口坏进处范行理围整高，核中或对资者定料对值试某，卷些审弯异核扁常与度高校固中对定资图盒料纸位试，置卷编.工保写况护复进层杂行防设自腐备动跨与处接装理地置，线高尤弯中其曲资要半料避径试免标卷错高调误等试高，方中要案资求，料技编试术写5、卷交重电保底要气护。设设装管备备置线4高、调动敷中电试作设资气高，技料课中并3术试、件资且中卷管中料拒包试路调试绝含验敷试卷动线方设技作槽案技术，、以术来管及避架系免等统不多启必项动要方高式案中，；资为对料解整试决套卷高启突中动然语过停文程机电中。气高因课中此件资，中料电管试力壁卷高薄电中、气资接设料口备试不进卷严行保等调护问试装题工置，作调合并试理且技利进术用行，管过要线关求敷运电设行力技高保术中护。资装线料置缆试做敷卷到设技准原术确则指灵：导活在。。分对对线于于盒调差处试动，过保当程护不中装同高置电中高压资中回料资路试料交卷试叉技卷时术调，问试应题技采，术用作是金为指属调发隔试电板人机进员一行，变隔需压开要器处在组理事在；前发同掌生一握内线图部槽 纸故内资障，料时强、，电设需回备要路制进须造行同厂外时家部切出电断具源习高高题中中电资资源料料，试试线卷卷缆试切敷验除设报从完告而毕与采，相用要关高进技中行术资检资料查料试和，卷检并主测且要处了保理解护。现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
《计算方法》实验报告材料
《计算方法》实验报告材料引言：计算方法是一门应用数学的基础课程，通过实验教学，能够帮助学生更好地理解和掌握各种数值计算的方法和技巧。
本次实验旨在通过编程实现一些常用的数值计算方法，并通过对比分析实验结果，验证方法的有效性和可行性。
实验一：插值算法插值算法是利用已知的数据点，构建一个连续函数以逼近数据的方法。
本次实验中使用的插值算法为拉格朗日插值和牛顿插值。
通过编程实现这两种算法，并选取若干个数据点进行测试，得到插值函数的结果。
通过比较原始数据和插值函数的结果，可以验证插值算法的准确性和可行性。
实验二：方程求解方程求解是数值计算中的一个重要问题，求解非线性方程、线性方程组和特征值问题等都需要采用相应的迭代方法。
本次实验中，我们实现了常用的牛顿迭代法和二分法，并选择数学问题进行求解。
通过比较实验结果和理论值的误差，可以验证求解方法的精确性和可行性。
实验三：数值积分数值积分是利用数值方法对定积分进行近似求解的过程。
本次实验中，我们实现了矩形法、梯形法和辛普森法等常用的数值积分方法，并选取若干函数进行数值积分的计算。
通过比较数值积分的结果和解析解或数值解的误差，可以验证数值积分方法的准确性和可行性。
实验四：常微分方程求解常微分方程求解是数值计算中的一个重要问题，常常涉及到物理、化学、生物等科学领域。
本次实验中，我们实现了欧拉方法和龙格-库塔方法等常用的常微分方程求解算法，并选取若干常微分方程进行求解。
通过比较数值解和解析解或数值解的误差，可以验证常微分方程求解方法的精确性和可行性。
实验五：线性方程组求解线性方程组求解是数值计算中的一个重要问题，常常涉及到矩阵的运算和迭代方法。
本次实验中，我们实现了高斯消元法和追赶法等常用的线性方程组求解算法，并选择一些矩阵进行求解。
通过比较数值解和解析解或数值解的误差，可以验证线性方程组求解方法的精确性和可行性。
结论：通过本次实验，我们掌握了插值算法、方程求解、数值积分、常微分方程求解和线性方程组求解等常用的计算方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
b.^3
b*b*b
1、 利用循环语句编写程序，求 1^2+2^2+3^2+…+100^2。
2、 利用选择语句编写函数文件 demo3 实现 sgn 函数功能，并计算 demo3(0) demo3(90) demo3(-12)
function f=demo3(x) if x>0
f=1; elseif x==0
计算方法实验报告
年级班级： 学 号： 姓 名： 日 期：
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
ans =
1 3 5 7 9 11 2 4 6 8 10 12 2、已知 a=[1 2 3;2 3 4;3 4 5]; b=[1 1 1;2 2 2;3 3 3]; 在计算机上运行下面 的命令，并写出计算结果。a.*b a*b a.\b a./b
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
计AB 控制语句 2、 熟悉数组运算 3、 MATLAB 图形处理功能 4、 MATLAB 程序初步设计 实验内容： 一、 数组运算 1、已知向量 a=1:12, 用 reshape 命令将其编程 2*6 矩阵。 >> a=1:12; >> reshape(a,2,6)
三、MATLAB图形处理（写出程序，给出图形） 1、在0≤x≤2区间内，绘制曲线 y=2e^(0.5x)cos(4x)
2、在[0, 20]之间，绘制三维曲线 x=sin(t); y=cos(t); z=tsin(t)cos(t); 图形要求：标题为“'Line in 3-D Space”，x轴，y轴, z轴分别添加标注 “X”“Y”“Z”。
3、已知 b=[1 1 1;2 2 2;3 3 3]; 在计算机上运行 b^3 并写出计算结果。
二、程序设计（写出程序，并给出计算结果）
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
f=0; else x<0
f=-1; end
demo3(0)
ans =
0 demo3(90)
ans =
1 demo3(90)
ans = -1
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。