数学实验
关注
分享
扫一扫
「已注销」
这个作者很懒，什么都没留下…
展开
-
用对分法求解方程x^3-x^2-x-1=0
原理： 每一次将区间压缩至原来的一半而且其中一个区间至少包含一个根，逐步缩短区间，直至最终区间最佳为止。 Matlab程序： function [y]=f(x) y=x^3-x^2-x-1; end function []=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 · 5390 阅读 · 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; end end 运行结果：原创 2015-01-17 14:42:39 · 1099 阅读 · 0 评论 -
对所给的斜率场求解方程，并将所给程序进行注释。
Matlab程序代码： clear s t x0 y0 a b syms s t f=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 · 931 阅读 · 0 评论 -
对所编函数输入x=[1:1:5];y=[5:1:9];z=[1:0.1:5];求出对应的拉格朗日插值函数
Matlab程序代码： %给出两组数据数据x,y求拉格朗函数 function[]=laggrange(x,y,z) hold on syms t plot(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 · 1861 阅读 · 0 评论