读书笔记(一)
这是第一部分——迭代 将代码复制到m文件中即可运行 % 迭代是计算的关键 % % 上键:调用先前的命令 % %下面这个“双%”表示一个section(就是英语听力那个section) %可以把程序分成好多片段,可以分块执行run section %% 长精度显示结果 format long % %% 浮点数通过与eps比较判断 % abs(x-y)< eps或者5e-5等小量 % 即为x=y % %% roots([多项式系数 空格隔开]) % 求多项式=0的根 x1 = roots([1 -1 -1]); %x^2 - x - 1 = 0的根 %% fsolve % fsolve(表达式变量名,不同方法参数-3~4) % fsolve(@自变量 表达式,区间) x2 = fsolve(@(x) sin(3*x),[1 4]); f = @(x) sin(3*x); x3 = fsolve(f,1); x4 = fsolve(f,-1); %% solve(带有符号变量的方程,符号变量) syms x x5 = solve(x + 1/x == 2.5,x); %% 迭代求解 clear clc x=3; while abs(x - (2.5 - 1/x)) > 5e-5 x = 2.5-1/x; end disp(x) %只能求出一个 %好在软件提供大量优化好的可以直接用的函数 %% 另一种迭代 clear clc x=3; for i = 1:8 x = 2.5-1/x; plot(i,x,'o'); axis([0 9 1 3]); hold on end