计算方法实习报告
一．给出一个有效的算法和无效的算法计算积分
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
#include
void main()
{ double y0,y1;
y0=1/4.0*log(5.0);
cout<<"y0="< for(int n=1;n<=10;n++) { y1=-1.0/4.0*y0+1.0/(4.0*n); cout<<"y"< y0=y1; if(n%3==0) cout< } } 其结果为： y0=0.402359 y1=0.14941 y2=0.0876475 y3=0.0614215 y4=0.0471446 y5=0.0382138 y6=0.0321132 y7=0.027686 y8=0.0243285 y9=0.0216957 y10=0.0195761 Press any key to continue 2．无效算法 利用递推公式y(n-1)=-4y(n)+1/n, 又由广义积分中值定理可得y(n)=1/((4n+1)(4ζ+1)), ζ∈(0,1),则1/(5(n+1)) 程序为：#include #include void main() { float y9,y10; y10=3.0 /55.0; cout<<"y10="< for(int n=9;n>=0;n--) { y9=-4.0*y10+1.0/(n+1); cout<<"y"< y10=y9; if(n%3==0) cout< } } 其结果为：y10=0.0545455 y9=-0.118182 y8=0.583838 y7=-2.21035 y6=8.98427 y5=-35.7704 y4=143.282 y3=-572.877 y2=2291.84 y1=-9166.86 y0=36668.4 3．心得体会 由有效算法与无效算法的结果可以知道：有效算法的误差的传递是逐步缩小的，而无效算法的误差的传递是逐步扩大的。因此有效算法是数值稳定的算法。 二．用牛顿法求下列方程lnx+x-2=0的根。 1．算法 给定初值X0，ε为根的容许误差，η为|F（X）|的容许误差，N为迭代次数的容许值。 (1).如果F（X0）的导数等于零或者迭代次数大于N，则算法失败，结束；否则执行2。 (2).计算X1=X0-F（X0）/F’（X0）。 (3).若|X1-X0|<ε,或者|F（X1）|<η,则输出X1，程序结束，否则执行4。 (4).令X0=X1，转向1。 2.程序 #include #include #define EPS 1e-8 #define ETA 1e-8 #define N 100 float f(float x) { return log(x)+x-2; } float f1(float x) { return 1+1/x; } void main() { float x1,x0=1.5,d; cout<<"x0="< for(int i=0;i<=N;) { if(f1(x0)==0) break; else x1=x0-f(x0)/f1(x0); d=(x1-x0); if(fabs(f(x1)) { cout< break; } else x0=x1; i++; cout<<"x"< while(i%3==0) cout< } cout<<"The root of the equation is x="< } 3.结果 x0=1.5 x1=1.55672 x2=1.55715 1.55715 The root of the equation is x=1.55715 4．心得体会 通过牛顿迭代算法的编程实习，我了解到牛顿法把非线性方程线性化，加快了收敛速度。 三．编写一个用牛顿前插公式计算函数值的程序，要求先输出差分表，再计算x点的函数值，并应用于下面的问题： xi 20 21 22 23 24 yi 1.30103 1.32222 1.34242 1.36173 1.38021 求x=21.4时的插值多项式的值。 1．算法 (1).输入n=4,xi,yi(i=0,1,2,3,4). (2).计算各阶差分f00,f10,f20,f30,f40。 (3).计算函数值：f00+f10*t+f20*t*(t-1)/2+f30*t*(t-1)*(t-2)/6+f40*t*(t-1)*(t-2)*(t-3)/2 4; 2．程序 #include float NewTon(float f00,float f10,float f20,float f30,float f40,float t) { return f00+f10*t+f20*t*(t-1)/2+f30*t*(t-1)*(t-2)/6+f40*t*(t-1)*(t-2)*(t-3)/24; } void main() { float f1[4],f2[3],f3[2],f4; float f0[5]={1.30103,1.32222,1.34242,1.36173,1.38021}; float x[5]={20,21,22,23,24}; cout<<" "<<"差分表"< for(int i=0;i<=4;i++) { cout<<" x"< } cout< for(i=0;i<=4;i++) { cout<<" "< } cout< for(i=0;i<=3;i++) { f1[i]=f0[i+1]-f0[i]; cout< } cout< for(i=0;i<=2;i++) { f2[i]=f1[i+1]-f1[i]; cout< } cout< for(i=0;i<=1;i++) { f3[i]=f2[i+1]-f2[i]; cout< } cout< f4=f3[1]-f3[0]; cout< float t=21.4-x[0]; float N4=NewTon(f0[0],f0[1],f0[2],f0[3],f0[4],t); cout<<"N(21.4)="< } 3．结果 差分表 x0=20 x1=21 x2=22 x3=23 x4=24 1.30103 1.32222 1.34242 1.36173 1.38021 1阶差分0.0211899 0.0202 0.01931 0.0184801 2阶差分-0.000989914 -0.000890017 -0.000829935 3阶差分9.98974e-005 6.00815e-005 4阶差分-3.98159e-005 N(21.4)=3.48268 4．心得体会 由于此题目选用的节点是等距的，所以可用等距节点的插值公式代替节点的任意分布的牛顿插值公式，这样就可以简化插值公式，同时可以避免作除法运算，而且程序算法的设计也相对比较简单。 四．用改进欧拉方法（取h=0.1）求解下列微分方程初值问题： y’=(x^2+y^2),x∈[0,1.0] y(0)=0 1．算法 解一阶常微分方程初值问题y’=f(x,y),y(x0)=y0,x∈[a,b],将区间[a,b]做n 等分，取步长h=(b-a)/n. 欧拉公式为：y(i+1)=y(i)+hf(x(i),y(i)) 梯形公式为：y(i+1)=y(i)+h[f(x(i),y(i))+ f(x(i+1),y(i+1))]/2 改进欧拉法，采用公式：?(i+1)=y(i)+h f(x(i),y(i)) y(i+1)=y(i)+h[f(x(i),y(i))+ f(x(i+1),?(i+1))]/2 2．程序 #include #define N 10 float f(float x,float y) { return x*x+y*y; } void main() { float x0=0,y0=0,a=0,b=1.0; float h=(b-a)/N; cout<<"x[0]="< 曲线拟合的数值计算方法实验 【摘要】实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合（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）借款成本——借款筹资的资金成本主要表现在利息方面。 计算公式：K=I×（1-T）；式中：K—借款成本率，I—借款年利率，T—所得税率 例：某企业从银行取得一笔三年期的借款，年利率为15%，每年付息一次，到期还本付息。企业所得税率为33%，则借款的成本率为：K=15%×(1-33%)=10.05% （2）债券资金成本——债券资金成本不仅考虑债券利息，还要考虑发行债券的筹资费率。计算公式：；公式中：K—债券资金成本率；L—债券各年利息额；T—所得税率；Q—债券发行总额；f—债券筹资费率。 例：某企业发行5年七债券，面值总额为5000万元，平价发行，已知债券票面利率为10%，发行费率为3%，所得税率33%，问该债券的资金成本率为多少？ 债券各年利息额=5000×10%=500（万元） 则，该债券的资金成本率为： 主权资金成本的计算： （1）普通股资金成本计算，假定股利每年均按固定数额支付，计算公式： ；--------（1）式中：K—普通股资金成本率；D—普通股股利；Q—普通股股本总额；f—筹资费率。如果股利股息逐年提高，则：；G—增加股息率 （2）优先股的资金成本计算，优先股每年的股利是一定的，所以，计算公式同上（1）。例：某企业发行总面额为400万元的优先股，其年股利率为14%，筹资费率为3%，则， 优先股资金成本率为： 综合资金成本计算：综合资金成本率亦即加权平均资金成本率，它是以各种资金来源占企业全部资金的比重为权数，对个别资金成本和主权资金成本进行加权平均计算的结果， 计算公式：；式中：K—综合资金成本率；Wi—某项来源的资金额占全部资金额的比重；Ki —某项资金来源的资金成本。 例：某企业共有资金额10000万元，其中，通过银行长期借款1000万元，借款年利率10%；发行长期债券5000万元，债券年利率12%，筹资费用率为2%；发行普通股票4000万元，第一年股利率为10%，以后逐年递增3%，其筹资费率为3%，已知企业所得税率为33%，则该企业综合资金成本计算如下： 长期借款资金成本率：K=I×（1-T）=10%×（1-33%）=6.7% 长期债卷资金成本率： 普通股的资金成本率： 三种资金在全部资金中的权重：长期借款10%；长期债卷50%；普通股40% 所以，综合资金成本率： 资本成本的计算 （1）企业发行长期债券，面值800万元，发行价格为700万元，筹资费率2%，债券年利率为10%，所得税率为25%。计算该债券的资本成本。 （2）企业向银行借款200万元。筹资费率1%，补偿性余额为10%，借款年利率为8%，所得税率为25%。计算该借款的资本成本。 （3）公司发行优先股150万元，年股利率为12.5%，预计筹资费用为5万元。计算该优先股的资本成本。 研硕16《化工数值方法及Matlab应用》试题 班级姓名成绩 1．（15分）数值计算方法的主要研究对象有哪些？其常用基本算法主要包括哪三个方面？举例说明Matlab在解决化工数值计算问题方面有什么样实用价值？答：（1）数值计算方法的主要研究对象为非线性方程求根，插值法、曲线拟合、数值积分、常微分方程（组）、初值问题求解、线性和非线性方程组求解。（2）基本算法包括①离散化方法：用差商代替导数、差分代替微分等，将连续的数学问题转化为离散问题。②逼近方法：用简单函数的值近似代替求解困难或形式未知的复杂函数的值。③迭代法：用一个固定公式反复计算，对较为粗糙的根的近似值进行加工直到满足精度要求的方法。 （3）Matlab在解决化工数值计算问题的实用价值有：数值计算和符号计算功能；图形功能；MATLAB语言；功能性和学科性工具箱。 2．（10分）数值计算中的“曲线拟合”，一般有哪些方法？请至少指出四种，并简述各自的基本特点。 答：（1）拉格朗日插值：，优点在于不要求数据点事等间隔的，缺点是数据点不易过多，当数据比较多时，差值函数有偏离原函数的风险； （2）牛顿插值法：它不仅克服了“增加一个节点时整个计算工作必须重新开始”的缺点，而且可以节省乘、除法运算次数。同时，在牛顿插值多项式中用到的差分与差商等概念，又与数值计算的其他方面有着密切的关系。 （3）牛顿迭代法：牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根，此时线性收敛，但是可通过一些方法变成超线性收敛。 （4）区间二分法：优点：算法简单，容易理解，且总是收敛的。缺点：收敛速度太慢，浪费时间，二分法不能求复根跟偶数重根。 （5）最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。 3. （15分）在298K 下，化学反应 2OF 2=O 2+2F 2 的平衡常数为0.410 atm ，如在298K 下将OF 2 通入容器，当t=0 时为1 atm ，问最后总压是多少？取计算精度为10-3 。 解：首先写出求解问题的数学方程式。 假设气体是理想气体，由反应的化学计量式可知， 22222F O OF += 设氧的分压为p ，平衡时有p 21- p p 2。 平衡时，有()410.02142 3=-p p 整理得 0410.064.1640.1423=-+-p p p 函数关系为 ()0410.064.1640.1423=-+-=p p p p f 实验一误差分析 实验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 姓名：万轩 实验二插值法 第四章资本成本 ●学习目标与要求 1.理解资金成本的概念、内容、作用。资本成本是指企业筹措和使用资金需要付出的代价。广义上，企业不论筹集和使用短期的资金还是长期的资金，都要付出代价。狭义的资本成本仅指筹集和使用长期资金的成本。资本成本包括用资费用和筹资费用两部分。资本成本的意义是：（1）资本成本是比较筹资方式，进行资本结构决策的重要依据。（2）资本成本是评价投资效益，比较投资方案的主要标准。（3）资本成本是评价企业经营管理业绩的客观依据。 2.掌握个别资本成本的计算方法。个别资本成本是指各种长期资本的成本，包括债务资本成本和权益资本成本。债务资本成本如长期借款成本、债券成本等，权益资本成本有优先股成本、普通股成本、留存收益成本等。 3.掌握综合资本成本的计算方法。综合资本成本一般是以各种资本的比重为权数，对个别资金成本进行加权平均确定的。权数的选择主要有账面价值、现行市价、目标价值和修正账面价值。 4.理解边际资本成本的内涵，掌握其计算方法。资本的边际成本是指资本每增加一个单位而增加的成本，是企业追加筹资的成本。企业追加筹资时有时可能只采取某一种筹资方式。在筹资数额较大，或在目标资本结构既定的情况下，往往通过多种筹资方式的组合来实现。这时，边际资本成本需要按加权平均法来计算，而且其资本比例必须以市场价值确定。 ●学习重点与难点 学习重点：正确理解各种资本成本的概念，掌握资本成本的计算方法。 学习难点：资本成本是财务管理中的一个重要概念，其应用非常广，因此，正确应用资本成本是学习的难点。 ●练习题 一、单项选择题 1.不需要考虑发行成本影响的筹资方式是（）。 A.发行债券 B.发行普通股 C.发行优先股 D.留存收益 2.从资本成本的计算与应用价值看，资本成本属于（）。 A.实际成本 B.计划成本 C.沉没成本 D.机会成本 3.在个别资本成本的计算中,不必考虑筹资费用影响因素的是（）。 给大家分享下考研公共课的一些经验。 英语： 我的英语基础：大一考过四级，大二上学期考过六级。但是考过六级后学英语就少了，所以说我的英语还是比较弱的。在考研准备期间，我背了蛋核英语微信推送的文章，这些文章大多比较短小，句子结构也比较简单，容易理解记忆，可能会有同学说背这个有什么用么？我觉得虽然阅读不会出这样的文章，但是这本书对于写作和培养语感还是很重要的，或者说背这些文章会对你的作文能力产生潜移默化的影响。 其次所用到的参考书就是英语历年考试真题，市面上有很多这样的书籍，我当时用的是《木糖英语真题手译》，当然单词不能忘记，用《一本单词》即可，不过里面只包括近10年的真题，因此我还把自1985年以来的考研英语真题都复印了拿来做。之后听说1985-1994年的题都太老了，不太复合现在考研的逻辑了，所以这些年份的题都可以不做，但1995年后的题还是值得一做的，起码可以复习一下语法。资料都找全了，剩下的就是做题了。我复习英语就是一遍一遍的做真题，分析句子结构，句型，逐字逐句的翻译。就这样英语真题大概总共做了5、6遍吧。 其实考研英语是有个规律的，完形填空20个题，肯定是5个A，5个B，5个C，5个D，印象中这个规律从未打破，这是在木糖英语考研微信中学到的。我在考试的时候基本就是先凭能力做，然后根据这个规律再改答案，结果完型做的很不错。阅读理解基本也是这个规律，但是也有例外，有可能不是5555，而是5546，,4556等等，而且一般来说，一篇阅读五个题目，不会出现三个相同选项的，如果出现了，你可要仔细看看了. 政治： 由于没有对过答案，不知道分数的具体分布，望请见谅。对于曾经的“文科尖子生”，我从来不认为政治是个问题。结果证明它真的不是一个问题。从大纲出来开始买书复习，大纲看了一遍。练习题买了李凡的《政治新时器》，做了几章。没有做过别的练习题。考前做了各种各样的押题卷的选择题，这里做选择题，如果时间允许，多多益善，并以此查缺不露。真题本身可能不太重要，但它给你带来的考场上的愉悦和放松的心情对应考还是大有裨益的。大题的话也是看《政 本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 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" 心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。 ：第一节 1.下面那一种方法不是函数的表示方法？（） A．分析法 B．图示法 C．表格法 D．解析法 答题：A. B. C. D. （已提交） 参考答案：D 1.设，则x的定义域为？（） A． B． C． D． 答题：A. B. C. D. （已提交） 参考答案：C 2.下面那一句话是错误的？（） A．两个奇函数的和是奇函数 B．两个偶函数的和是偶函数 C．两个奇函数的积是奇函数 D．两个偶函数的积是偶函数 答题：A. B. C. D. （已提交） 参考答案：A 2.多选：可以看做是哪些基本初等函数的复合或有限次四则运算步骤组成？（） A． B． C． D． 答题：A. B. C. D. >>（已提交） 参考答案：ABCD 3.函数定义中包括哪两个要素？（） A．定义域 B．值域 C．对应法则 D．对称性 答题：A. B. C. D. >>（已提交） 参考答案：AC 4.函数与是相等的。（） 答题：对. 错. （已提交） 参考答案：× 5.函数与是相等的。（） 答题：对. 错. （已提交） 参考答案：× 第二节 1.某厂为了生产某种产品，需一次性投入10000元生产准备费，另外每生产一件产品需要支付3元，共生产了100件产品，则每一件产品的成本是？（） A．11元 B．12元 C．13元 D．14元 答题：A. B. C. D. （已提交） 参考答案：C 2.某产品每日的产量是件，产品的总成本是元，每一件的售价为元，则每天的利润为多少？（） A．元 B．元 C．元 D．元 答题：A. B. C. D. （已提交） 参考答案：A 第三节 1.的反函数是？（） A． B． C． D． 答题：A. B. C. D. （已提交） 参考答案：C 2.的反函数是？（） A． B． C． D． 答题：A. B. C. D. （已提交） 参考答案：B 3.下面关于函数哪种说法是正确的？（） A．它是多值、单调减函数 B．它是多值、单调增函数 C．它是单值、单调减函数 D．它是单值、单调增函数 实验报告 学院（系）名称： 主程序部分列选主元部分 实验结果： 一．列主元消去法 输入各个数据，最终使用列选主元法，得到结果为：x1=x2=x3=1二．高斯-赛德尔迭代法 输入各个数据，输出每一步迭代数据，最终结果为：x1=0.285716，附录（源程序及运行结果） 一．列主元高斯消去法 #include 资本成本的计算 借款的资本成本 （筹资费很小时可以略去不计） 筹资费率）（借款成本所得税税率）（年利率借款成本-1-1???= K =）（）（f -1L T -1i L ??=f -1）T -1（i =）T -1（i 债券的资本成本 筹资费率）（发行价所得税税率）（票面利率债券面值-1-1???= K =） （）（f -1P T -1i B ? 优先股资本成本 筹资费率）（发行价优先股每年股利-1?= K =） （f -1P D 普通股资本成本 股利固定增长率筹资费率） （股价第一年的股利 +?= -1K = g f 1P D 1+-）（ 留存收益资本成本 计算留存收益成本的方法主要有三种： （1）股利增长模型法 股利固定增长率股价 第一年的股利 += K =g P D 1+ （2）资本资产定价模型 无风险利率）（平均收益率贝塔系数无风险利率-?+=K =)R R R m f f β-?+（ （3）风险溢价 风险溢价债券成本+=K =C P B R K + 加权平均资本成本 债券的资金成本总资金债券资金留存收益的资金成本总资金留存收益普通股的资金成本 总资金 普通股资金 借款的资金成本总资金借款资金?+?+?+?= ωK = ∑=n 1 i j j K W 边际贡献总额·Q P V S M C ?-=-=）（b 息税前净利润·a -M EBIT = 净利润·I EBIT -=净利润 普通股每股收益·股数 N I EBIT EPS -= 杠杆原理的计算 杠杆：一个因素发生较小变动，导致其它因素发生较大变动。 经营杠杆系数·a -=M M DOL 财务杠杆系数·T D I EBIT EBIT DFL -- -= 1 复合杠杆系数·DFL DOL DTL ?= DOL. EBIT Q S ?→?a / DFL. EPS EBIT I ?→? DTL. EPS Q S I ?→?,a / 项目投资的计算 建设期资本化利息原始总投资投资总额+= 建设期资本化利息固定资产投资额固定资产原值+= 该年回收额 该年利息费用该年摊销额该年折旧该年净利润经营期净现金流量++++=【直线法计提折旧】 折旧年限 净残值 建设期资本化利息固定资产投资额折旧年限净残值固定资产原值折旧--+== 复习题(一) 一、填空题： 1、求方程0.5x2 101x 1 0的根，要求结果至少具有6位有效数字。已知 V10203 101.0099，贝卩两个根为x1 _____________________________ ， X2 ________________________________ .(要有计算过程和结果) 4 1 0 A A 1 4 1 2、0 1 4，则A的LU分解为。 1 2 A 3、 3 5，贝卩(A) ____________ ，A __________ . 4、已知f(1)「Q f(2)「2，f(3) =3，则用抛物线(辛卜生)公式计算求 3 得1 f(x)dx -------------------- ，用三点式求得f (1) ________________ . 5、f(1) 1，f(2) 2，f(3) 1，则过这三点的二次插值多项式中x2的系数 为_____ ，拉格朗日插值多项式为 _________________________ . 二、单项选择题： 1、Jacobi迭代法解方程组Ax b的必要条件是( ). A. A的各阶顺序主子式不为零 B. (A) 1 C a ii 0,i 1,2, ,n D|| A 1 2、设f(x) 3x99 5x 7，均差f[1,2,22, ,299]=(). D. 3 4、三点的高斯求积公式的代数精度为 （ ). A.3 B. -3 C. 5 D.0 2 2 3 A 0 5 1 3、设 0 0 7 ，则 (A )为( ). A. 2 B. 5 C. 7 分别用拉格朗日插值法和牛顿插值法求 f （x ）的三次插值多项式P 3（x ），并 求f （2）的近似值（保留四位小数）. 4、 取步长h 0.2,用预估-校正法解常微分方程初值问题 y 2x 3y y （0） 1 （0 x 1） 5、 已知 A. 2 B.5 C. 3 D. 4 5、幕法的收敛速度与特征值的分布 A.有关 B.不一定 C. 无关 三、计算题: 1、用高斯-塞德尔方法解方程组 4X ! 2X 2 X 3 11 X 1 4X 2 2X 3 18 2X ! X 2 5X 3 22 （°） /c c c\T ，取 x （°，°，°），迭 四次（要求按五位有效数字计算 ). 1 2、求A 、B 使求积公式 1 f （X ）dX A[f( 1) f (1)] 1 B [f (2)f （2）］ 的代数精 度尽量高，并求其代数精度；利用此公式求 I 21dx 1 x （保留四位小 数）。 3、已知 3. 什么是融资成本，它包含哪几个部分？简述其计算方法并举例说明 令狐采学 资金成本是指筹集资金所付出的代价。一般由两部分组成： （1）资金占用费：指使用资金过程中向资金提供者支付的费用，包括借款利息、债券利息、优先股利息、普通股红利及权益收益 （2）资金筹集费：指资金筹集过程中发生的各种费用，包括律师费、资信评估费、证券印刷费、发行手续费、担保费、承诺费、银团贷款管理费等 资金成本通常以资金成本率K来反映 融资成本计算方法： 1.借款资金成本计算 （1）静态计算方法 理想情况下 K d=i 现实社会 K d=i(1-t) 如果考虑银行借款筹资费率和第三方担保费率，则K d=（i+V d）（1-t）/(1-f) 其中i表示贷款年利率，t为所得税率，f为银行借款筹资 费率，V d为担保费率 V d=V/P n X100% V为担保费总额，P为贷款总额，n为担保年限 例题4-1 某企业为建设新项目需向银行贷款400万元，偿还期5年，年利率为10%，筹资费率为2%。同时，按照银行要求提供第三方担保，担保费为70万元，担保期5年。项目投产后应纳所得税税率为25%，请计算该项目贷款的资金成本。 解： （2）动态计算方法 把各种费用按一定的折现率折现到某一时间点，然后再按静态法计算 例4-2 项目期初向银行借款100万元，年利率6%，期限为3年，到期一次还清借款，筹资费率为5%，项目投资当年即生产即生产并盈利，所得税率为33%，试计算项目借款税后资金成本。 解：融资现金流量分析： 期初（0年）实际投资额：100-100*5%=95（万元） 在1,2,3年末，考虑抵税效果后，各年缴纳利息为： 100*6%*（1-33%）=4.02（万元） 在3年末需要偿还本金100万元。 现在我们按公式（4-1）计算资金成本率K，它是使实际 华南理工大学数学分析2011-2013考研解答 1. ($12'$) 求极限 $\dps{\lim_{n\to\infty}\sqrt{n}\sex{\sqrt[4]{n^2 +1}-\sqrt{n+1}}}$. 解答: $$\beex \bea \mbox{原极限} &=\lim_{x\to 0}\sqrt{\frac{1}{x}}\sex{\sqrt[4]{\frac{1}{x^2}-1}-\sqrt{\frac{1}{x}-1}}\\ &=\lim_{x\to 0} \frac{\sqrt[4]{1+x^2}-\sqrt{1+x}}{x}\\ &=\lim_{x\to 0} \sez{\frac{1}{4}(1+x^2)^{-\frac{3}{4}}\cdot \frac{1}{2}(1+x)^{-\frac{1}{2}}}\\ &=-\frac{1}{2}. \eea \eeex$$ 2. ($12'$) 确定函数项级数$\dps{\sum_{n=1}^\infty \frac{x^2}{n}}$ 的收敛域, 并求其和函数. 解答: 由$a_n=1/n$ 知收敛半径为$R=1$. 又$\dps{\sum_{n=1}^\infty \frac{x^2}{n}}$ 当 $x=-1$ 时收敛, 当 $x=1$ 时发散, 而收敛域为 $[-1,1)$. 另外, 在收敛域范围内, $$\bex \sum_{n=1}^\infty \frac{x^2}{n} =\sum_{n=1}^\infty\int_0^xt^{n-1}\rd t =\int_0^x \sum_{n=1}^\infty t^{n-1}\rd t =\int_0^x \frac{1}{1-t}\rd t=-\ln (1-x). \eex$$ 3. ($12'$) 设函数$f\in C^2(\bbR)$, 且$$\bex f(x+h)+f(x-h)-2f(x)\leq 0,\quad\forall\ x\in \bbR,\quad \forall\ h>0. \eex$$ 证明: 对 $\forall\ x\in\bbR$, 有 $f''(x)\leq0$. 证明: 由$$\bex 0\geq \lim_{h\to 0}\frac{f(x+h)+f(x-h)-2f(x)}{h^2} =\lim_{h\to 0}\frac{f'(x+h)-f'(x-h)}{2h}=f''(x) \eex$$ 即知结论. 4. ($12'$) 设$\beta>0$ 且$$\bex x_1=\frac{1}{2}\sex{2+\frac{\beta}{2}},\quad x_{n+1}=\frac{1}{2}\sex{x_n+\frac{\beta}{x_n}},\ n=1,2,3,\cdots. \eex$$ 试证数列 $\sed{x_n}$ 收敛, 并求其极限. 证明: (1) $$\bex x_n=\frac{1}{2}\sex{x_{n-1}+\frac{\beta}{x_{n-1}}} \geq \sqrt{\beta},\quad n=2,3,\cdots. \eex$$ (2) 设 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
资本成本练习题
2021华南理工大学基础数学考研真题经验参考书
太原理工大学数值计算方法实验报告
华南理工大学经济数学随堂练习答案
数值计算方法实验5
资金成本公式(整理版)
现代数值分析复习题
融资成本的定义与计算方法
华南理工大学数学分析-考研解答
(完整版)数值计算方法上机实习题答案
数值分析实验报告总结