第04章 函数与方程.docVIP

函数和方程 1、求多项式的根，并分析误差大小 （1） 解： A=[1,1,1];roots(A) ans = -0.5000 + 0.8660i -0.5000 - 0.8660i y=polyval(A,ans) ％验证 y = 1.0e-015 * 0.3331 0.3331 （2） 解： A=[3 0 -4 0 2 -1];x=roots(A) x = -0.9479 + 0.3845i -0.9479 - 0.3845i 1.0000 0.4479 + 0.3435i 0.4479 - 0.3435i y=polyval(A,x) ％验证 y = 1.0e-013 * -0.0144 - 0.1138i -0.0144 + 0.1138i -0.0888 0.0044 - 0.0083i 0.0044 + 0.0083i （3） 解： A=[5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0];x=roots(A) y=polyval(A,x) (4) 解：p1=[2 3];p2=conv(p1,p1); p=conv(p2,p1); %(2x+3)3 n=length(p); for (i=1:n-1) q(i)=p(i); end q(n)=p(n)-4; %(2x+3)3-4 x=roots(q) 结果：x = -1.8969 + 0.6874i -1.8969 - 0.6874i -0.7063 验证： y=polyval(q,x) y = 1.0e-014 * -0.7105 - 0.6217i -0.7105 + 0.6217i 0 2、求方程的正根 解： f=inline(x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x,x); fplot(f,[1,10]);grid on fzero(f,2) ans = 2.1155 第一章习题4：解方程在区间（－2，2）内的根 解： f=inline(x^4-2^x,x); fplot(f,[-2,2]);grid on x1=fzero(f,-1),x2=fzero(f,1) x1 = -0.8613 x2 = 1.2396 4、（超越方程）作出在[-0.1,0.1]内的图，并求附近的解。 x=-0.1:0.01:0.1;y=x.*sin(1./x);plot(x,y) 或者 fplot(x.*sin(1./x),[-0.1,0.1]) 求附近的解：通过下列函数文件，找最接近原点的值 function xi=fun(k0,k1) xi=k0; for x0=k0:1/1000:k1 x1=fzero(x*sin(1/x),x0); if abs(x1)abs(xi) xi=x1; end end 解非线性方程组在原点附近的根 解： f=inline([9*x(1)^2+36*x(2)^2+4*x(3)^2-36,x(1)^2-2*x(2)^2-20*x(3),16*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2],x); [x,f,h]= fsolve(f,[0,0,0]) x = 0.1342 0.9972 -0.0985 f = 1.0e-008 * 0.7690 -0.0418 -0.1054 h = 1 6、求解下列方程组在区域内的解 f=inline([x(1)-0.7*sin(x(1))+0.2*cos(x(2)),x(2)-0.7*cos(x(1))-0.2*sin(x(2))],x); [x,f,h]= fsolve(f,[0.5,0]) x = -0.4442 0.7715 f = 1.0e-011 * -0.0947 0.1749 h = 1 7、（椭圆的交点） close;clear; x=-2:0.01:5;y=sqrt(5-(x-2).^2)+3-2*x;plot(x,y);axis([-1 5 -6 6]);hold on;grid on; y=-sqrt(5-(x-2).^2)+3-2*x;plot(x,y);hold on; y=sqrt(4-2*(x-3).^2)*3;plot(x,y);hold on;