文档视界 最新最全的文档下载
当前位置:文档视界 › 北京理工大学数值计算方法大作业数值实验1

北京理工大学数值计算方法大作业数值实验1

北京理工大学数值计算方法大作业数值实验1
北京理工大学数值计算方法大作业数值实验1

数值实验(一)

(一)书P14/4

分别将区间[10,10]

分为100,200,400等份,利用mesh或surf命令画出二元函数的三维图形。

Z=

【MATLAB求解】

[x,y]=meshgrid(-10:0.1:10);

a=exp(-abs(x));

b=cos(x+y);

c=1./(x.^2+y.^2+1);

z=a+b+c;

mesh(x,y,z);

[x,y]=meshgrid(-10:0.05:10); a=exp(-abs(x));

b=cos(x+y);

c=1./(x.^2+y.^2+1);

z=a+b+c;

mesh(x,y,z);

[x,y]=meshgrid(-10:0.025:10); a=exp(-abs(x));

b=cos(x+y);

c=1./(x.^2+y.^2+1);

z=a+b+c;

mesh(x,y,z);

(二)书P7/1.3.2数值计算的稳定性(I)C语言程序—不稳定解

取=ln1.2,按公式(n=1,2,…)

#include

#include

#include

void main()

{

float m=log(6.0)-log(5.0),n;

int i;

printf("y[0]=%-20f",m);

while(i<20)

{

n=1/i-5*m;

printf("y[%d]=%-20f",i,n);

m=n;

i++;

if(i%3==0)printf("\n");

}

getch();

}

(II)C语言程序—稳定解

按公式(k=n,n-1,n-2,…)

#include

#include

#include

void main()

{

float m=(1/105.0+1/126.0)/2,n;

int i;

printf("y[50]=%-20f",m);

while(i>1)

{

n=1/(5.0*i)-m/5.0;

printf("y[%d]=%-20f",i-1,n);

m=n;

i--;

if(i%3==0)printf("\n");

}

getch();

}

北京理工大学远程教育学院 操作系统 2019-2020学年

(93) 北京理工大学远程教育学院2019-2020学年第二学期 《操作系统》期末试卷(B卷) 教学站学号姓名成绩 应用题(每题20分,共100分) 1.批处理操作系统、分时操作系统和实时操作系统的特点各是什么? 2. 动态分区分配的基本分配思想是什么?常用的分配算法有哪四种,简述其含义。 3.为什么要引入缓冲区? 4.解释产生死锁的四个必要条件是什么,当预防死锁时可以通过破坏这四个条件的任意一个来实现,请分别论述破坏每个条件会产生什么影响,哪些条件能破坏哪些条件不能破坏? 5.某工厂有一个可以存放设备的仓库,总共有8个位置可以存放8台设备。生产部门生产的每一台设备都必须入库。销售部门可以从仓库提出设备供应客户。设备的出库和入库都必须借助运输工具。现在只有一套运输工具,每次只能运输一台设备,系统共使用三个信号量,S代表互斥信号量,表示运输工具;S1和S2均为同步信号量,S1表示仓库中可以存放设备的空闲位置,S2表示仓库中已经被设备占用了的位置。请设计一个能协调工作的自动调度管理系统,并利用记录型信号量写出解决此问题的程序代码,请注明信号量的初值。

(编号) 北京理工大学远程教育学院2019-2020学年第二学期 《操作系统》期末试卷(B卷)答题纸 教学站学号姓名成绩 1.答:1)批处理操作系统的用户脱机使用计算机,作业是成批处理的,系统内多道程序并发执行,交互能力差。 2)分时操作系统可以让多个用户同时使用计算机,人机交互性较强,具有每个用户独立使用计算机的独占性,系统响应及时。 3)实时操作系统能对控制对象作出及时反应,可靠性高,响应及时,但是资源利用率低。 2. 答:算法思想:将空闲分区链以地址递增的顺序连接;在进行内存分配时,从链首开始顺序查找,直到找到一块分区的大小可以满足需求时,按照该作业的大小,从该分区中分配出内存,将剩下的空闲分区仍然链在空闲分区链中。 第一种:首次适应算法(First fit) 按空闲分区依地址递增次序链接,分配内存时按顺序查找,放入第一个匹配到的空闲分区,会造成内部碎片,有着较大的浪费 第二种:最佳适应算法(Best fit) 将空闲分区按内存大小递增的顺序链接起来,分配内存时按照顺序放入第一个匹配的空闲分区。 第三种:最坏适应算法(Worst fit) 将空闲分区按容量递减的顺序链接起来,分配内存时放入第一个匹配的空闲分区,即最大的分区,造成内部碎片

数值计算实验课题目

数值实验课试题 本次数值实验课结课作业,请按题目要求内容写一篇文章。按题目要求 人数自由组合,每组所选题目不得相同(有特别注明的题目除外)。试题如下: 1)解线性方程组的Gauss 消去法和列主元Gauss 消去法(2人)/*张思珍,巩艳华*/ 用C 语言将不选主元和列主元Gauss 消去法编写成通用的子程序,然后用你编写的程序求解下列84阶的方程组 ???? ?????? ? ??=??????????? ????????????? ? ?1415151515768 168 168 168 1681684 8382321 x x x x x x 参考书目: 1.《计算机数值方法》,施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》,徐树方、高立、张平文编 3.《数值分析简明教程》,王能超编 2)解线性方程组的平方根法(4人)/*朱春成、黄锐奇、张重威、章杰*/ 用C 语言将平方根法和改进的平方根法编写成通用的子程序,然后用你编写的程序求解对称正定方程组b Ax =,其中 (1)b 随机的选取,系数矩阵为100阶矩阵 ?????? ???? ? ? ?101 1101 1101 1101 1101110 ; (2)系数矩阵为40阶的Hilbert 矩阵,即系数矩阵A 的第i 行第j 列元素为 1 1-+= j i a ij ,向量b 的第i 个分量为∑=-+ = n j i j i b 1 1 1. 参考书目: 1.《计算机数值方法》,施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》,徐树方、高立、张平文编

3.《数值分析简明教程》,王能超编 3)三对角线方程组的追赶法(3人)/*黄佳礼、唐伟、韦锡倍*/ 用C 语言将三对角线方程组的追赶法法编写成通用的子程序,然后用你编写的程序求解如下84阶三对角线方程组 ???? ?????? ? ??=??????????? ????????????? ? ?1415151515768 168 168 168 16816 84 8382321 x x x x x x 参考书目: 1.《计算机数值方法》,施吉林、刘淑珍、陈桂芝编 2.《数值分析简明教程》,王能超编 4)线性方程组的Jacobi 迭代法(3人)/*周桂宇、杨飞、李文军*/ 用C 语言将Jacobi 迭代法编写成独立的子程序,并用此求解下列方程组, 精确到小数点后5位 ???? ? ??=????? ??????? ? ?-149012 2111221 3 2 1 x x x 参考书目: 1.《计算机数值方法》,施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》,徐树方、高立、张平文编 3.《数值分析简明教程》,王能超编 5)线性方程组的Gauss-Seidel 迭代法(3人)/*张玉超、范守平、周红春*/ 用C 语言将Gauss-Seidel 迭代法编写成独立的子程序,并用此求解下列方程组,精确到小数点后5位 ???? ? ??=????? ??????? ? ?--39721 1111112 3 2 1 x x x 参考书目: 1.《计算机数值方法》,施吉林、刘淑珍、陈桂芝编 2.《数值线性代数》,徐树方、高立、张平文编 3.《数值分析简明教程》,王能超编 6)解线性方程组的最速下降法法(2人)/*赵育辉、阿热孜古丽*/ 用C 语言将最速下降法编写成通用的子程序,然后用你编写的程序求解对称

数值分析大作业-三、四、五、六、七

大作业 三 1. 给定初值 0x 及容许误差 ,编制牛顿法解方程f (x )=0的通用 程序. 解:Matlab 程序如下: 函数m 文件:fu.m function Fu=fu(x) Fu=x^3/3-x; end 函数m 文件:dfu.m function Fu=dfu(x) Fu=x^2-1; end 用Newton 法求根的通用程序Newton.m clear; x0=input('请输入初值x0:'); ep=input('请输入容许误差:'); flag=1; while flag==1 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)

while flag1==1 && m<=10^3 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)=ep flag=0; end end fprintf('最大的sigma 值为:%f\n',sigma); 2.求下列方程的非零根 5130.6651()ln 05130.665114000.0918 x x f x x +??=-= ?-???解:Matlab 程序为: (1)主程序 clear clc format long x0=765; N=100; errorlim=10^(-5); x=x0-f(x0)/subs(df(),x0); n=1; while nerrorlim n=n+1; else break ; end x0=x; end disp(['迭代次数: n=',num2str(n)]) disp(['所求非零根: 正根x1=',num2str(x),' 负根x2=',num2str(-x)]) (2)子函数 非线性函数f function y=f(x) y=log((513+0.6651*x)/(513-0.6651*x))-x/(1400*0.0918); end

曲线拟合的数值计算方法实验

曲线拟合的数值计算方法实验 【摘要】实际工作中,变量间未必都有线性关系,如服药后血药浓度与时间的关系;疾病疗效与疗程长短的关系;毒物剂量与致死率的关系等常呈曲线关系。曲线拟合(curve fitting)是指选择适当的曲线类型来拟合观测数据,并用拟合的曲线方程分析两变量间的关系。曲线直线化是曲线拟合的重要手段之一。对于某些非线性的资料可以通过简单的变量变换使之直线化,这样就可以按最小二乘法原理求出变换后变量的直线方程,在实际工作中常利用此直线方程绘制资料的标准工作曲线,同时根据需要可将此直线方程还原为曲线方程,实现对资料的曲线拟合。常用的曲线拟合有最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束。 关键词曲线拟合、最小二乘法拟合、幂函数拟合、对数函数拟合、线性插值、三次样条插值、端点约束 一、实验目的 1.掌握曲线拟合方式及其常用函数指数函数、幂函数、对数函数的拟合。 2.掌握最小二乘法、线性插值、三次样条插值、端点约束等。 3.掌握实现曲线拟合的编程技巧。 二、实验原理 1.曲线拟合 曲线拟合是平面上离散点组所表示的坐标之间的函数关系的一种数据处理方法。用解析表达式逼近离散数据的一种方法。在科学实验或社会活动中,通过 实验或观测得到量x与y的一组数据对(X i ,Y i )(i=1,2,...m),其中各X i 是彼此不同的。人们希望用一类与数据的背景材料规律相适应的解析表达式,y=f(x,c)来反映量x与y之间的依赖关系,即在一定意义下“最佳”地逼近或 拟合已知数据。f(x,c)常称作拟合模型,式中c=(c 1,c 2 ,…c n )是一些待定参 数。当c在f中线性出现时,称为线性模型,否则称为非线性模型。有许多衡量拟合优度的标准,最常用的一种做法是选择参数c使得拟合模型与实际观测值在

数值分析大作业三 四 五 六 七

大作业 三 1. 给定初值 0x 及容许误差 ,编制牛顿法解方程f (x )=0的通用程序. 解:Matlab 程序如下: 函数m 文件:fu.m function Fu=fu(x) Fu=x^3/3-x; end 函数m 文件:dfu.m function Fu=dfu(x) Fu=x^2-1; end 用Newton 法求根的通用程序Newton.m clear; x0=input('请输入初值x0:'); ep=input('请输入容许误差:');

flag=1; while flag==1 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)

while flag==1 sigma=k*eps; x0=sigma; k=k+1; m=0; flag1=1; while flag1==1 && m<=10^3 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)=ep flag=0;

end end fprintf('最大的sigma 值为:%f\n',sigma); 2.求下列方程的非零根 5130.6651()ln 05130.665114000.0918 x x f x x +?? =-= ?-???解: Matlab 程序为: (1)主程序 clear clc format long x0=765; N=100; errorlim=10^(-5); x=x0-f(x0)/subs(df(),x0); n=1;

17春北理工《操作系统》在线作业

2017秋17春北理工《操作系统》在线作业 一、单选题(共20 道试题,共40 分。) 1. 操作系统中有一类程序,其执行过程不可中断,该程序叫做:( ) A. 过程 B. 子程序 C. 原语 D. 模块机制 正确答案: 2. 操作系统向用户提供了使用计算机的程序一级的接口为() A. 原语 B. 进程管理 C. 命令接口 D. 系统调用 正确答案: 3. 请求分页系统管理中,若把页面的尺寸增加一倍,程序顺序执行时,其缺页中断次数一般会:( ) A. 增加 B. 减少 C. 不变 D. 可能增加,也可能减少 正确答案: 4. 以下关于主存空间的说法中正确的是() A. 主存储器的空间分成三部分:系统区、用户区和缓冲区 B. 操作系统与硬件的接口信息、操作系统的管理信息和程序等存放在主存储器的系统区 C. 所有的程序存放在用户区 D. 存储管理是对主存储器空间的各部分如系统区、用户区等进行管理 正确答案: 5. 在可变式分区分配方案中,最佳适应算法是将空闲区在空闲区表中按()次序排列 A. 容量递增 B. 容量递减 C. 地址递增 D. 地址递减 正确答案: 6. 一个功能强的文件系统,向用户提供更加灵活的文件物理结构是:( ) A. 连续结构 B. 串联结构

C. 索引结构 D. 三者都不对 正确答案: 7. UNIX系统中,()是实现把一个进程的输出连接到另一个进程的输入功能的机制。 A. 普通文件 B. 特殊文件 C. 目录文件 D. 管道文件 正确答案: 8. 虚拟存储器的最大容量() A. 为内外存容量之和 B. 由计算机的地址结构决定 C. 是任意的 D. 由作业的地址空间决定 正确答案: 9. 以下叙述中,不正确的是() A. 采用动态重定位,在必要时可以改变装入的作业在主存中的存放区域 B. 采用动态重定位的系统支持“程序浮动” C. 采用静态重定位的系统不支持“程序浮动” D. 采用动态和静态重定位的系统都支持“程序浮动” 正确答案: 10. 吞吐量是指:( ) A. 单位时间内完成的信息量 B. 操作系统响应进程命令需要的信息量 C. 完成作业或进程所需要的信息量 D. 都不对 正确答案: 11. 计算机系统的二级存储包括() A. 主存储器和辅助存储器 B. ROM和RAM C. 超高速缓存和内存储器 D. CPU寄存器和主存缓冲区 正确答案: 12. 最佳适应算法通常是将空闲区按()排列 A. 地址大到小 B. 地址小到大 C. 空间大到小 D. 空间小到大 正确答案: 13. 适合多道程序运行的最简单的方案是( ) A. 分页式存储管理 B. 固定分区式存储管理 C. 分段式存储管理

数值分析实验报告1

实验一误差分析 实验1.1(病态问题) 实验目的:算法有“优”与“劣”之分,问题也有“好”与“坏”之别。对数值方法的研究而言,所谓坏问题就是问题本身对扰动敏感者,反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中,如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决,当然一般要付出一些代价(如耗用更多的机器时间、占用更多的存储空间等)。 问题提出:考虑一个高次的代数多项式 显然该多项式的全部根为1,2,…,20共计20个,且每个根都是单重的。现考虑该多项式的一个扰动 其中ε(1.1)和(1.221,,,a a 的输出b ”和“poly ε。 (1(2 (3)写成展 关于α solve 来提高解的精确度,这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。 实验过程: 程序: a=poly(1:20); rr=roots(a); forn=2:21 n form=1:9 ess=10^(-6-m);

ve=zeros(1,21); ve(n)=ess; r=roots(a+ve); -6-m s=max(abs(r-rr)) end end 利用符号函数:(思考题一)a=poly(1:20); y=poly2sym(a); rr=solve(y) n

很容易的得出对一个多次的代数多项式的其中某一项进行很小的扰动,对其多项式的根会有一定的扰动的,所以对于这类病态问题可以借助于MATLAB来进行问题的分析。 学号:06450210 姓名:万轩 实验二插值法

数值计算方法大作业

目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20)

4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28)

数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码: Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例:求f(x)=e2x-6x=0在x=0.5附近的根(ep=10-10)

1.2牛顿法 程序代码: Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例:求56的值。(ep=10-10)

北理工应用文在线作业答案

北京理工大学应用文写作在线作业答案 北理工应用文在线作业答案 第1阶段第一阶段 第2阶段第二阶段 计划的标题同许多事务文书一样,不必写明时限 写总结不一定要按照完成工作的时间先后顺序来写 调查报告可以用于向上级机关汇报工作,供领导决策参考 对未来一定时期的任务作出预想性安排的文种是 总结的开头包括的内容,下列不准确的一项是 调查报告在格式上没有固定的要求,一般包括 计划的种类很多,而且从不同的角度可以对其进行不同的分类,如按其内容分,可将其分为

综合性计划和 总结的正文一般包括基本情况、成绩收获、____________、经验体会等几部分内容 关于计划的主要作用,以下叙述正确的是 总结主体的主要内容包括 条款式合同适用于工程承包、科技合作、合作生产、技术引进等内容比较的经济合同。 招标、投标文书最为突出的两大特点是竞争性和公开性。 我国《广告法》中所称的广告包括商业广告和非商业广告。 下列计划标题拟定恰当的一项是 写作调查报告,在表达方式上,要做到 一篇演讲稿的结尾以极富鼓动性的言辞号召人们为某种目的、某种理想而行动起来。这种结尾的方式叫 经济活动分析报告的标题各项内容中不能省去的一项是

是审计机构或审计人员在完成某一项审计工作后,向委托者或授权者提交的情况书面报告。投标书一般是由_____________设计并送给投标单位的。 审计报告写作应注意的事项 第3阶段第三阶段 个人请柬应一人一柬,夫妻也不可合写一柬。 撰写欢迎词要大量选择感情色彩浓烈,感染力量强大的形容词、比喻词、象征词。 启事的标题可以只用事由表示。 消息的第一自然段或开头的一两句话,一般被称作(),它将消息最重要、最新鲜的事实概括出来,并吸引读者。 通讯写作首先要注意 不管是欢迎词表达“有朋自远方来,不亦乐乎”的愉悦心情,还是欢送词表达亲朋远行的依依惜别之情,都具有的特点是( )。 着重记述社会变化、风土人情和建设状况,并在报纸上常以“巡礼”、“侧记”等形式出现的新闻体裁是____________。 广播稿主要使用语言来影响听众,所以_________是其最大特色。 消息有不同的划分方法,如从写作的角度来划分,可分为()四类。 请柬结尾的礼貌用语有 一般书信常用于个人之间的交往,也可以是个人写给单位或集体的 申请书内容比较单纯,一般一事一书 介绍信是用来介绍联系接洽事宜的一种应用文体,它只具有介绍的作用 中央电视台的新闻联播属于 下列材料不适合写贺信的是 下面关于悼词写作不正确的一项是 书信的种类很多,按使用目的和范围可分一般书信和( )两大类 下列语句叙述有错误的有 报告要反映工作实践,是可以运用某些修辞手法写得生动些,但有些修辞手法是不能使用的。

太原理工大学数值计算方法实验报告

本科实验报告 课程名称:计算机数值方法 实验项目:方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点:行勉楼 专业班级: ******** 学号: ********* 学生姓名: ******** 指导教师:李誌,崔冬华 2016年 4 月 8 日

y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为:f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法: // 方程求根(割线法).cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"

心得体会 使用不同的方法,可以不同程度的求得方程的解,通过二分法计算的程序实现更加了解二分法的特点,二分法过程简单,程序容易实现,但该方法收敛比较慢一般用于求根的初始近似值,不同的方法速度不同。面对一个复杂的问题,要学会简化处理步骤,分步骤一点一点的循序处理,只有这样,才能高效的解决一个复杂问题。

数值计算方法实验5

实验报告 学院(系)名称: 主程序部分列选主元部分

实验结果: 一.列主元消去法 输入各个数据,最终使用列选主元法,得到结果为:x1=x2=x3=1二.高斯-赛德尔迭代法 输入各个数据,输出每一步迭代数据,最终结果为:x1=0.285716,附录(源程序及运行结果) 一.列主元高斯消去法 #include #include void print(double a[3][3],int n,double b[3]){ printf("输出矩阵:\n"); for(int i=0;ifabs(d)){ d=a[i][k]; l=i; } i++; } printf("选出主元:%lf\n",d); if(d==0) printf("矩阵奇异!\n"); else if(l!=k){ for(int j=k;j

数值分析大作业

数值分析报大作业 班级:铁道2班 专业:道路与铁道工程 姓名:蔡敦锦 学号:13011260

一、序言 该数值分析大作业是通过C语言程序编程在Microsoft Visual C++ 6.0编程软件上运行实现的。本来是打算用Matlab软间来计算非线性方程的根的。学习Matlab也差不多有一个多月了,感觉自己编程做题应该没什么问题了;但是当自己真心的去编程、运行时才发现有很多错误,花了一天时间修改、调试程序都没能得到自己满意的结果。所以,我选择了自己比较熟悉的C程序语言来编程解决非线性的求值问题,由于本作业是为了比较几种方法求值问题的收敛速度和精度的差异,选择了一个相对常见的非线性函数来反映其差异,程序运行所得结果我个人比较满意。编写C语言,感觉比较上手,程序出现问题也能比较熟练的解决。最终就决定上交一份C程序语言编程的求值程序了!

二、选题 本作业的目的是为了加深对非线性方程求根方法的二分法、简单迭代法、、牛顿迭代法弦截法等的构造过程的理解;能将各种方法的算法描述正确并且能够改编为程序并在计算机上实现程序的正确合理的运行,能得到自己满意的结果,并且能调试修改程序中可能出现的问题和程序功能的增减修改。本次程序是为了比较各种方法在求解同一非线性方程根时,在收敛情况上的差异。 为了达到上面的条件我选择自己比较熟悉的语言—C语言来编程,所选题目为计算方程f(x)=x3-2x-5=0在区间[2,3]内其最后两近似值的差的绝对值小于等于5 ?的根的几种方法的比较。 110- 本文将二分法、牛顿法、简单迭代法、弦截法及加速收敛法这五种方法在同一个程序中以函数调用的方式来实现,比较简洁明了,所得结果能很好的比较,便于分析;发现问题和得出结论。

北京理工大学现代远程教育学院

北京理工大学现代远程教育学院 毕业设计(论文) 层次专升本 专业工商管理 姓名顾强 毕业时间 2003年 7月 ________________________________________________________________ 1

北京理工大学 现代远程教育学院 毕业设计任务书 专业名称工商管理 班号(教学站)北京船舶工业管理干部学院 姓名顾强 毕业设计题目 浅析企业合同管理中的风险控制 指导教师陈琪 负责人签字陈琪 2003 年6月8日 ________________________________________________________________ 2

内容和要求:在日常的社会经济活动中,无论是企业还是个人都可能面临风险。所以说风险无处不在,作为企业主要经济活动之一的合同管理也不例外。 但是,在我国社会主义市场经济体制不断完善的过程中,随着社会经济活动的日益活跃,由于信用体制的不健全,企业在生产经营活动中所面临的风险越来越多,甚至风险损失使一些企业到了难以为继的地步。所以,企业唯有从内部管理做起,提高自身的风险防范能力,才能保持健康的发展。 由于风险管理的概念引入我国的时间不长,在生产型企业的合同管理中还远未得到普及和应用。而且,由于在不同的行业或应用领域,风险具有不同的表现形式和特征,无论是金融、投资风险管理的成熟经验,还是目前热门的项目管理均不能在企业的合同管理中照搬套用。为此,本文在对风险与风险管理做了一般概述的基础上,根据生产型企业的特点,结合所学知识和实践经验,重点剖析了企业在合同管理过程中所面临的种种风险及其表现形式,并从实际应用的角度出发,在权衡风险管理成本等因素的前提下,对如何在合同管理活动中加强风险的防范措施进行了系统、全面的探讨。提出了优化业务流程、完善风险管理制度、建立风险管理机制等风险防范和应对措施。希望借此能帮助 ________________________________________________________________ 3

北京理工大学网络教育期末考试大学英语二 3

大学英语二-2005-0003 (题目数量:29 总分:100.0) 1.单选题(题目数量:20 总分:40.0) 1. She’s been reported missing and we’re checking her movements for the two weeks before she ## . A.vanished B.survived C.observed D.revealed 答案:A 2. —Tom told me that you collect stamps. —##. A.at do you mean? B.I don't think so. C.Why did he tell you? D.Yes, I do. Do you? 答案:D 3. “Is there anything in your personal life which would ## the party if it became public?” A.embarrass B.disturb C.waken D.involve 答案:A

4. They provide information on methods to improve the ## training programs for older workers, including helping older workers to find productive employment. A.formation B.efficiency C.selection D.dignity 答案:B 5. If you don’t read, then ## you will forget how to read. A.mysteriously B.unfortunately C.positively D.eventually 答案:D 6. —Could you tell me where the nearest parking lot is? —## A.Don't ask me. B.Parking is very difficult. C.Sorry, I'm a stranger here. D.Sorry, you can't park here 答案: C 7. Burns can be caused when skin that isn’t used to the sun is ## to short bursts of strong sunlight. A.exposed B.snapped C.wrinkled D.revealed 答案:A

(完整版)数值计算方法上机实习题答案

1. 设?+=1 05dx x x I n n , (1) 由递推公式n I I n n 1 51+-=-,从0I 的几个近似值出发,计算20I ; 解:易得:0I =ln6-ln5=0.1823, 程序为: I=0.182; for n=1:20 I=(-5)*I+1/n; end I 输出结果为:20I = -3.0666e+010 (2) 粗糙估计20I ,用n I I n n 51 5111+- =--,计算0I ; 因为 0095.05 6 0079.01020 201 020 ≈<<≈??dx x I dx x 所以取0087.0)0095.00079.0(2 1 20=+= I 程序为:I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); end I 0I = 0.0083 (3) 分析结果的可靠性及产生此现象的原因(重点分析原因)。 首先分析两种递推式的误差;设第一递推式中开始时的误差为000I I E '-=,递推过程的舍入误差不计。并记n n n I I E '-=,则有01)5(5E E E n n n -==-=-Λ。因为=20E 20020)5(I E >>-,所此递推式不可靠。而在第二种递推式中n n E E E )5 1(5110-==-=Λ,误差在缩小, 所以此递推式是可靠的。出现以上运行结果的主要原因是在构造递推式过程中,考虑误差是否得到控制, 即算法是否数值稳定。 2. 求方程0210=-+x e x 的近似根,要求4 1105-+?<-k k x x ,并比较计算量。 (1) 在[0,1]上用二分法; 程序:a=0;b=1.0; while abs(b-a)>5*1e-4 c=(b+a)/2;

数值分析实验报告总结

数值分析实验报告总结 随着电子计算机的普及与发展,科学计算已成为现代科 学的重要组成部分,因而数值计算方法的内容也愈来愈广泛和丰富。通过本学期的学习,主要掌握了一些数值方法的基本原理、具体算法,并通过编程在计算机上来实现这些算法。 算法算法是指由基本算术运算及运算顺序的规定构成的完 整的解题步骤。算法可以使用框图、算法语言、数学语言、自然语言来进行描述。具有的特征:正确性、有穷性、适用范围广、运算工作量少、使用资源少、逻辑结构简单、便于实现、计算结果可靠。 误差 计算机的计算结果通常是近似的,因此算法必有误差, 并且应能估计误差。误差是指近似值与真正值之差。绝对误差是指近似值与真正值之差或差的绝对值;相对误差:是指近似值与真正值之比或比的绝对值。误差来源见表 第三章泛函分析泛函分析概要 泛函分析是研究“函数的函数”、函数空间和它们之间 变换的一门较新的数学分支,隶属分析数学。它以各种学科

如果 a 是相容范数,且任何满足 为具体背景,在集合的基础上,把客观世界中的研究对象抽 范数 范数,是具有“长度”概念的函数。在线性代数、泛函 分析及相关的数学领域,泛函是一个函数,其为矢量空间内 的所有矢量赋予非零的正长度或大小。这里以 Cn 空间为例, Rn 空间类似。最常用的范数就是 P-范数。那么 当P 取1, 2 ,s 的时候分别是以下几种最简单的情形: 其中2-范数就是通常意义下的距离。 对于这些范数有以下不等式: 1 < n1/2 另外,若p 和q 是赫德尔共轭指标,即 1/p+1/q=1 么有赫德尔不等式: II = ||xH*y| 当p=q=2时就是柯西-许瓦兹不等式 般来讲矩阵范数除了正定性,齐次性和三角不等式之 矩阵范数通常也称为相容范数。 象为元素和空间。女口:距离空间,赋范线性空间, 内积空间。 1-范数: 1= x1 + x2 +?+ xn 2-范数: x 2=1/2 8 -范数: 8 =max oo ,那 外,还规定其必须满足相容性: 所以

北京理工大学网络教育期末考试大学英语二 1

大学英语二-2005-0001 (题目数量:29 总分:100.0) 1.单选题(题目数量:20 总分:40.0) 1. —Sorry, I couldn’t come to the party. I was sick that day. —##. A.I don’t know that B.That’s all right C.Yes, we’ll have another party D.No, the party wasn’t held 答案:B 2. But it is ## that publishers, to stay in business, need to bring out a large number of titles every year. A.spontaneous B.boring C.elementary D.evident 答案:D 3. The fight lasted but few minutes, when the enemy soldiers ## before us for twelve miles. A.vanished B.retreated C.approached D.proceeded 答案:B

4. If people are made aware that the better the service, the more they earn, then it really is an ## to work hard and perform well. A.incentive B.atmosphere https://www.docsj.com/doc/f02931831.html,munity D.gratitude 答案:A 5. Eye play an important role in body balance and it is important to ## as much information from the environment as possible. A.absorb B.involve C.invent D.conquer 答案:A 6. —Let me help you carry the suitcase. —##. A.It's OK. I can manage. B.It's not very light. C.I can help you with it. D.Put it down on the ground. 答案:A 7. Mr. Trump would no t## details of his plan for a new factory to make mobile phones in China. A.expose B.reveal C.indicate D.enclose 答案:B

上海大学_王培康_数值分析大作业

数值分析大作业(2013年5月) 金洋洋(12721512),机自系 1.下列各数都是经过四舍五入得到的近似值,试分别指出它 们的绝对误差限, 相对误差限和有效数字的位数。 X1 =5.420, x 2 =0.5420, x 3=0.00542, x 4 =6000, x 5=50.610? 解:根据定义:如果*x 的绝对误差限 不超过x 的某个数位的半个单位,则从*x 的首位非零数字到该位都是有效数字。 显然根据四舍五入原则得到的近视值,全部都是有效数字。 因而在这里有:n1=4, n2=4, n3=3, n4=4, n5=1 (n 表示x 有效数字的位数) 对x1:有a1=5, m1=1 (其中a1表示x 的首位非零数字,m1表示x1的整数位数) 所以有绝对误差限 143 11 (1)101022 x ε--≤ ?=? 相对误差限 31() 0.510(1)0.00923%5.4201 r x x x εε-?= == 对x2:有a2=5, m2=0 所以有绝对误差限 044 11 (2)101022 x ε--≤ ?=? 相对误差限 42() 0.510(2)0.00923%0.54202 r x x x εε-?= == 对x3:有a3=5, m3=-2 所以有绝对误差限 235 11 (3)101022 x ε---≤ ?=? 相对误差限 53() 0.510(3)0.0923%0.005423 r x x x εε-?= == 对x4:有a4=0, m4=4 所以有绝对误差限 4411(4)1022 x ε-≤?= 相对误差限 4() 0.5 (4)0.0083%6000 4 r x x x εε= = = 对x5:有a5=6, m5=5 所以有绝对误差限 514 11(5)101022 x ε-≤ ?=? 相对误差限 45() 0.510(5)8.3%600005 r x x x εε?= ==

数值分析实验报告-Sor法分析

数值分析实验报告 一、 实验目的 1、会使用Sor 法求解一个线性方程组 2、熟悉matlab 语言并结合原理编程求方程组 3、改变ω的值观察实验结果 4、会分析实验结果 二、实验题目 编制Sor 迭代格式程序进行求解一个线性方程组的迭代计算情况,运行中要选用不同的松弛因子ω进行尝试 三、 实验原理 Jacobi 迭代和seidel 迭代对具体的线性方程组来说,逼近*x 的速度是固定不变的,遇到收敛很慢的情况时就显得很不实用。 Sor 法是一seidel 迭代为基础,并在迭代中引入参数ω以增加迭代选择的灵活性,具体为: ! 用seidel 迭代算出的,)()1()()1(k k J k k J x x x x x -=?++相减得到差向量与再用参数ω乘之再加上 )1()()()1()1()()()1(++++-=?+=k J k k k k k k x x x x x x x x ωωω,即的下一步迭代作为,由seidel 迭代的公式可以得到Sor 法的迭代格式为 n i x a x a b a x x k j n i j ij k j i j ij i ii k i k i ,2,1),()1()(1)1(11)()1( =--+-=∑∑+=+-=+ω ω 式中ω称为松弛因子。 四、 实验内容 用matlab 编程得到Sor 法求线性方程组的算法为: function [x,n]=SOR(A,b,x0,w,eps,M) if nargin==4

eps= ; M = 200; elseif nargin<4 error return : elseif nargin ==5 M = 200; end if(w<=0 || w>=2) error; return; end D=diag(diag(A)); %求A的对角矩阵L=-tril(A,-1); %求A的下三角阵( U=-triu(A,1); %求A的上三角阵B=inv(D-L*w)*((1-w)*D+w*U); f=w*inv((D-L*w))*b; x=B*x0+f; n=1; %迭代次数 while norm(x-x0)>=eps x0=x; x =B*x0+f; n=n+1; if(n>=M) (

相关文档
相关文档 最新文档