有关最近matlab计算的记录

1. matlab画指定半径的圆

之前画圆是

file=dir('./*.txt');   
%clear all
    s=load(file(jj).name);
    x=s(:,2); 
    y=s(:,3);

set(0,'defaultfigurecolor','w') 
set(0,'DefaultFigureVisible', 'off')
plot(x,y,'blue.','markersize',90);

甚至一昧的追求使用markersize来实现半径，markersize单位是dot是1/72inch，简单来说，100markersize是一个格

最终使用

rectangle('Position',[overlap_aa-r,overlap_bb-r,2*r,2*r],'Curvature',[1,1],'linewidth',1,'EdgeColor', 'b'),axis equal

来画圆，画多个圆可使用for循环。rectangle是画矩形的函数，把曲率设置成1，画出来的是圆。参考链接是这里

2. 使用solve和fsolve对方程进行求解

一开始我使用的是solve，过程如下

syms xt yt m n x1 y1
assume(a-R<=x1<=a+R)
assume(b-R<=y1<=b+R)

eq1 = 。。。==1;
eq2 = 。。。==x1;
eq3 = 。。。==y1;
eq4 = 。。。==1;

[xt, yt,x1,y1] = solve(eq1, eq2,eq3,eq4,xt, yt,x1,y1);
res=~isempty(x1);

使用assume可以对自变量的范围进行限定，结果也让人满意，代码也比较简单，但是速度太慢了

导师说，可以使用fsolve，我的最终代码如下

options = optimset('MaxFunEvals',1e3);
[res_x,fval,exitflag]=fsolve(@fxy,[(overlap_a+R)*100,(overlap_b+R)*100,(overlap_a+R)*100,(overlap_b+R)*100]',options);
res_x=res_x./100;

function q=fxy(pp)
xt=pp(1);yt=pp(2);
x1=pp(3);y1=pp(4);

global overlap_etat overlap_aa overlap_bb overlap_a overlap_b R p;
 q(1) =。。。;
 q(2) =。。。;
 q(3) =。。。;
 q(4) =。。。;
end

注意，solve等式是eq=x^2==1，而fsolve是eq=x^2-1。

fsolve使用的是最小二乘法和迭代，公式中[(overlap_a+R)*100,(overlap_b+R)*100,(overlap_a+R)*100,(overlap_b+R)*100]这部分是开始计算的初始值，因为我要计算的数值太小了，就全乘以100，这样放大误差，可以使最后的结果更完美。fsolve不能像solve使用assume那样限制自变量的范围，但是（只是针对我使用来说） 比如自变量范围是[0,2]，在其他范围内求不出来更精确的解（意思是有解，但是精确度不好），我就设定初始值是2。（？这一步我好像设置成1会更好？）

optimset('MaxFunEvals',1e3)是设置了迭代次数，一般是400次迭代，选取误差最小的。但是这里我设置成了1000次

[res_x,fval,exitflag]中

res_x是求解的数值。

fval是函数值 我感觉更像是误差值

exitflag1、一阶最bai优性条件满足容du许范围2、X的变化小于zhi容许范围3、目标函数的变化小于容许范围4、重要搜索方向小于规定的容许范围并且约束违背小于options.TolCon5、重要方向导数小于规定的容许范围并且约束违背小于options.TolCon0、到达最大迭代次数或到达函数评价-1、算法由输出函数终止-2、无可行点

fsolve的速度真的是非常快了。

补充链接在这里