- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一、编写函数作图的程序,通过图形求出方程的近似解。
1、 clear
syms x y
x=-10:0.1:10;
y=2.*x.*sin(x)-3;
plot(x,y);
grid
gtext(y=2.*x.*sin(x)-3)
2、 clear
syms x y1 y2
x=0:0.1:10;
y1=2.*sin(x);
y2=3./x;
plot(x,y1,x,y2);
grid
gtext(y1=2sin(x),y2=3/x)
3、 clear
syms x y
x=-2:0.1:4;
y=4.*x.^5-8.*x.^4-26.*x.^3+30;
plot(x,y);
grid
gtext(y=4*x^5-8*x^4-26*x^3+30)
x=solve(4*x^5-8*x^4-26*x^3+30=0,x);
x1=double(x)
x1 =
1.0000
3.7117
-1.9244
-0.3936 + 0.9461i
-0.3936 - 0.9461i
二、用逐步搜索的方法求解。
function [k,r]=zhubuss(a,b,h,tol)
X=a:h:b;Y=funs(X);
n=(b-a)/h+1;m=0;
X(n+1)=X(n);Y(n+1)=Y(n);
for k=2:n
X(k)=a+k*h;
Y(k)=funs(X(k));
sk=Y(k)*Y(k-1);
if sk=0
m=m+1;
r(m)=X(k);
end
xielv=(Y(k+1)-Y(k))*(Y(k)-Y(k-1));
if (abs(Y(k))tol)(xielv=0)
m=m+1;
r(m)=X(k);
end
end
(1)function y=funs(x)
y=2.*x.^3+2.*x.^2-3.*x-3
[k,r]=zhubuss(-2,2,0.001,0.0001)
k =
4001
r =
-1.2240 -1.0000 -1.0000 -0.9990 1.2250
(2)function y=funs(x)
y=sin(cos(2.*x.^3))
[k,r]=zhubuss(-2,2,0.001,0.0001)
k =
4001
r =
-1.9190 -1.7640 -1.5770 -1.3300 -0.9220
0.9230 1.3310 1.5780 1.7650 1.9200
三、用二分法求解。
1、 clear
syms x y
x=-4:0.1:4;
y=x.^3-x+4;
plot(x,y)
grid;gtext(y=2*x^3+2*x^2-3*x-3)
2、 x=-4:1:4,y=x.^3-x+4
x =
-4 -3 -2 -1 0 1 2 3 4
y =
-56 -20 -2 4 4 4 10 28 64
3、function [k,x,wucha,yx]=erfen(a,b,abtol)
a(1)=a;b(1)=b;
ya=fun(a(1));yb=fun(b(1));
if ya*yb0
disp(ya*yb0,重新调整区间端点a,b.)
return
end
max1=-1+ceil((log(b-a)-log(abtol))/log(2));
for k=1:max1+1
a;ya=fun(a);
b;yb=fun(b);
x=(a+b)/2;yx=fun(x);
wucha=abs(b-a)/2;k=k-1;
[k,a,b,x,wucha,ya,yb,yx]
if yx==0
a=x;b=x;
elseif yb*yx0
b=x;yb=yx;
else
a=x;ya=yx;
end
if b-aabtol
return
end
end
k=max1;x;wucha;yx=fun(x);
function y=fun(x)
y=x.^3-x+4
[k,x,wucha,yx]=erfen(-2,-1,0.001)
k =
9
x =
文档评论(0)