a=rand(100,1); %生成点的x坐标
b=rand(100,1)+1;%生成点的Y坐标
plot(a,b,'.') %画点
y)
dotx=t-y,y]=ode45(',1)
用plot(t;dy=x-y',3],你要哪点值都可以,即使用插值得到的也是近似值;myfun'exp(x)
-
1
把x=1;,')
得到
y=x
+
2/,[0;
%%%%%%%%%%%%%%%%m-file
以上要单存一个文件在工作路径下%%%%%%%%%%%%%在command
窗口
[t,',直接在图上看就可以了,因为你选择ode45已经是近似了;x',要想要精确值建议求解析解如下
y=dsolve(',y)
绘图后;y(0)=1',2%%%%%%%%%%%%%%%%m-file
以下要单存一个文件在工作路径下%%%%%%%%%%%%%function
dotx=myfun(t;
n=100000;
x=rand(1,n);
y=rand(1,n)*2;
plot(x,y,'.')
%%%计算相同面积内的点数
dx=0.1;
dy=0.2;
[xx,yy]=meshgrid(0:dx:1,0:dy:2);
[noy,nox]=size(xx);
nx=ceil((x+eps)/dx);
ny=ceil((y+eps)/dy);
nn=accumarray([nx(:),ny(:)],ones(n,1));
nn=[nn,zeros(size(nn,1),nox-size(nn,2));zeros(noy-size(nn,1),nox)];
figure
pcolor(xx,yy,nn);
shading flat
colorbar %如果均匀的话,这里只有一个颜色,或色彩相当接近
max(max(nn(1:end-1,1:end-1)))-min(min(nn(1:end-1,1:end-1))) %这个数最好是零