matlab求QR分解，QL分解，LQ分解和RQ分解
摘要：Matlab的命令里只有QR分解，如何求QL分解和LQ,RQ分解？LQ分解： [Q 1,R 1]= qr (H'); H=(Q 1*R 1)'=R 1' *Q 1'=L*Q; 所以下三角阵L=R 1',正交阵Q=Q 1'QL分解：[Q 1,R 1]= qr ( inv(H' ) ); H'=(Q 1
阅读全文
posted @ 2023-10-22 21:57
随笔分类 - MATLAB
LTI系统，已知系统输入和输出，求脉冲响应
摘要：https://www.ilovematlab.cn/forum.php?mod=viewthread&tid=298897&ordertype=1&_dsign=961ef8cf 问题描述：一个LTI系统，输入序列为x(n)，输出序列是y(n),均已知。求系统的脉冲响应。我的思路是这样的：设系统响
阅读全文
posted @ 2023-05-17 09:38
matlab生成指定均值向量和协方差矩阵的多维正态分布样本
摘要：产生一个协方差矩阵为R的n维随机正态分布的一组样本，matlab没有现成的函数，不过我们可以通过一个线性变换来实现。 我们知道，matlab产生的n维正态样本中的每个分量都是相互独立的，或者说，它的协方差矩阵是一个数量矩阵mI，如：X = randn(10000,4);产生10000个4维分布的正态
阅读全文
posted @ 2021-12-15 20:07
matlab追加写入txt
摘要：save([path, 'data.txt'], 'data', '-ascii', '-append')
阅读全文
posted @ 2021-10-21 23:06
matlab histfit(data,nbins,dist)直方图拟合
摘要：rng default; % For reproducibility r = normrnd(10,1,100,1); figure histfit(r) rng default; % For reproducibility r = betarnd(3,10,100,1); figure histf
阅读全文
posted @ 2021-09-18 16:01
matlab normspec叠加图
摘要：p = normspec(specs,mu,sigma) matlab自带的这个函数每次绘制都会打开一个新的句柄，无法绘制两个分布的叠加。在原函数基础上修改获得如下图类似的效果。 normspec两个分布叠加 修改 function [p,h] = mynormspec(specs,mu,sigma
阅读全文
posted @ 2021-09-18 15:45
PCA Biplot含义
摘要：PCA biplot = PCA score plot + loading plot matlab PCA分析命令： [coeff,score,latent,tsquared,explained,mu] = pca(X); 前两个主成分的PCA Biplot 如上图所示，PCA Biplot包含两部
阅读全文
posted @ 2021-09-03 10:21
matlab和python中进行拉丁超立方抽样（逆变换法）
摘要：什么是拉丁超立方抽样法？它和蒙特卡罗模拟有什么关系？ 逆变换方法(The Inverse Transform Method)采样 拉丁超立方抽样 拉丁超立方抽样(LHS)是一种从多维分布中生成参数值的近随机样本的统计方法。抽样方法常用于构建计算机实验或进行蒙特卡罗积分。在统计抽样的上下文中，当(且仅
阅读全文
posted @ 2021-09-01 22:14
matlab PCA分析
摘要：pca( ) 采用matlab自带的函数pca()进行主成分分析 [coeff, score, latent, tsquared, explained, mu] = pca(x) 假设数据x为n行p列的多变量数据，n为观测次数，p为变量维度。 coeff：为PCA变换系数，也称为loadings。
阅读全文
posted @ 2021-08-31 23:10
matlab计算马氏平方距离
摘要：马氏距离计算一个点到一个参考分布的距离。马氏距离的优点是不受各个维度特征尺度的影响（尺度无关）。 d2 = mahal(Y, X) d2 = mahal(Y,X) returns the squared Mahalanobis distance of each observation in Y to
阅读全文
posted @ 2021-08-31 10:13
matlab计算排列组合
摘要：perms(x) perms(x) 给出向量x的所有排列。 perms(1:3) ans = 3 2 1 3 1 2 2 3 1 2 1 3 1 3 2 1 2 3 nchoosek(x,m) nchoosek(x,m)给出了从包含n个元素的向量x中选取m个元素的组合。 nchoosek(1:4,3
阅读全文
posted @ 2021-08-29 15:47
matlab函数tight_subplot控制图像的边界（margin），subplot的间距（gap）
摘要：tight_subplot实现比matlab自带的subplot更紧凑的子图绘制。 function [ha, pos] = tight_subplot(Nh, Nw, gap, marg_h, marg_w) % tight_subplot creates "subplot" axes with
阅读全文
posted @ 2021-08-26 22:43
matlab获取数据直方图和经验累计分布CDF
摘要：histogram获取数据直方图 通过BinWidth参数调整箱子的宽度 BinWidth = 0.01 BinWidth = 0.1 cdfplot获取经验累积分布
阅读全文
posted @ 2021-08-26 17:44
matlab计算数据分位点prctile
摘要：prctile计算数据集的百分位数（分位点） Y = prctile(X,p) 根据区间 [0,100] 中的百分比 p 返回数据向量或数组 X 中元素的百分位数。 a = rand(1, 10); sort(a) prctile(a, 70) a = 0.3452 0.2825 0.7467 0.
阅读全文
posted @ 2021-08-26 17:33
matlab散点矩阵图
摘要：1. gplotmatrix Matrix of scatter plots by group 按组划分的散点图矩阵 gplotmatrix(x,[],group) 创建 x 中数据的散点图矩阵，按 group 中的分组变量分组。 结果图中的每一组单独的轴都包含一个 x 列与另一列 x 的散点图。
阅读全文
posted @ 2021-08-25 22:44
matlab绘图充满图窗
摘要：通常以下命令能满足要求： set(gca,'LooseInset', get(gca,'TightInset')) 但有时上述命令会覆盖一些元素，此时手动控制边距： 宽度方向空白区域0.03， 高度方向空白区域0.05 set(gca,'looseInset',[0 0 0.03 0.05]);
阅读全文
posted @ 2021-04-29 10:41
matlab遍历文件夹
摘要：遍历文件加下指定后缀的文件 通过fullfile + dir命令完成 dirOutput = dir(fullfile(path,'*.fig')); % 遍历文件夹下所有.fig文件 提取文件名 fileNames = {dirOutput.name}; % 得到带后缀的文件名cell 将文件名与
阅读全文
posted @ 2021-04-27 14:49
matlab更改legend中marker的大小
摘要：Ps：更简单的方式 可以手动添加一组自定义的数据（X, Y），把这组数据的X或Y坐标与自己的数据错开。通过句柄对随便添加的自定义的数据设置markersize。然后通过xlim或ylim把手动添加的点遮掉即可。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
阅读全文
posted @ 2021-04-26 15:40
matlab信号小波分解与重构入门
摘要：%% 1. 利用MATLAB产生分解与重构滤波器组 % [Ld, Hd, Lr, Hr] = wfilters(wn); % wfname:小波名 % Ld:分解低通滤波器h0[-n]; % Hd:分解高通滤波器h1[-n]; % Lr:分解低通滤波器h0[-n]; % Hr:分解高通滤波器h1[-
阅读全文
posted @ 2021-04-09 16:55
matlab绘制分组boxplot图
摘要：matlab采用boxplot绘制分组箱形图时不是很方便，这里给出一种实现方法。 假设ratio1~ratio9每个矩阵的形状均为[N, L]，其中，L为一组数据中有几个类别，N为每个类别的数据量。共有9组数据（ratio1~ratio9）。 MRL2 = {ratio1, ratio2, rati
阅读全文
posted @ 2021-04-02 16:34
