数值分析方法(基于Matlab实现)
潇潇雨歇_
www
展开
-
数据的可视化绘图(Matlab实现)
1.在区间[-4,4]上绘制函数y=sin x 图形代码:>> x=-4:0.1:4;>> y=sin(x);>> plot(x,y)>> gird on %显示网格线 输出如图:2.绘制三维曲线图像 根据下列函数绘图: | x=t sin t, | y=t cos t, 0 | z=t,Matlab代码:原创 2017-03-12 12:02:36 · 4158 阅读 · 0 评论 -
Matlab 数值计算泰勒公式
1.利用n阶泰勒多项式展开来估算函数f(x)=e^x在指定点x处的近似值。要求绝对误差在最大阶数MAXN以内达到给定精度EPS。 taylorcal函数:function taylorcal()eps=0.00001;MAXN=20;x=input('存储输入的浮点数/向量');f1=zeros(1,length(x));for j=1:length(x) xzs=1原创 2017-04-09 12:03:24 · 16371 阅读 · 0 评论 -
Matlab 数值计算----二分法求非线性方程组
bisect.m函数(二分法的实现算法)function[x_star,index,it]=bisect(fun,a,b)%求解非线性计算方程的二分法,其中,fun(x)为需要求根的函数;%a,b为初始区间的端点;%ep为精度,当(b-a)/2<ep时,算法能终止计算,%缺省值为1e-5;%当x_star迭代成功时,输出方程的根%当x_start迭代失败时,输出两端点的值;%in原创 2017-04-09 12:04:01 · 8228 阅读 · 0 评论 -
Matlab 数值计算----迭代法计算非线性方程组在指定区间的根
iterate.m函数function [x_star,index,it]=iterate(phi,x,ep,it_max)if nargin<4 it_max=100;endif nargin<3 ep=1e-5;endindex=0;k=1;while k<it_max x1=x;x=feval(phi,x); if abs(x-x1)<ep ind原创 2017-04-09 12:04:47 · 4657 阅读 · 0 评论 -
Matlab 数值计算----斯特芬森加速迭代法
steffensen.m函数function[x_star,index,it]=steffensen(phi,x,ep,it_max)%斯特芬森加速迭代方法% x为初始点% ep为精度,当| x(k)-x(k-1) |<ep时,终止计算,缺省值为1e-5% it_max为最大迭代次数% x_star为当迭代成功时,输出方程的根% 当迭代失败时,输出最后的迭代值;% index原创 2017-04-09 12:05:28 · 9626 阅读 · 1 评论 -
Matlab 数值计算----牛顿迭代法
Newton.m函数function [x_star,index,it] = Newton(fun,x,ep,it_max)%求解非线性方程的牛顿法%第一个分量是函数值,第二个分量是导数值% x为初始点% ep为精度,当 | x(k)-x(k-1) |<ep时,终止计算,缺省值为1e-5% it_max为最大迭代次数,缺省值为100% x_star为当迭代成功时,输出方程的根%原创 2017-04-09 12:05:41 · 59688 阅读 · 3 评论 -
Matlab 数值计算----牛顿法解非线性方程组
Newtons.m函数function[x_star,index,it]=Newtons(fun,x,ep,it_max)%求解非线性方程组的牛顿法,其中,fun(x)为需要求根的函数%第一个分量是函数值,第二个分量是导数值%x为初始向量(列向量)% ep为精度,当 | |x(k)-x(k-1)| |<ep时,终止计算,缺省值为1e-5% it_max为最大迭代次数,缺省值为100原创 2017-04-09 12:05:52 · 15640 阅读 · 2 评论