1、一、 科学计算的算法及其举例应用和利用 MATLAB 自带函数实现科学计算法：弦截法非线性方程实例: 采用弦截法求方程 lgx+ 2 在区间1，4上的一个根。x1.算法说明弦截法的算法过程如下：（1） 过两点(a,f(a),(b,f(b) 作一直线，它与 x 轴有一个交点，记为 x1（2） 如果 f(a)f(x1)0)disp(两端点函数值乘积大于0！);return;elsetol=1;fa=subs(sym(f),findsym(sym(f),a);fb=subs(sym(f),findsym(sym(f),b);root=a-(b-a)*fa/(fb-fa);while(toleps)调
2、用 Secant 进行运算输出结果 r =1.8773结束输入条件：sqrt(x)+log(x)-2积分区间：（1,4）开始r1=root;fx=subs(sym(f),findsym(sym(f),r1);s=fx*fa;if(s=0)root=r1;elseif(s0)root=b-(r1-b)*fb/(fx-fb);elseroot=a-(r1-a)*fa/(fx-fa);endendtol=abs(root-r1);endend实例源代码：r=Secant(sqrt(x)+log(x)-2,1,4)运行结果：r =1.8773二、 分析单自由度阻尼系统的阻尼系数对其固有振动模态的影响1
3、 流程图2.源程序代码：结束绘制 x,t 二维图象subplot(1,2,1)j=1j10?plot(t,x(j,:)j=j+2开始输入 wn， x0， v0 ，tfj=1zeta=0.1*jwd=wn 12A= /wd(0+0)2+(0)2)a=atan2(wd*x0,v0+zeta*wn*x0)j 10?t=0j=j+1ttf ?t=t+tf/1000x(j,:)=A* * sin(wd*t+a)subplot(1,2,2),mesh(x)NoYesYesNoNoYes3.运行结果：三Appollo 卫星的运行轨迹1.流程图2.源程序代码：函数文件：开始输入初值t,y=ode45(appollo,0,20,x0,options)绘制图形 plot(y(:,1),y(:,3)结束命令文件：3.运行结果
