matlab练习程序(数据主方向)(代码片段)

tiandsp tiandsp     2022-12-13     378

关键词:

计算方法还是svd,这里把两个轴显示出来了。

过去用svd做过pca,就是在此基础上多一步运算,把原始数据旋转到主轴对应的坐标系上即可。

matlab代码如下:

clear all;
close all;
clc;
 
mu = [0 0];
sigma = [1 1.3;1.3 3];
data = mvnrnd(mu,sigma,1000);
 
me = mean(data);
sigma =(data-me)*(data-me) / (length(data)-1);
 
plot(data(:,1),data(:,2),.);
hold on;
[u,s,v] = svd(sigma);
 
plot([0 u(1,1)*5], [0 u(2,1)*5]);  %乘个5好看
plot([0 u(1,2)*5], [0 u(2,2)*5]);
 
axis([-6 6 -6 6])
  

结果:

技术图片

matlab|如何绘制二维散点主方向直方图(代码片段)

简简单单画个二维散点图主方向上的直方图:工具函数输入nx2大小的数组,返回以下四个图形对象:H1:‘数据点’H2:‘主方向’H3:‘数据范围’H4:‘统计图’function[H1,H2,H3,H4]=PCA2Hist(X,varargin)%H1:'数据点'%H2:'主方... 查看详情

matlab|如何绘制二维散点主方向直方图(代码片段)

简简单单画个二维散点图主方向上的直方图:工具函数输入nx2大小的数组,返回以下四个图形对象:H1:‘数据点’H2:‘主方向’H3:‘数据范围’H4:‘统计图’function[H1,H2,H3,H4]=PCA2Hist(X,varargin)%H1:'数据点'%H2:'主方... 查看详情

matlab练习程序(神经网络识别mnist手写数据集)(代码片段)

记得上次练习了神经网络分类,不过当时应该有些地方写的还是不对。这次用神经网络识别mnist手写数据集,主要参考了深度学习工具包的一些代码。mnist数据集训练数据一共有28*28*60000个像素,标签有60000个。测试数据一共有28*2... 查看详情

matlab练习程序(神经网络识别mnist手写数据集)(代码片段)

记得上次练习了神经网络分类,不过当时应该有些地方写的还是不对。这次用神经网络识别mnist手写数据集,主要参考了深度学习工具包的一些代码。mnist数据集训练数据一共有28*28*60000个像素,标签有60000个。测试数据一共有28*2... 查看详情

sql刷题day18----sql汇总数据专项练习(代码片段)

博主昵称:跳楼梯企鹅创作初心:本博客的初心为与技术朋友们相互交流,每个人的技术都存在短板,博主也是一样,虚心求教,希望各位技术友给予指导。博主座右铭:发现光,追随光,成为... 查看详情

matlab练习程序(最小二乘多项式拟合)(代码片段)

最近在分析一些数据,就是数据拟合的一些事情,用到了matlab的polyfit函数,效果不错。因此想了解一下这个多项式具体是如何拟合出来的,所以就搜了相关资料。这个文档介绍的还不错,我估计任何一本数值分析教材上讲的都... 查看详情

matlab练习程序(椭球拟合)(代码片段)

这次我们来拟合一个椭球,之前也拟合过空间的椭圆,不过当时只用了五个点,方程组应该是欠定的,看看就好。要拟合椭球,首先设定椭球一般方程:根据这个方程和已有的空间椭球点数据,利用最小二乘就能得到上面九个参... 查看详情

matlab练习程序(修正指数曲线拟合)(代码片段)

对于一般的指数曲线如:y=a*e^(k*t),可以先对两边求对数得到:log(y)=log(a)+k*t这样的曲线,然后用最小二乘来计算系数。但是对于修正指数曲线如:y=k+a*b^t这样的函数,没法直接求对数然后用最小二乘,因为有一个常数项k,这里... 查看详情

matlab练习程序(快速搜索随机树rrt)(代码片段)

...小于一定阈值,搜索结束,根据前驱列表画出搜索路径。matlab代码如下:main.m:clearall;closeall;clc;img=imread(‘map.png‘);%空间地图imshow(img);holdon;[h,w]=size(img);p=ginput();%选取起始与结束位置plot(p(:,1),p(:,2),‘r.‘);pc=p(1,:);%随机节点 查看详情

matlab练习程序(简单多边形的核)(代码片段)

...计算下一个点,循环第2步直到遍历所有点。结果如下:matlab代码如下: clearall;closeall;clc;n=20;p=rand(n,2);p=createSimplyPoly(p);n=n+1;p(n,:)=p(1,:);maxX=max(p(:,1));minX=min(p(:,1));maxY=max(p(:,2));minY=min(p(:,2));core=[minXminY;%初始化核minXmaxY;maxXmaxY;maxXm... 查看详情

matlab练习程序(贝塞尔曲线)(代码片段)

下面三个公式分别是一次、二次和三次贝塞尔曲线公式:通用的贝塞尔曲线公式如下:可以看出,系数是由一个杨辉三角组成的。这里的一次或者二次三次由控制点个数来决定,次数等于控制点个数-1。实现的效果如下:代码如... 查看详情

matlab练习程序(渲染三原色)(代码片段)

这里我用的空间是x向右为正,y向下为正,z向屏幕里面为正。相当于标准右手系绕x轴旋转了180度。将三个点光源放在r=[0.3,0,0.5];g=[0.3,-0.5*cos(pi/6),-0.5*sin(pi/6)];b=[0.3,0.5*cos(pi/6),-0.5*sin(pi/6)];这三个位置上,向四周发射光线,取光线到... 查看详情

matlab利用主成分分析进行综合评价排名(代码片段)

文章目录问题提出数据问题的求解数据预处理编写程序求解问题提出利用主成分分析对2007年安徽省17个地区的经济发展进行分析,给出排名;此时能否只用第一主成分进行排名,为什么?数据数据采用安徽省各地... 查看详情

matlab练习程序(局部加权线性回归)(代码片段)

通常我们使用的最小二乘都需要预先设定一个模型,然后通过最小二乘方法解出模型的系数。而大多数情况是我们是不知道这个模型的,比如这篇博客中z=ax^2+by^2+cxy+dx+ey+f这样的模型。局部加权线性最小二乘就不需要我们预先知... 查看详情

matlab练习程序(三种方法解最小二乘)(代码片段)

解最小二乘的方法有很多,这里给出常见的三种方法实现。一是一般方法,之前博客一般都用这种方法。二是svd分解法,之前有用过svd(见这里,这里,这里和这里)解其他问题,但是没用来解过最小二乘。三是qr分解法,这个好... 查看详情

sql刷题day19----sql联结表专项练习(代码片段)

博主昵称:跳楼梯企鹅创作初心:本博客的初心为与技术朋友们相互交流,每个人的技术都存在短板,博主也是一样,虚心求教,希望各位技术友给予指导。博主座右铭:发现光,追随光,成为... 查看详情

sql刷题day12----sql汇总数据专项练习(代码片段)

博主昵称:跳楼梯企鹅博主主页面链接:博主主页传送门博主专栏页面连接:专栏传送门--网路安全技术创作初心:本博客的初心为与技术朋友们相互交流,每个人的技术都存在短板,博主也是一样,... 查看详情

crackme练习——acid_burn(代码片段)

前言首先十分感谢整理这些CrackMe的吾爱大佬,真是方便了我这样的小白。希望最后这160个CrackMe都可以做出,我会认真做好分析过程的。准备系统:Windows7SP1x64ultimate工具:PEiD、OD分析运行程序,看下需要做些什么。程序先弹出了... 查看详情