matlab生成均匀部分散点图,应用halton序列生成均匀散点图

前言

一门课的作业要用RBF-DQ方法计算流场。简单来讲就是要在无网格的条件下用高精度格式实现流场的仿真计算。生成散点图时,用蒙特卡洛的方法得到的散点图分布不是很均匀,于是想到一种叫做halton的序列。本篇将运用halton序列生成均匀的散点图。

halton序列原理

点下面的链接看halton序列的介绍:

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

生成散点图代码实现

为了方便展示,这里用matlab进行编程。实际上halton序列实现起来非常简单,可以在任何编程语言上进行。

首先,建立halton函数,输入为序数和基数。

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;

然后,在主程序中,设置好x和y方向的两个基数,这里分别为2和3(最好是互质的,那样效果更好)

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

结果展示

基数为2和3时:

d6c7b84f572386c4b758a14345bbae83.png

基数为3和5时:

472e6b7d1e2bb957252ec9c5a899e39a.png

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/445764.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

php分页页数过多隐藏,织梦DEDE中内容页分页数目太多时隐藏方法

在dede程序的文章或图集模型,内容页分页过多,假如分页80页,那么dede默认的内容分页标签{dede:pagebreak/}会把所有的分页都显示出来,这样会撑破网页模板,即使使用css控制可以解决这个撑破模板问题,用overfl…

nginx php 跨域访问权限,nginx + php 实现跨域请求填坑笔记

最近自己构建了一个轻量级的 MVC 框架,现在将该MVC用于生产环境的项目中,目前因为项目中有跨域的请求,所以需要做一些CORS授权,但在实际使用出现了一些问题,目前已解决,故做下记录。服务器:cent…

python 内存溢出能捕获吗,Python——写自己的异常(形成强大的捕捉和处理能力)...

"""1、编写自己的异常时,通常是继承Exception2、如果继承的是BaseException,那么就不会被通用的except Exception捕捉到3、自己定义的异常Python不会自动触发,需要自己触发,用raise进行触发4、python只会自动触发它…

解决动态规划问题4步曲

概述 (确定状态)确定问题状态 提炼最后一步子问题转化 (求得方程)转移方程,把问题方程化(设初置界)按照实际逻辑设置初始条件和边界情况(确序再解)确定计算顺序并求解 …

php ajax队列,AJAX请求队列实现

这篇文章主要为大家详细介绍了AJAX请求队列的实现代码,具有一定的参考价值,感兴趣的小伙伴们可以参考一下AJAX在使用的过程中会遇到一个问题,当用户短时间内执行了多个异步请求的时候,如果前一个请求没完成,将会被取消执行最新的一…

php Spreadsheet 导出,PhpSpreadsheet 导出Excel

/*** Excel 助手* sudo composer require phpoffice/phpspreadsheet*/namespace CommonUtil;use PhpOfficePhpSpreadsheetSpreadsheet;use PhpOfficePhpSpreadsheetWriterXlsx;use PhpOfficePhpSpreadsheetStyleAlignment;use PhpOfficePhpSpreadsheetStyleColor;class ExcelUt…

php 不同时区时间转换,在PHP中将DateTime字符串转换为不同的时区

好吧,我有以下代码$from "Asia/Manila";$to "UTC";$org_time new DateTime("2012-05-15 10:50:00");$org_time $org_time->format("Y-m-d H:i:s");$conv_time NULL;$userTimezone new DateTimeZone($from);$gmtTimezone new…

php iis ajax 无效,IIS7中Ajax.AjaxMethod无效的原因及解决方法

最近做用Ajax.AjaxMethod方法的时候,在asp.net的服务器下一切正常,用iis的时候,js中总是cs类找不到,我就郁闷了,折腾了大半天,终于找到错误原因了。因为我发布网站用的是iis7,所以在web.config位…

查看oracle监听服务器,处理Oracle 监听文件listener.log问题

如果连接时候变得较慢 查看Oracle日志记录,可能是因为此文件太大,超过2G,需要定期清理,(如果多用户,记得用root,可能没权限)查看listener.log?find / -name listener.log经查看,竟然…

oracle添加伪列,Oracle伪列 - jifengtang的个人空间 - OSCHINA - 中文开源技术交流社区...

在oracle10g和下,伪列包括如下内容:lHierarchical Query Pseudocolumns分级查询是oracle提供的递归查询语法,在这里不做展开。只有在分级查询下,才可以使用以下伪列:1.CONNECT_BY_ISCYCLE Pseudocolumn2.CONNECT_BY_IS…

美国oracle球场,美国体育馆考察——美国体育产业是如何盈利的?

体育是美国一项较高利润的产业,其发展规模、发展水平和效益都是世界一流的。美国体育馆考察,主要考察美国体育产业的盈利模式和体育赛事的赞助模式以及球馆的运营管理,并对比中美体育产业的差异,从中获得先进的体育产业运营思维&a…

php集成环境怎么打开,PHP集成开发环境PhpStorm快速入门指南(二):打开一个项目...

PhpStorm是一个轻量级且便捷的PHP IDE,其旨在提高用户效率,可深刻理解用户的编码,提供智能代码补全,快速导航以及即时错误检查。可随时帮助用户对其编码进行调整,运行单元测试或者提供可视化debug功能。PhpStorm 2019.…

如何查询oracle最近报警信息,教你怎样用Oracle方便地查看报警日志错误

在网上查了几天的资料,尝试综合清除告警日志内容及建外部表的方式来解决这一问题。一:备份并清除告警日志内容将每天的告警日志备份好,然后进行清除。1:备份报警日志在$ORACLE_HOME/SID/bdump/ 目录下,按日期备份alert_ORACLE_你…

计算机分php,计算机按照处理数据的形态分类,可以分为什么?

计算机按照处理数据的形态分类,可以分为:1、数字计算机,是以数字形式的量值在机器内部进行运算和存储的电子计算机;2、模拟计算机,是根据相似原理,用一种连续变化的模拟量作为被运算的对象的计算机&#xf…

2.oracle物理结构,oracle实验2oracle物理结构管理

oracle实验2oracle物理结构管理 (6页)本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦!9.9 积分实验2 oracle物理存储结构管理、实验目的1. 掌握物理结构的创建和修改方法2. 掌握表空间的存储参数设置方…

linux mount 查看挂载目录,Linux下使用mount来挂载设备到目录

一般情况下直接mount 设备路径 目录路径,就可以了。umount 设备名,就可以卸载这个设备了使用lsblk -f可以查看挂载的设备,以及这些设备的文件系统。roottao-PC:/boot# lsblk -fNAME FSTYPE LABEL UUID MOUNTPOINTsda├─sda1├─sda2 vfat SY…

centos7是哪种版本Linux,centos7怎么查看系统版本是不是7.2 7.5 7.6

CentOS的版本号信息一般存放在配置文件当中,在CentOS中,与其版本相关的配置文件中都有centos关键字,该文件一般存放在/etc/目录下,所以说我们可以直接在该文件夹下搜索相关的文件。其中存放其版本配置信息的文件为“centos-releas…

linux6.0 安装教程,CentOS 6.0安装步骤

1.安装引导选择安装或升级现有系统(Install or upgrade an existing system):这个选项是默认的。 选择此选项,安装到您的计算机使用CentOS的图形安装程序的系统。2.检测光盘介质可以选择skip跳过3.选择安装过程中的语言这里选择chinese中文简…

LeetCode - Easy - 14. Longest Common Prefix

Topic String Description https://leetcode.com/problems/longest-common-prefix/ Write a function to find the longest common prefix string amongst an array of strings. If there is no common prefix, return an empty string “”. Example 1: Input: strs […

linux 虚函数调用性能,C++对象布局及多态实现探索之虚函数调用

我们再看看虚成员函数的调用。类C041中含有虚成员函数,它的定义如下:struct C041{C041() : c_(0x01) {}virtual void foo() { c_ 0x02; }char c_;};执行如下代码:C041 obj;PRINT_DETAIL(C041, obj)PRINT_VTABLE_ITEM(obj, 0, 0)obj.foo();C0…