科学计算与编程程序.doc

１.求解非线性方程的所有解：fzero； 题：求解y=x.*sin(1./x)-0.2*exp(-x)在[0,2]的解。 1.写出函数： function y=myfun01(x) y=x.*sin(1./x)-0.2*exp(-x); end 2.再新建一个文件： clc close all clear all figure fplot(@myfun01,[0,2]) grid on （画图） x1=fzero(@myfun01,[0.2,0.4]) (找出零点所在的区间) myfun01([x1]) （验算是否是方程的解） ２.使用弦截法求解非线性方程的解； 题3-1为例 步骤： 1.用function写出求解函数； function y=myfunction(x) y=0.5*exp(x/3)-sin(x); end 2.用function创建弦截法求解非线性方程的程序secant函数： secant函数如下： function r=secant(f,a,b,n) fa=f(a); fb=f(b); for k=1:n m=a-(b-a)*fa/(fb-fa); fm=f(m); if(fm==0) break elseif(fa*fm0) b=m;fb=fm; else a=m;fa=fm; end end r=a-(b-a)*fa/(fb-fa); 3.在主文件下（新建一个文件做主文件）画图确定方程解的个数及解所在区间 clc close all clear all f=@myfunction; figure(1) fplot(f,[0,10]) grid on x1=secant(f,0,1,100) x2=secant(f,1,2,100) ３.使用fsolve求解非线性方程组的解； P75习题3-4 程序如下： 1、先建一个myfun文件 function F=myfun01(X) x=X(01); y=X(02); F(1)=x.*exp(x*y+0.8)+exp(y^2)-3; F(2)=x^2-y^2-0.5*exp(x*y); 2、在建Edit文件编写主程序 clc clear all close all [fsolve(@myfun01,[1 1]) fsolve(@myfun01,[1 -1]) fsolve(@myfun01,[-1 1])] 求出三个区间的三个解 ４.求解线性方程组的解（恰定方程组、超定方程组）； 第四题：（4-1到4-6都可以用这个方法写出来） 本题为4-4： clc clear all close all A=[1,1,2;-1,2,1;1,-1,2]; b=[4,4,-1]; x=pinv(A)*b ５.使用多项式进行曲线拟合； P101例1 拟合一次多项式： clc clear all close all x=[1,1.2,1.4,1.6,1.8]; y=[1,1.0954,1.1832,1.2649,1.3416]; plot(x, y,’r.’,’markersize’, 25) grid on 先画图，再接着写下面的 x=x’; A=[X.^1;X.^0]’ ; b=y’ ; a=pin(A)*b 或a=A\b hold on 写到这运行程序，第一个数是x的系数，第二个是常数项 fplot(‘0.4264*x+0.5801’,[0.8,2]) 图线区间 拟合二次多项式： 第5.4题 1、程序： clc close all clear all x=[0.5 1.0 1.5 2.0 2.5 3.0]; y=[1.75 2.45 3.81 4.80 8.00 8.60]; figure plot(x,y,r.,markersize,25) grid on A=[x.^2;x.^1;x.^0]; b=y; a=A\b hold on 先写到这里，运行，得到三个数，第一个是二次项系数，以此类推，写出下式。 fplot(0.4900*x^2+1.2501*x+0.8560,[0.5 3]) 运行得到连线图。 ６.使用经验公式进行曲线拟合： ａ幂函数：y=xα 以题5-5为例 程序： clc close all clear all x=[2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0]; y=[40.1 54.3 73.1 98.8 133.4 180.2 243.0 328.0 442.9 597.7 806.8]; figure(2) plot(x,y,.r,markersize,30) grid on hold on X=log