前言

一門課的作業(yè)要用RBF-DQ方法計算流場。簡單來講就是要在無網(wǎng)格的條件下用高精度格式實現(xiàn)流場的仿真計算。生成散點圖時，用蒙特卡洛的方法得到的散點圖分布不是很均勻，于是想到一種叫做halton的序列。本篇將運用halton序列生成均勻的散點圖。

halton序列原理

點下面的鏈接看halton序列的介紹：

https://blog.csdn.net/Amber_amber/article/details/47421053

生成散點圖代碼實現(xiàn)

為了方便展示，這里用matlab進(jìn)行編程。實際上halton序列實現(xiàn)起來非常簡單，可以在任何編程語言上進(jìn)行。

首先，建立halton函數(shù)，輸入為序數(shù)和基數(shù)。

function halton=halton(index,base)

result=0;

f=1/base;

i=index;

while(i>0)

result=result+f*mod(i,base);

i=floor(i/base);

f=f/base;

end

halton=result;

然后，在主程序中，設(shè)置好x和y方向的兩個基數(shù)，這里分別為2和3(最好是互質(zhì)的，那樣效果更好)

index=1000;

xy=zeros(index,2);

base1=3;

base2=2;

for i=1:index

xy(i,1)=halton(i,base1);

xy(i,2)=halton(i,base2);

end

結(jié)果展示

基數(shù)為2和3時：

基數(shù)為3和5時：

總結(jié)

以上是生活随笔為你收集整理的matlab生成均匀部分散点图,应用halton序列生成均匀散点图的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。