佛山科学技术学院数值分析第二次上机报告
佛山科学技术学院
上机报告
课程名称数值分析
上机项目曲线拟合
专业班级08数学与应用数学（师范）姓名植伟良学号2008214133
一.上机目的
1.建立数学模型，用最少二乘法计算值
2.编写拟合程序求值，描绘出拟合直线或者曲线
3．熟悉运用编写程序
二.上机内容
三.上机方法与步骤
x=[0.10.20.40.50.60.70.80.91];
f=[1.9783.286.167.347.669.589.489.3011.2];
aa=polyfit(x,f,2);
y=polyval(aa,x);
plot(x,f,'r+',x,y,'k')
xlabel('x');
ylabel('y');
gtext('y=s1(x)')
练习：在最小二乘准则下对下面一组数据作二次多项式拟合
四.上机结果
运用程序编写结算出结果如下，运用最少二乘法可以得到拟合曲线。
东南大学数值分析上机题答案
数值分析上机题 第一章 17.（上机题）舍入误差与有效数 设∑=-= N j N j S 2 2 11 ，其精确值为）111-23（21+-N N 。 （1）编制按从大到小的顺序1 -1 ···1-311-21222N S N +++=，计算N S 的通用 程序； （2）编制按从小到大的顺序1 21 ···1)1(111 222-++--+ -=N N S N ，计算N S 的通用程序； （3）按两种顺序分别计算210S ,410S ,610S ,并指出有效位数（编制程序时用单精度）； （4）通过本上机题，你明白了什么？ 解： 程序： （1）从大到小的顺序计算1 -1 ···1-311-21222N S N +++= ： function sn1=fromlarge(n) %从大到小计算sn1 format long ; sn1=single(0); for m=2:1:n sn1=sn1+1/(m^2-1); end end （2）从小到大计算1 21 ···1)1(111 2 22 -++--+-= N N S N function sn2=fromsmall(n) %从小到大计算sn2 format long ; sn2=single(0); for m=n:-1:2 sn2=sn2+1/(m^2-1); end end （3） 总的编程程序为： function p203()
clear all format long; n=input('please enter a number as the n:') sn=1/2*(3/2-1/n-1/(n+1));%精确值为sn fprintf('精确值为%f\n',sn); sn1=fromlarge(n); fprintf('从大到小计算的值为%f\n',sn1); sn2=fromsmall(n); fprintf('从小到大计算的值为%f\n',sn2); function sn1=fromlarge(n) %从大到小计算sn1 format long; sn1=single(0); for m=2:1:n sn1=sn1+1/(m^2-1); end end function sn2=fromsmall(n) %从小到大计算sn2 format long; sn2=single(0); for m=n:-1:2 sn2=sn2+1/(m^2-1); end end end 运行结果：
佛山科学技术学院(佛山大学)-兽医学硕士生导师—马春全简介
佛山科学技术学院导师信息采录表姓 名 马春全贴相片（登记照）性 别 男 出生年月 1962年10月 院系名称 生命科学学院动物医学系 专业职称 兽医教授 最高学历/学位 硕士研究生、博士 导师类别 硕导 研究方向 动物免疫病理学 导师说明 （主要简介个人重大事项，如二级教授岗位、国务院特殊津贴、全国优秀教师、省部级成果奖等） 导师简介一、著作 1. 《猪病混合感染鉴别诊断与防治彩色图谱》中国农业大学出版社2009.5 主编 2.《猪病学》（第三版）中国农业出版社2010.3主编 3.《猪病类症鉴别诊断与防治彩色图谱》中国农业出版社2011.5主编.二、论文 1. 雉鸡自然感染高致病性禽流感的病理组织学观察, 中国预防兽医学报, 2010,32(5): 379-382. 通讯作者 2. 双重PCR检测猪产肠毒素性大肠埃希菌耐热肠毒素基因, 黑龙江畜牧兽医2010,11:123-124. 通讯作者 3 .同位素检测卵黄抗体抑制猪大肠杆菌粘附试验 畜牧与兽医2011-5通讯作者 三、研究项目 1. 主持2010广东省生猪现代产业技术创新体系猪病防控岗位专家项目。经费30万元 2. 主持2011广东省生猪现代产业技术创新体系猪病防控岗位专家项目.经费30万元 3. 主持2012广东省生猪现代产业技术创新体系猪病防控岗位专家项目.经费30万元 4.主持广东高校畜禽疾病诊断防控技术开发中心”项目。经费70万元 四、科研奖励 1.主持《DataProc畜禽诊断及检测预警信息平台》 2010年1月经广东省科学技术厅组织专家鉴定为“国内同类研究领先水平”；2011年获得佛山市科学技术奖励二等奖 联系方式手机： E-mail：个人网页： 注：此表内容对外公布，业绩、获奖情况请填写有代表性项目，一般应为近五年（2007年1月之后）成果。
数值分析实验报告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 姓名：万轩 实验二插值法
佛山科学技术学院自命题考试科目考试大纲
2017年全国硕士研究生入学考试 佛山科学技术学院自命题考试科目考试大纲 （科目名称：结构力学科目代码：807） 一、考查目标 结构力学是佛山科学技术学院土木工程学科硕士学位研究生入学考试科目之一。该科目主要考查考生是否具备与土木工程有关的结构力学基本知识以及综合分析计算能力，以判别考生是否具备开展土木工程学科相关学术领域高水平、创新性科学研究的潜力。从而为国家培养具有较强分析问题和解决实际问题能力，并具有一定创新意识和创新能力的高层次专门技术人才。 该课程具体考查要求有： 1.掌握机动分析的原理及其应用； 2.掌握静定结构内力和位移的计算方法； 3.掌握超静定结构内力和位移的各种计算方法； 4.掌握静定结构影响线的绘制方法及其应用； 5.掌握结构动力学理论与计算。 二、考试形式与试卷结构 （一）试卷成绩及考试时间 本试卷满分为150分，考试时间180分钟。 （二）答题方式 答题方式为闭卷、笔试。 （三）试卷内容结构 各部分内容所占分值为： 1.机动分析(5~10分) 2.静定结构的受力分析和位移计算(20~28分) 3.力法(20~28分)
4.位移法(20~28分) 5.力矩分配法(5~10分) 6.影响线及其应用(15~20分) 7.矩阵位移法(20~25分) 8.结构动力学(20~28分) （四）试卷题型结构 1.是非题：5小题，共15分； 2.选择题：5小题，共20分； 3.填空题：5小题,共25分； 4.计算题：6小题,共90分。 （五）主要参考书目 李廉锟主编《结构力学》（第5版），高等教育出版社，2010年。 三、考查范围 1.机动分析 （1）自由度和联系； （2）平面几何不变体系的基本组成规则及其运用。 2.静定结构的受力分析 （1）静定梁、静定平面刚架的计算及内力图的绘制； （2）三铰拱的受力特点。 （3）计算静定平面桁架杆件内力的结点法和截面法； （4）组合结构的组成特点和内力计算； （5）静定结构受力特性。 3.静定结构的位移计算 （1）单位荷载法； （2）静定结构在荷载作用下的位移计算； （3）图乘法； （4）静定结构在非荷载因素（支座移动、温度变化）作用下的位移计算。
数值分析上机题目详解
第一章 一、题目 设∑ =-= N N j S 2 j 2 1 1，其精确值为)11 123(21+--N N 。 1) 编制按从大到小的顺序1 1 13112122 2-+??+-+-=N S N ，计算S N 的通用程序。 2) 编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ，计算S N 的通用程序。 3) 按两种顺序分别计算64210,10,10S S S ,并指出有效位数。（编制程序时用单精度） 4) 通过本次上机题，你明白了什么？ 二、通用程序 N=input('Please Input an N (N>1):'); AccurateValue=single((0-1/(N+1)-1/N+3/2)/2); Sn1=single(0); for a=2:N; Sn1=Sn1+1/(a^2-1); end Sn2=single(0); for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1); end fprintf('The value of Sn (N=%d)\n',N); fprintf('Accurate Calculation %f\n',AccurateValue); fprintf('Caculate from large to small %f\n',Sn1); fprintf('Caculate from small to large %f\n',Sn2); disp('____________________________________________________')
三、结果 从结果可以看出有效位数是6位。 感想：可以得出，算法对误差的传播有一定的影响，在计算时选一种好的算法可以使结果更为精确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象，容易产生较大的误差，求和运算从小数到大数所得到的结果才比较准确。
数值分析实验报告
数值分析实验报告 姓名：周茹 学号： 912113850115 专业：数学与应用数学 指导老师：李建良
线性方程组的数值实验 一、课题名字：求解双对角线性方程组 二、问题描述 考虑一种特殊的对角线元素不为零的双对角线性方程组（以n=7为例） ?????????? ?????? ? ???? ?d a d a d a d a d a d a d 766 55 44 3 32 211??????????????????????x x x x x x x 7654321=?????????? ? ???????????b b b b b b b 7654321 写出一般的n （奇数）阶方程组程序（不要用消元法，因为不用它可以十分方便的解出这个方程组） 。 三、摘要 本文提出解三对角矩阵的一种十分简便的方法——追赶法，该算法适用于任意三对角方程组的求解。 四、引言 对于一般给定的d Ax =，我们可以用高斯消去法求解。但是高斯消去法过程复杂繁琐。对于特殊的三对角矩阵，如果A 是不可约的弱对角占优矩阵，可以将A 分解为UL ，再运用追赶法求解。
五、计算公式（数学模型） 对于形如????? ?? ????? ??? ?---b a c b a c b a c b n n n n n 111 2 2 2 11... ... ...的三对角矩阵UL A =，容易验证U 、L 具有如下形式： ??????? ????? ??? ?=u a u a u a u n n U ...... 3 3 22 1 ， ?? ????? ? ?? ??????=1 (1) 1132 1l l l L 比较UL A =两边元素，可以得到 ? ?? ??-== = l a b u u c l b u i i i i i i 111 i=2, 3, ... ,n 考虑三对角线系数矩阵的线性方程组 f Ax = 这里()T n x x x x ... 2 1 = ,()T n f f f f ... 2 1 = 令y Lx =，则有 f Uy = 于是有 ()?????-== --u y a f y u f y i i i i i 1 1 11 1 * i=2， 3, ... ,n 再根据y Lx =可得到
佛山科学技术学院-期末总复习-学长整理-终极版-计算机组成原理
1.二进制代码10000000对应原码真值：-0。补码真值：-128。 反码真值：-127 2.计算机系统中采用补码运算的目的是为了简化计算机的设 计。 3.定点机 a)小数定点机：范围–(1 – 2-n) ~(1 – 2-n)。 b)整数定点机：范围–(2n - 1) ~(2n - 1)。 4.在定点运算器中，无论采用双符号位还是单符号位，必须有 溢出判断电路，它一般用异或门来实现。 5.在定点数运算中产生溢出的原因是______。 A. 运算过程中最高位产生了进位或借位 B. 参加运算的操作数超过了机器的表示范围 C. 运算的结果的操作数超过了机器的表示范围 D. 寄存器的位数太少，不得不舍弃最低有效位 6.下溢指的是______。 A. 运算结果的绝对值小于机器所能表示的最小绝对值 B. 运算的结果小于机器所能表示的最小负数 C. 运算的结果小于机器所能表示的最小正数 D. 运算结果的最低有效位产生的错误 7.存储单元是指________。 A. 存放一个二进制信息位的存储元 B. 存放一个机器字的所有存储元集合
C. 存放一个字节的所有存储元集合 D. 存放两个字节的所有存储元集合 8.扩展操作码是__________。 A. 操作码字段外辅助操作字段的代码 B. 操作码字段中用来进行指令分类的代码 C. 指令格式中的操作码 D. 一种指令优化技术，不同地址数指令可以具有不同的操作 码长度 9.异步控制方式常用于__________作为其主要控制方式。 A. 在单总线结构计算机中访问主存与外设时 B. 微型机的CPU控制中 C. 组合逻辑控制的CPU中 D. 微程序控制器中 10.在一个微周期中__________。 A. 只能执行一个微操作 B. 能执行多个微操作，但它们一定是并行操作的 C. 能顺序执行多个微操作 D. 只能执行相斥性的操作 11.在CPU中跟踪指令后继地址的寄存器是__________。B A. 主存地址寄存器 B. 程序计数器 C. 指令寄存器
2014级硕士研究生数值分析上机实习报告
2014级硕士研究生数值分析上机实习（第一次） 姓名：学号：学院： 实习题目：分别用二分法和Newton迭代法求方程x3■ 2x210x-20=0的根.实习目的：掌握两种解法，体会两种解法的收敛速度. 实习要求：用C程序语言编程上机进行计算，精确到8位有效数字. 报告内容： 1.确定实根的个数以及所在区间 2.将最后两次计算结果填入下表（保留8位数字）： 3.实习过程中遇到哪些问题？如何解决？有何心得体会?
4.两种解法的计算程序（此页写不下时可以加页）：
2014级硕士研究生数值分析上机实习（第二次）姓名：学号：学院： 实习题目：计算8阶三对角矩阵A=tridiag（0.235, 1.274, 0.235）的行列式.实习目的：掌握计算行列式的方法. 实习要求：首先选择一种算法，然后用C程序语言编程上机进行计算.报告内容： 1.简单描述所采用的算法: 2?计算结果： A 3.实习过程中遇到哪些问题？如何解决？有何心得体会?
4.写出C语言计算程序（此页写不下时可以加页）：
2014级硕士研究生数值分析上机实习（第三次） 姓名：学号：学院： 分别用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组实习题目： 2lx + 9.8y+ 3.4z= 6.7 <2.7x + 1.8y+ 7.2z= 2.4 8.6x + 1.5y + 3.4z = 1.9 实习目的：感受两种迭代法的收敛速度. 首先构造收敛的Jacobi迭代法和Gauss-Seidel迭代法，然后用实习要求： C程序语言编程上机进行求解，初始值均取为0,精确到4位小 数. 报告内容： 1.写出收敛的Jacobi迭代法和Gauss-Seidel迭代法:
数值计算实验报告
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
佛山科学技术学院2013---2014学年第一学期
佛山科学技术学院2013 ---20 14 学年第一学期 就业指导课程教学《职业体验报告 专业：网络工程班级1班姓名：郭雪城学号: 成绩: 《关于联诚地产业务员兼职的职业体验报告》 刖言 .体验时间.地点 2013年10月8日，季华三路智博E区 二.体验目的 1.通过本次实践使我能够把职业生涯规划课学到的知识融汇到实践上，更好的实现理论和实践的结合，为我以后工作和学习奠定初步的知识。 2.通过本次实践，使我能够亲身感受到有一个学生步入社会，和这个社会的发展同步成为职业人的过程。 3.通过本次实践对我完为人处世起到很大借鉴作用。 三.体验意义 纸上得来终觉浅，绝知此事要躬行。没有亲身的经历，就无法清楚地了解当今就业形势的真实状况，就只能让自己继续隔绝于社会之外。社会的复杂多样，变幻莫测，是在书本、网络和学校里无法感受和洞悉的到的。 正文 具体体验过程及描述 在这学期职业生涯规划课的最后一堂接到老师的通知后，我开始一脸的茫然，现在叫我去进行职业体验，以我现在的情况，又没有什么扎实的职业技能和工作经验，我该何去何从? 因为我家离学校较远吧，而且学校也就在佛山，所以就在佛山这里找了份家 教。本来打算在这里再找找有没有其他的工作可以当做职业体验的，到最后如果
没有的话，就把家教当做职业体验。因为我觉得家教实在算是比较简单的，不上职 也算业体验。整个学期过来，家教结束也没找到，发觉找工作也是比较不容易 的。面临即将到了交作业的时候，同学们陆陆续续的到学校了，2号晚上，旁边 寝室的一个同学跑过来问我们明天做不做兼职，是联诚地产的一个宣传活动，当时想想自己的职业体验总算也是有着落了。第二天6点多，我们怀着兴奋的心情做上公交车离开了学校，经过了一个半小时的车程到达了我们的目的地一一季华三路联诚地产智博E 区。由于负责人还没到，所以我们无所事事的等了会。 最后，一位年纪比我们大三四岁的小伙子走了出来，交代了一下我们今天 的任务，以及相关的注意事项。大概9点多我们坐上他的小面包车出发，前往工作地点。 10点多了，负责人和各个职员做好交涉，然后要求我们搭建好充气拱门，由于拱门太大，有些销售厅前面搭建不下，只好向外移动，基本到街道上了。发觉这工作真的很不容易，前前后后搭了6个充气拱门，时间都快到1点了。早饭还没有着落。搭建结束后，负责人叫我们自己找饭馆吃饭去，我们人生地不熟，找个饭馆又花了20几分钟。大概是10：50负责人给了我们一个锣，5个小喇叭, 5个智博E区的宣传板，便开始分配我们下午的下午的工作--------- 游街宣传。要求我们排成一列纵队，然后敲锣，对着喇叭喊：买房子，选智博E区！刚开始时真的很难喊出口，感觉自己一个大学生现在却像个街头小贩一样，特别是当街上的人群投来那异样的目光的时候，感觉特别的丢脸，所以一趟下来也没怎么喊。 第二趟时负责人一直在旁边督促着，要求我们扯开嗓子喊，想想豁出去了，反正别人又不认得我。突破了这层心理障碍，之后就开始一帮人继续开始宣传去了。 中午1点多，我们吃完饭回到智博E区休息，可怜我们只能在那坐一下，聊天。我和其中一个职员在闲聊，我说，你在房地产工作一段时间了吧？那个哥哥 说，他在智博E区工作有一年多了。“你的工作真好，我都羡慕死了，工作轻松, 工资又高。”“开玩笑，我工作好？我一个月才3000多而已，都不够我自己用，况且我还要补贴家里头。你今早也看到我被老板骂了一顿，我老是受他的气，哎! 谁叫我当年没听老师的教诲，没有好好地读书，现在想找份好工作都不行……” 听了那位哥哥的话，我开始思考自己的大学生活了，四年是那么的短啊，佛山科学技术学院--我的大学将从这里开始，我一定要好好珍惜大学生活。 那天的天气晴雨交加，各个身上衣服淋湿了又晒干了，到最后都很落魄。到了4点多，一群城管来了，要求我们拆了那个充气拱门，说是未经允许，妨碍了 交通。然后一帮人又收拾好充气拱门，负责人还被罚了1000多，然后他把火气发到我们头上，说我们没有看到城管通知他。发泄完之后，要求我们各自找车回去，到了最后回到学校将近9点，去吃了晚饭，拖着疲惫的身躯结束了第一天的职业体验。 虽然感觉很累人，但是之后几天几个人还是坚持下去了。我们收获了很多日常的生活中学不到的东西。
《数值分析》上机实验报告
数值分析上机实验报告 精选
《数值分析》上机实验报告 1.用Newton 法求方程 X 7-X 4+14=0 在（0.1，1.9）中的近似根（初始近似值取为区间端点，迭代6次或误差小于0.00001）。 1.1 理论依据： 设函数在有限区间[a ，b]上二阶导数存在，且满足条件 {}α?上的惟一解在区间平方收敛于方程所生的迭代序列 迭代过程由则对任意初始近似值达到的一个中使是其中上不变号 在区间],[0)(3,2,1,0,) (') ()(],,[x |))(),((|,|,)(||)(|.4;0)(.3],[)(.20 )()(.110......b a x f x k x f x f x x x Newton b a b f a f mir b a c x f a b c f x f b a x f b f x f k k k k k k ==- ==∈≤-≠>+K 令 )9.1()9.1(0)8(4233642)(0)16(71127)(0)9.1(,0)1.0(,1428)(3 2 2 5 333647>?''<-=-=''<-=-='<>+-=f f x x x x x f x x x x x f f f x x x f Θ 故以1.9为起点 ?? ?? ? ='- =+9.1)()(01x x f x f x x k k k k 如此一次一次的迭代，逼近x 的真实根。当前后两个的差<=ε时，就认为求出了近似的根。本程序用Newton 法求代数方程(最高次数不大于10)在（a,b ）区间的根。
1.2 C语言程序原代码： #include
数值分析实验报告1
实验一 误差分析 实验（病态问题） 实验目的：算法有“优”与“劣”之分，问题也有“好”与“坏”之别。对数值方法的研究而言，所谓坏问题就是问题本身对扰动敏感者，反之属于好问题。通过本实验可获得一个初步体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出：考虑一个高次的代数多项式 )1.1() ()20()2)(1()(20 1∏=-=---=k k x x x x x p 显然该多项式的全部根为1,2,…,20共计20个，且每个根都是单重的。现考虑该多项式的一个扰动 )2.1(0 )(19=+x x p ε 其中ε是一个非常小的数。这相当于是对（）中19x 的系数作一个小的扰动。我们希望比较（）和（）根的差别，从而分析方程（）的解对扰动的敏感性。 实验内容：为了实现方便，我们先介绍两个Matlab 函数：“roots ”和“poly ”。 roots(a)u = 其中若变量a 存储n+1维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,,,+n a a a ，则输出u 的各分量是多项式方程 01121=+++++-n n n n a x a x a x a 的全部根；而函数 poly(v)b =
的输出b 是一个n+1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“poly ”是两个互逆的运算函数。 ;000000001.0=ess );21,1(zeros ve = ;)2(ess ve = ))20:1((ve poly roots + 上述简单的Matlab 程序便得到（）的全部根，程序中的“ess ”即是（）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。 如果扰动项的系数ε很小，我们自然感觉（）和（）的解应当相差很小。计算中你有什么出乎意料的发现表明有些解关于如此的扰动敏感性如何 （2）将方程（）中的扰动项改成18x ε或其它形式，实验中又有怎样的现象 出现 （3）（选作部分）请从理论上分析产生这一问题的根源。注意我们可以将 方程（）写成展开的形式， ) 3.1(0 ),(1920=+-= x x x p αα 同时将方程的解x 看成是系数α的函数，考察方程的某个解关于α的扰动是否敏感，与研究它关于α的导数的大小有何关系为什么你发现了什么现象，哪些根关于α的变化更敏感 思考题一：（上述实验的改进） 在上述实验中我们会发现用roots 函数求解多项式方程的精度不高，为此你可以考虑用符号函数solve 来提高解的精确度，这需要用到将多项式转换为符号多项式的函数poly2sym,函数的具体使用方法可参考Matlab 的帮助。
数值分析实验报告模板
数值分析实验报告模板 篇一：数值分析实验报告(一)(完整) 数值分析实验报告 1 2 3 4 5 篇二：数值分析实验报告 实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。利用二分法求解给定非线性方程的根，在给定的范围内，假设f(x,y)在[a,b]上连续，f(a)xf(b) 直接影响迭代的次数甚至迭代的收敛与发散。即若x0 偏离所求根较远，Newton法可能发散的结论。并且本实验中还利用利用改进的Newton法求解同样的方程，且将结果与Newton法的结果比较分析。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。掌握二分法的原理，验证二分法，在选对有根区间的前提下，必是收
敛，但精度不够。熟悉Matlab语言编程，学习编程要点。体会Newton使用时的优点，和局部收敛性，而在初值选取不当时，会发散。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b) Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式xk?1?xk?f(xk) f'(xk) 产生逼近解x*的迭代数列{xk}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 xk?1?xk?rf(xk) 'f(xk) 其中r为要求的方程的根的重数，这就是改进的Newton 法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x);
数值分析上机题(matlab版)(东南大学)
数值分析上机题(matlab版)(东南大学)
数值分析上机报告
第一章 一、题目 精确值为)1 1 123(21+--N N 。 1) 编制按从大到小的顺序 1 1 131121222-+??+-+-= N S N ，计算S N 的通用程序。 2) 编制按从小到大的顺序 1 21 1)1(111222-+??+--+-= N N S N ，计算S N 的通用程序。 3) 按两种顺序分别计算6 42 10,10, 10S S S ,并指出有效位 数。（编制程序时用单精度） 4) 通过本次上机题，你明白了什么？ 二、通用程序 clear N=input('Please Input an N (N>1):'); AccurateValue=single((0-1/(N+1)-1/N+3/2)/2); Sn1=single(0); for a=2:N; Sn1=Sn1+1/(a^2-1); end Sn2=single(0); for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1); end fprintf('The value of Sn using different algorithms (N=%d)\n',N); disp('____________________________________________________') fprintf('Accurate Calculation %f\n',AccurateValue); fprintf('Caculate from large to small %f\n',Sn1); fprintf('Caculate from small to large %f\n',Sn2);
(完整版)佛山科学技术学院-期末总复习-学长整理-终极版-操作系统
什么是进程？什么是线程？进程与线程有何区别? 答： （1）进程是具有独立功能程序在某个数据集合上的一次执行过程。（2分） （2 ）线程是进程内的一个执行实体或执行单元。（2分） （3）进程和线程的区别： （a）不同进程的地址空间是独立的，而同一进程内的线程共享同一地址空间。一个进 程的线程在另一个进程内是不可见的。 （b）在引入线程的操作系统中，进程是资源分配和调度的单位，线程是处理机调度和分配的单位，资源是分配给进程的，线程只拥有很少资源，因而切换代价比进程切换低。（2 分） 1）请叙述分段式存储管理的基本原理，通过绘制地址变换原理图说明分段式存储管理中逻 辑地址到物理地址的变换过程。（10分） 2）请叙述SPOOLING系统的概念和组成。（10分） 1）请叙述段式存储管理的基本原理，结合动态段式管理说明地址变换过程。 答：段式存储管理的基本思想是：把程序按照内容或过程（函数）关系分成段，每个段由自 己的名字。一个用户作业或进程所包含的段对应一个二维线性虚拟空间，即一个二维虚拟存 储器。段式存储管理程序以段为单位分配内存，然后通过地址映射机构把段式虚拟存储器地 址转换成实际的内存物理地址。段式存储管理把那些机场访问的段驻留内存，而把那些不经 常访问的段存放到外存，需要的时候再调到内存。 段式管理技术通过段表实现内存分配和回收工作。下面以进程某条指令访问数据为例说明地 址变换过程。 当进程的某条指令发出数据访问指令，系统根据指令中给定的虚拟地址，虚拟地址中包括段号和段内偏移，首先通过段表地址寄存器查找该段号的内存物理起始地址，如果该段不在内 存则将该段调入内存，然后将得到的段的内存起始地址加上指令中的页内偏移，至此得到该 指令访问数据的物理地址。参见下图。 2）请叙述SPOOLING系统的概念和组成。 答：为了缓和CPU的高速与I/O设备的低速之间的矛盾，引入了脱机输入、脱机输出，使外 围设备操作与CPU数据处理同时进行，实现了联机情况下的同时外围设备操作的技术称为SPOOLING其组成主要有三部分：1）输入输出井：在磁盘上开辟的用于缓存进程数据的输入、输出存储空间；2）输入、输出缓冲区：在内存开辟的用于缓和CPU与磁盘速度不匹配的矛 盾；3）输入、输出处理进程：用于处理输入、输出过程。组成如下图：
(完整版)哈工大-数值分析上机实验报告
实验报告一 题目：非线性方程求解 摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。 前言：（目的和意义） 掌握二分法与Newton法的基本原理和应用。 数学原理： 对于一个非线性方程的数值解法很多。在此介绍两种最常见的方法：二分法和Newton法。 对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。重复运行计算，直至满足精度为止。这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式 产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。另外，若将该迭代公式改进为 其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。 程序设计： 本实验采用Matlab的M文件编写。其中待求解的方程写成function的方式，如下 function y=f(x); y=-x*x-sin(x); 写成如上形式即可，下面给出主程序。 二分法源程序： clear %%%给定求解区间 b=1.5; a=0;
%%%误差 R=1; k=0;%迭代次数初值 while (R>5e-6) ; c=(a+b)/2; if f12(a)*f12(c)>0; a=c; else b=c; end R=b-a;%求出误差 k=k+1; end x=c%给出解 Newton法及改进的Newton法源程序：clear %%%% 输入函数 f=input('请输入需要求解函数>>','s') %%%求解f(x)的导数 df=diff(f);
佛山科学技术学院校级科研成果奖励及评选办法
佛科院科…2004? 8号 佛山科学技术学院 校级科研成果奖励及评选办法 （试行） 第一条：为了促进我校教师及专业技术人员开展科研活动，促进我校科研水平和效益的提高，为地方经济建设和社会发展多作贡献，参照有关文件精神，结合我校的实际，特制定我校科研成果奖励及评选办法。 第二条：奖励类别、等级、名额、时限 1.类别：校级科研成果奖分为优秀学术论著、优秀学术论文、优秀自然科学应用成果、优秀社会科学应用成果四个类别。 2.等级和名额：各类别校级科研成果奖均设一、二、三等奖。分为自然科学优秀论著（50～60篇）、社会科学优秀论著（40～50篇）、优秀自然科学应用成果（12～18项）、优秀社会科学应用成果（12～15项）等四类。实行精神奖励与物质奖励相结合，以精神奖励为主的原则，对获奖项目授予奖状、证书和奖金。 3.时限：校级科研成果奖由原每两年评选一次调整为每三年评选一次（一般在学校每三年一次的科研大会上奖励）。 第三条：奖励范围和申报条件
1．奖励范围：项目第一完成单位（或第一完成人）必须为我校（教职工），论文、著作必须署名我校；成果必须是调入我校后完成的，并经科研主管部门立项，项目经费进入学校财务管理（自然科学类入账经费五万元以上，人文社科类五千元以上，校级课题经费不受此限制）；教材和教研成果不列入科研成果奖的评选范围。 2．申报条件 所有参评成果，必须是遵循坚持为社会主义服务、为地方经济建设和科技发展服务的方针，具有较好的科学性、创新性，并符合下列条件之一的成果： （一）科研成果通过科研主管部门鉴定（或结题）、自然科学类取得“国内先进”以上评价。 （二）科研成果取得较大的经济效益，并具有广阔的应用前景（有相关财务部门证明）。 （三）科研成果应用取得较大的社会效益（具有县级以上党政主管部门的证明）。 （四）在国内外权威以上刊物上发表的论文，或被SCI、EI、SSCI、ISTP、ISSHP索引，或新华文摘摘录、人大复印资料全文转载。 （五）由国内外权威出版社出版的学术专著，达到一定发行量，并多次被其他作者引用（肯定性地引用或作为参考文献）。 （六）已技术转让的专利，其转让金已按规定上交学校。 （七）有较大社会影响的艺术作品，并入选省级以上创作竞赛成果。 第四条：等级评定方式及条件 学校届时组织专家评审小组对申报奖项进行评审。 专家评审小组将严格遵循全面地、公正地评价申报项目的科学性、学术性、先进性、创新性及其经济效益和社会效益等的原则评审各审报奖项，各等级奖项的基本条件如下：
数值分析2016上机实验报告
序言 数值分析是计算数学的范畴，有时也称它为计算数学、计算方法、数值方法等，其研究对象是各种数学问题的数值方法的设计、分析及其有关的数学理论和具体实现的一门学科，它是一个数学分支。是科学与工程计算（科学计算）的理论支持。许多科学与工程实际问题（核武器的研制、导弹的发射、气象预报）的解决都离不开科学计算。目前，试验、理论、计算已成为人类进行科学活动的三大方法。 数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。现在面向数值分析问题的计算机软件有：C,C++,MATLAB,Python,Fortran等。 MATLAB是matrix laboratory的英文缩写，它是由美国Mathwork公司于1967年推出的适合用于不同规格计算机和各种操纵系统的数学软件包，现已发展成为一种功能强大的计算机语言，特别适合用于科学和工程计算。目前，MATLAB应用非常广泛，主要用于算法开发、数据可视化、数值计算和数据分析等，除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。 本实验报告使用了MATLAB软件。对不动点迭代，函数逼近（lagrange插值，三次样条插值，最小二乘拟合），追赶法求解矩阵的解，4RungeKutta方法求解，欧拉法及改进欧拉法等算法做了简单的计算模拟实践。并比较了各种算法的优劣性，得到了对数值分析这们学科良好的理解，对以后的科研数值分析能力有了极大的提高。
目录 序言 (1) 问题一非线性方程数值解法 (3) 1.1 计算题目 (3) 1.2 迭代法分析 (3) 1.3计算结果分析及结论 (4) 问题二追赶法解三对角矩阵 (5) 2.1 问题 (5) 2.2 问题分析（追赶法） (6) 2.3 计算结果 (7) 问题三函数拟合 (7) 3.1 计算题目 (7) 3.2 题目分析 (7) 3.3 结果比较 (12) 问题四欧拉法解微分方程 (14) 4.1 计算题目 (14) 4.2.1 方程的准确解 (14) 4.2.2 Euler方法求解 (14) 4.2.3改进欧拉方法 (16) 问题五四阶龙格-库塔计算常微分方程初值问题 (17) 5.1 计算题目 (17) 5.2 四阶龙格-库塔方法分析 (18) 5.3 程序流程图 (18) 5.4 标准四阶Runge-Kutta法Matlab实现 (19) 5.5 计算结果及比较 (20) 问题六舍入误差观察 (22) 6.1 计算题目 (22) 6.2 计算结果 (22) 6.3 结论 (23) 7 总结 (24) 附录
数值分析实验报告
学生实验报告实验课程名称 开课实验室 学院年级专业班 学生姓名学号 开课时间至学年学期
if(A(m,k)~=0) if(m~=k) A([k m],:)=A([m k],:); %换行 end A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); %消去end end x=zeros(length(b),1); %回代求解 x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end y=x; format short;%设置为默认格式显示,显示5位 （2）建立MATLAB界面 利用MA TLAB的GUI建立如下界面求解线性方程组： 详见程序。 五、计算实例、数据、结果、分析 下面我们对以上的结果进行测试，求解：
? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - - - 7 2 5 10 13 9 14 4 4 3 2 1 13 12 4 3 3 10 2 4 3 2 1 x x x x 输入数据后点击和，得到如下结果： 更改以上数据进行测试，求解如下方程组： 1 2 3 4 43211 34321 23431 12341 x x x x ?? ???? ?? ???? ?? ???? = ?? ???? - ?? ???? - ???? ?? 得到如下结果：