数学实验
「已注销」
这个作者很懒,什么都没留下…
展开
-
用对分法求解方程x^3-x^2-x-1=0
原理:每一次将区间压缩至原来的一半而且其中一个区间至少包含一个根,逐步缩短区间,直至最终区间最佳为止。Matlab程序:function [y]=f(x)y=x^3-x^2-x-1;endfunction []=bina(a,b)while abs(a-b)>0.1 if (f(a)*f(b))<0 c=(a+b)/2; if (f(a原创 2015-01-17 14:43:43 · 5305 阅读 · 1 评论 -
编写SQRT.m文件以实现求一个数的平方根
方法:实验代码:function y=SQRT(x,n)z=x;for i=1:n y=(1/2)*(x+z/x); if (x-y)>=(10^(-5)) || (y-x)>=(10^(-5)) x=y; endend 运行结果:原创 2015-01-17 14:42:39 · 1073 阅读 · 0 评论 -
对所给的斜率场求解方程,并将所给程序进行注释。
Matlab程序代码:clear s t x0 y0 a bsyms s tf=sin(s)*sin(t);a=16.0;b=16.0;x0=-8;y0=-8;m=40;n=40;h1=a/m;h2=b/n;hold % 将新旧图形保持在同一个图中for i=1:m % i 从1每次以1为距离差一直加到m s=x0+(i-1)*h原创 2015-01-17 14:52:25 · 902 阅读 · 0 评论 -
对所编函数输入x=[1:1:5];y=[5:1:9];z=[1:0.1:5];求出对应的拉格朗日插值函数
Matlab程序代码:%给出两组数据数据x,y求拉格朗函数function[]=laggrange(x,y,z)hold onsyms tplot(x,y,'r:o')m=length(x)n=length(z) La=0; for i=1:m L=1; for j=1:m if i~=j原创 2015-01-17 14:56:57 · 1843 阅读 · 0 评论