1、1线性方程组数值解法与非线性方程求解淮 海 工 学 院实 验 报 告 书课程名称: 数学实验 实验名称:线性方程组的数值解法与非线性方程求解班 级 数学091 姓 名: 耿萍 学号: 090911107 日 期: 2018.4.27 地点数学实验室指导教师: 曹卫平 成绩:数 理 科 学 系1. 实验目的:(1) 掌握线性方程组的常用数值解法,包括高斯消去法、LU分解法以及校正法。(2) 体验数值计算的时间复杂度和计算规模的关系。(3) 加深对数值计算误差的理解。(4) 学习使用迭代法等算法,求解非线性方程。(5) 学习如何使用MATLAB解非线性方程组和方程组。2. 实验内容:、列出求各负载
2、电流I1,I2,In的方程;2)设R1=R2=Rn=R,r1=r2=rn=r,在r=1,R=6,V=18,n=10的情况求I1,I2,In及总电流I0。 2)种群的的繁殖与稳定收获:种群的数量因素因繁殖而增加,因自然死亡而减少,对于人工饲养的种群比如家畜)而言,为了保证稳定的收获,各个年龄的种群数量应维持不变。种群因雌性个体的繁殖而改变,为方便起见一下种群数量均指其中的雌性。种群年龄记作bk每个雌性个体一年繁殖的数量),自然存活率记作sk=1-dk,dk为一年的死亡率),收获量记作hk,则来年年龄k的种群数量xk应为x1=cigmabkxk,xk+1=skxk-hkk=1,2,3,n-1)。要
3、求各个年龄的种群数量每年维持不变就是要使xk=xk若bk,sk已知,给定收获量hk,建立求各个年龄的稳定种群数量xk的模型用矩阵、向量表示)2)设n=5,b1=b2=b5=0,b3=5,b4=3,s1=s4=0.4,s2=s3=0.6,如要求h1h5为500,400,200,100,100,求x1x5.3)要使h1h5均为500,如何达到? 3) 1)小张夫妇以按揭的方式贷款买了1套价值为20万的房子,首付了5万元。每月还款1000元,15年还清。问贷款利率是多少?2)某人欲贷款50万元购房,他咨询了两家银行,第一家银行开出的条件是每月还4500元,15年还清;第二家银行开出的条件是每年450
4、00元,20年还清。从利率方面看,哪家银行较优惠简单地假设年利率=月利率*12) 计算序列ykk=0,1,2,),其中b取1.3,2.8,3.2,3.5,3.55,3.7,任意取y00y0,1),观察其收敛性。3. 实验步骤:记r1rn上的电流为iin。设电源负极为电势为0,电阻R1上对应节点电压为V1,对于任意节点,根据KCL定律列出方程:而,可得:k=2,3,n-1;k=1时,为与上式形式一致,化为k=m n=10。 %由题目要求设定A11=sparse(1:n-1,1:n-1,-1,n,n。 %定义A的对角元素,除A12=sparse(n,n,-0.5,n,n。 %定义(n,nA1=A1
5、1+A12。 %对角元素A2=sparse(1:n-1,2:n,0.5,n,n。 %输入A的上次对角元素A3=sparse(2:n,1:n-1,0.5,n,n。 %输入A的下次对角元素A=A1+A2+A3。b1=0.5*ones(n,1。 %b的除第一项元素b2=sparse(1,1,18,n,1。 %b的第一项元素b=b1-b2。R=Ab得到结果:R = 26.0000 17.0000 9.0000 2.0000 -4.0000 -9.0000 -13.0000 -16.0000 -18.0000 -19.0000所以各阻值为(R1,R2,R10=(26,17,9,2,-4,-9,-13,-
6、16,-18,-19总电阻R0即输入等效电阻)为,又得到要使各年龄种群数量每年维持不变即,依题意得用矩阵形式表示原方程组为:,2代入题中数据,3要使h1h5均为500,则h变为:输入程序: format bank。A1=0.4,-1,0,0,0 0,0.6,-1,0,0 0,0,0.6,-1,0 0,0,0,0.4,-1 -1,0,5,3,0。h1=500,400,200,100,0。 x1=A1h1 format bank。A2=0.4,-1,0,0,0 0,0.6,-1,0,0 0,0,0.6,-1,0 0,0,0,0.4,-1 -1,0,5,3,0。h2=500,500,500,500,
7、0。x2=A2h2format bank。A3=0.6,-1,0,0,0 0,0.8,-1,0,0 0,0,0.8,-1,0 0,0,0,0.6,-1 -1,0,1,2,0。h3=500,500,500,500,0。x3=A3h3 结果:x1 = 8481.01 2892.41 1335.44 601.27 140.51x2 = 10981.01 3892.41 1835.44 601.27 -259.49x3 = 13467.74 7580.65 5564.52 3951.61 1870.97从x1可以看出,第5年龄段:x5=140.5100=h5 ,说明收获量h5可以达到100。从x2可以
8、看出,x5为-259.49,但种群数量不可能为负数,在本题所给条件下,无法使h1h5均为500。从x3可以看出,x5=1870500=h5,说明收获量h5可达到500,从而h1h5均可达到500。由题目已知条件,假设第i月月初待还贷款为,贷款月利率为r,则可列出:=150000=*(1+r-1000=1000/r+(-1000/r2 记第一家银行月利率为s,第二家银行年利率为t,则:=4500/s+(-4500/r输入程序:1)r=fzero(inline(1000/r+(150000-1000/r*(1+r180,12)r1=fzero(inline(4500/s+(500000-4500/
9、s*(1+s180,1r2=fzero(inline(45000/t+(500000-45000/t*(1+t20,1if r1。else disp(第二家月利率小。end实验结果:r =0.0021r1 =0.0059r2 =0.0639第二家月利率小=0.5for k=1:(n-1 y(k+1=1.3*y(k*(1-y(k。end。y。y = Columns 1 through 12 0.5000 0.3250 0.2852 0.2650 0.2532 0.2458 0.2410 0.2378 0.2356 0.2341 0.2331 0.2324 Columns 13 through 20 0.2319 0.2316 0.2313 0.2312 0.2310 0.2310 0.2309 0.23094. 实验数据记录及分析或程序及运行结果):结果已在上一步中给出)通过这次实验,是我对线性方程和非线性方程有了更进一步的认识,并且对于它们的解法也有了更深入的了解。它们可以解决很多实际问题。所以熟练的掌握它们对于解决现实中的问题有很大的帮助。申明:所有资料为本人收集整理,仅限个人学习使用,勿做商业用途。
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2