matlab练习.doc

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 Ds1421 学 院 电气 年级 2012级电气 专业班 3 学 生 姓 名 黄桂川 学 号 开 课 时 间 2013 至 2014 学年第 1 学期 总 成 绩 教师签名 数 学 与 统 计 学 院 制 开课学院、实验室： 实验时间 ： 年 月 日 课程 名称 数学实验 实验项目 名 称 方程求解 实验项目类型 验证 演示 综合 设计 其他 指导 教师 龚劬 成 绩 实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习3．针对实际问题，试建立数学模型，并求解。function f=y(x) f=x^2-x; 第二个： function g=findroot(x) m=x(1);n=x(2); while abs(n-m)1e-10 if (y((m+n)/2)*y(m))0 n=(m+n)/2; k=(m+n)/2; else m=(m+n)/2; k=(m+n)/2; end end 窗口输入：x=[0.5 2];findroot(x) 答案： ans = 1.0000 总结：在输入x=[0.5 2]时，要保证在0.5到2之间有零点，否则无法满足abs(y(n-m))1e-10从而导致电脑死机 5. 设非线性方程组为 其中已知，随机产生数据后，用fsolve解这个方程组。 6. 已知平面上的五条直线（aix + biy + ci = 0, i = 1,...,5），如何判断它们是否交于一点。 求出两直线的交点带入其他直线中验证 function l=fa(v) syms x y a1 a2 a3 a4 a5 a6; [x,y]=solve(a1*x+y*a2+a3=0,a4*x+y*a5+a6=0,x,y); x0=subs(x,[a1 a2 a3 a4 a5 a6],[A(1,1) A(1,2) A(1,3) A(2,1) A(2,2) A(2,3)]) y0=subs(y,[a1 a2 a3 a4 a5 a6],[A(1,1) A(1,2) A(1,3) A(2,1) A(2,2) A(2,3)]) f3=A(3,1)*x0+A(3,2)*y0+A(3,3); f4=A(4,1)*x0+A(4,2)*y0+A(4,3); f5=A(5,1)*x0+A(5,2)*y0+A(5,3); if f3~=0|f4~=0|f5~=0 l=不交于一点; else l=交于一点; end end 在窗口输入： A=[1 2 1;1 3 1;1 4 1;1 5 1;1 1 1];fa(A) 答案：x0 = -1 y0 = 0 ans = 交于一点 7. 使用fsolve计算方程组的解时，为验证初值是否对解有影响，采用随机产生的100组随机数作为初始值，依次进行求解。 建立一个名为f的M文件，用于存放输入函数组 function Fangcheng=f(x) Fangcheng=[x(1)+x(2)^2-13;log(2*x(1)+x(2))-x(1)^x(2)+2]; 主函数 A=rand(100,2)*5; for n=1:100 a=[A(n,1) A(n,2)] f1=fsolve(f,a); x1=f1(1) x2=f1(2) end 答案： 第几次 a x1 x2 1 2.4400 3.4123 1.4880 3.3929 2 2.4892 1.7906 1.4880 3.3929 3 4.6799 4.9346 1.4880 3.3929 5 1.9464 0.4199 1.4880 3.3929 6 0.5857 1.2515 1.4880 3.3929 7 1.2021 4.0567 1.4880 3.3929 8 3.4245 0.4221 12.5709 0.6551 9