Matlab协方差矩阵分解法生成随机场

Matlab协方差矩阵分解法生成随机场

相关系数矩阵

在这里插入图片描述

% function out=cohesion(x,y,mu,theta)
% end
% xy=load('F:\Research-OUC\基于机器许学习模型的海底斜坡可靠度研究\基于comsol的斜坡稳定性分析\comsol网格操作\grid_operate-matlab.mphtxt');
% x=xy(:,1);
% y=xy(:,2);
X=linspace(0,10,20);
Y=linspace(0,5,10)
[xxx,yyy]=meshgrid(X,Y);x=reshape(xxx,200,1);
y=reshape(yyy,200,1);% 自相关距离
Lh = 10; % 根据需要设置
Lv = 5; % 根据需要设置
% 计算自相关函数矩阵
n = length(x);
autocorr_matrix = zeros(n, n);
for i = 1:nfor j = 1:nautocorr_matrix(i, j) = exp(-((x(i) - x(j))^2 / Lh) - ((y(i) - y(j))^2 / Lv));end
end
% 对相关函数进行chelosky分解
R=chol(autocorr_matrix,'upper');
random_m=randn(length(R));
G=dot(random_m,R);
mu=1e7;
theta=3e6;
thetaln=sqrt(log(1+(theta/mu)^2));
muln=log(mu)-thetaln^2/2;Rf=exp(thetaln*G+muln);
%组合成多列
E_modulus=[x,y,reshape(Rf,200,1)];
% save('E_field.txt', 'E_modulus', '-ascii');% contour(x,y,Rf)
% 创建网格
xlin = linspace(min(x), max(x), 100); % 创建 x 轴的网格点
ylin = linspace(min(y), max(y), 100); % 创建 y 轴的网格点
[X, Y] = meshgrid(xlin, ylin); % 创建网格% 插值
Z = griddata(x, y, Rf, X, Y, 'natural'); % 使用 natural 插值方法适用于不规则数据% 绘制等值线图
figure;
pcolor(X, Y, Z); % 绘制等值线,20 表示等值线的数量
colorbar; % 显示颜色条
title('不规则数据的等值线图');
xlabel('X 轴');
ylabel('Y 轴');%对自相关函数求解特征值和特征函数% % Nkl = 50;
% % [eigenVec,eigenVal]=eigs(autocorr_matrix,Nkl);
% [eigenVec,eigenVal]=eigs(autocorr_matrix);% % % 保存自相关函数矩阵到文件
% % save('autocorr_matrix.txt', 'autocorr_matrix', '-ascii');
% 
% % 绘制自相关函数矩阵图
% figure;
% imagesc(autocorr_matrix);
% colorbar;
% title('自相关函数矩阵');
% xlabel('Index');
% ylabel('Index');

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

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

相关文章

游戏开发面试题1

C#,泛型、反射? C# 泛型是指可以让开发者在声明和使用类和方法时指定类型参数的编程技术。它可以减少代码重复,使程序的可维护性更高,并且能够更好地管理内存。 反射是指.NET程序集中类、方法、属性和字段的可编程性。它允许在不知…

等保测评推动哈尔滨数字化转型中的安全保障

在数字经济的浪潮下,哈尔滨作为东北老工业基地的核心城市,正积极推动数字化转型,以创新技术驱动产业升级和经济发展。网络安全等级保护测评(简称“等保测评”)作为国家网络安全战略的重要组成部分,为哈尔滨…

泛型(generics)

目录 一、泛型概述 二、定义泛型类 三、定义泛型方法 四、List接口中的泛型是如何定义的 五、泛型通配符 六、 泛型通配符的上限和下限 七、可变参数 八、可变参数的使用 九、泛型擦除机制 一、泛型概述 Java泛型是JDK5中引入的一个新特性,提供了编译时类型…

Adobe Premiere Pro 2024 v24.5 (macOS, Windows) - 专业视频编辑软件

Adobe Premiere Pro 2024 v24.5 (macOS, Windows) - 专业视频编辑软件 Acrobat、After Effects、Animate、Audition、Bridge、Character Animator、Dimension、Dreamweaver、Illustrator、InCopy、InDesign、Lightroom Classic、Media Encoder、Photoshop、Premiere Pro、Adob…

在Ubuntu 16.04上安装和配置Redis的方法

前些天发现了一个巨牛的人工智能学习网站,通俗易懂,风趣幽默,忍不住分享一下给大家。点击跳转到网站。 简介 Redis 是一个内存中的键值存储,以其灵活性、性能和广泛的语言支持而闻名。在本指南中,我们将演示如何在 U…

【Selenium配置】WebDriver安装浏览器驱动(ChromeEdge)

【Selenium配置】WebDriver安装浏览器驱动(Chrome&Edge) 文章目录 【Selenium配置】WebDriver安装浏览器驱动(Chrome&Edge)Chrome确认Chrome版本下载对应driver把解压后的chromedriver文件放在chrome安装目录下&#xff0…

算法简介:什么是算法?——定义、历史与应用详解

引言 在现代计算机科学中,算法是一个核心概念。无论是编程还是数据分析,算法都扮演着至关重要的角色。在这篇博客中,我们将深入探讨算法的定义、历史背景以及它在计算机科学中的地位和实际应用。 什么是算法? 算法是解决特定问题…

Linux 定期自动修改密码,防止过期

背景 😂 Jenkins 中使用 ssh 密码登录,实现项目前端部署。 😂 ssh 账号(假设叫做 ft)每隔三个月就会过期,就会导致前端部署报错,提示“Your password is expired” 😂 每次修改密码,需要申请权限,审批,耗时较长。=》希望自动修改密码,防止过期 分析 1、模拟生…

表格布局的概念与属性

表格布局的概念与属性 表格布局(TableLayout)是以行、列的形式来管理控件的,类似与表格。如图所示,是一个表格布局。TableLayout继承自LinearLayout,支持LinearLayout所支持的全部属性,默认为垂直方向的Li…

在线制作网页PHP源码+IAPPv3源码

好久以前刚学PHP的时候搞的,现在翻出来并修复了一下就分享出来了,希望可以给初学者一点借鉴学习例子,虽然不是很好哈,但拿来耍耍应该还算可以吧。 使用教程 环境是PHP就行,直接把源码上传到服务器(虚拟主…

React-Native中关于图片问题知识总结

系列文章目录 React-Native环境搭建(IOS)React-Native项目 — 关于IOS知识储备React-Native项目工程搭建(开发模板搭建)React-Native项目矢量图标库(react-native-vector-icons)React-Native项目 — 自定义…

Transformer-LSTM预测 | Matlab实现Transformer-LSTM多变量时间序列预测

Transformer-LSTM预测 | Matlab实现Transformer-LSTM多变量时间序列预测 目录 Transformer-LSTM预测 | Matlab实现Transformer-LSTM多变量时间序列预测效果一览基本介绍程序设计参考资料 效果一览 基本介绍 1.Matlab实现Transformer-LSTM多变量时间序列预测,Transf…

JavaWeb开发基础7个Web术语解析

7个Web术语 Website: static vs dynamic HTTP HTTP Requests GET vs POST Servlet Container Server: Web vs Application Content Type Website: static vs dynamic 网站内容包括文本、图片、音频、视频,通过URL来访问。网站分为静态网站和动态网站。 静态网…

【计算机组成原理实验】——实验 MIPS 指令系统和 MIPS 体系结构

实验 MIPS 指令系统和 MIPS 体系结构 一、实验目的 了解和熟悉指令级模拟器。 熟练掌握 MIPSsim 模拟器的操作和使用方法。 熟悉 MIPS 指令系统及其特点,加深对 MIPS 指令操作语义的理解。 熟悉 MIPS 体系结构。 二、实验平台 实验平台采用指令级和流水线操作…

通用的职位招聘小程序ui模板

蓝色简单的校园招聘,行业招聘,职位招聘手机小程序页面模板。包含:职位列表、职位详情、基本信息填写、登录、个人主页、消息页面等功能 通用的职位招聘小程序ui模板

洛谷 P3613 学习用map代替大大大数组的好题

题目链接:P3613 【深基15.例2】寄包柜 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn) 题目截图: 题意分析: 非常简单的存入和取出操作 唯一的 “难点” 在于 数组开不到 a[100007][100007],会暴内存 非常巧妙的引入 map 来解决…

Springboot 3.x - Reactive programming

一、Preliminary Knowledge Functional InterfaceLambda expressionStream API Intermediate operation filter:Used to filter elements in a streammap:One-to-one conversionflatMap:One-to-many conversiondistinct、sorted、peek、limi…

详解Java垃圾回收(GC)机制

一、为什么需要垃圾回收 如果不进行垃圾回收,内存迟早都会被消耗空,因为我们在不断的分配内存空间而不进行回收。除非内存无限大,我们可以任性的分配而不回收,但是事实并非如此。所以,垃圾回收是必须的。 二、哪些内…

【ARMv8/v9 GIC 系列 1.8 -- PE 中断处理的前期评估】

请阅读【ARM GICv3/v4 实战学习 】 文章目录 Interaction of group and individual interrupt enablesLPIs的启用Interaction of group and individual interrupt enables 在ARM GICv3和GICv4架构中,GICD_*和GICR_*寄存器组 决定了处理器元素(PE)能够识别的最高优先级 pend…

【Java】了解异常

初始异常 我们平时应该已经接触过一些 “异常” 了,这里列举一些例子。 算术异常: 数组下标越界异常: 访问空指针异常: 所谓异常指的就是程序在 运行时 出现错误时通知调用者的一种机制。 异常的基本用法 捕获异常 try{ 有可能…