锦到说MATLAB
高开低走。
这个作者很懒,什么都没留下…
展开
-
数学实验第四课:不定方程求解与综合实验
蒙特卡罗求积分值近似求解0到1上e^x的积分a=0;b=1;m=20000%迭代20000次s=0;H=exp(1);%设置为落在曲边题型内的点数for i=1:mxi=rand();yi=H*rand();if yi<exp(xi)s=s+1;end;end;fprintf('%s%g\n','exp(x)在[0,1]上的积分约等于',H*(b-a)*s/m)这样算的结果比较接近e-1...原创 2020-08-05 15:34:27 · 224 阅读 · 0 评论 -
数学实验第三课:线性映射的迭代与特征向量的计算
线性映射的迭代列出线性方程组A=[0.5,0.3,0.3;0.3,0.6,0.3;0.2,0.1,0.4];n=10;x(:,1)=[200,300,300]';%按200,300,300来分,但最后的结果与初始值无关for k=1:n%循环x(:,k+1)=A*x(:,k);endround(x)ans = 1 至 5 列 200 280 296 299 300 300 330 339 342 343 300 190原创 2020-07-26 23:30:39 · 1165 阅读 · 0 评论 -
数学实验第二课:函数迭代
迭代定理在迭代函数连续的条件下,如果迭代序列收敛,则一定收敛于方程x=f(x)的根x*,这个根为函数的不动点。不动点处的一阶导数小于1,则称不动点x为稳定的,若一阶导数等于0,则称不动点为超稳定,在超稳定附近,迭代过程x(n+1)=f(x(n))收敛到x*是非常快的。用迭代法近似计算计算三次根号2过程分析首先取f(x)=2/x^2因为f(2(1/3))=2(1/3)在该点的导数=2>1,不满足要求>> syms x;>> y=diff(2/x^2);原创 2020-07-26 15:39:04 · 3108 阅读 · 0 评论 -
数学实验第一章:MATLAB基础
基本操作命令命令该命令的功能命令该命令的功能cd显示或改变工作目录hold图形保持命令clear清除内存变量pack整理内存碎片clf清除图形窗口path显示搜索目录diary日志文件命令quit退出MATLABdir显示当前目录下文件save保存内存变量到指定文件disp显示变量或文字内容type显示文件内容常用标点符号标点符号定义标点符号定义;区分行,取消运行显示灯·小数点以原创 2020-06-30 10:49:23 · 5091 阅读 · 1 评论