在遗忘边缘试探的matlab函数

深藏阁楼爱情的钟 2023-02-19 01:28 104阅读 0赞

### matlab使用 ###

*  文件
 *  矩阵
 *  结构体/cell?
 *   *  cell
 *  图像
 *  点云

注意：mtlab里面的访问都是从1开始的

# 文件 #

1.  读取txt

load ('center.txt');

1.  批量读取文件

for i=0:4
        str=strcat ('tree3_cluster_',int2str(i),'.ply');%文件名
        pc1{ i+1}=pcread(str);%数组索引必须为正整数或逻辑值，所以只能加1
    end

1.  输出  
    把矩阵保存成txt  
    如果东西不多的话，直接打开工作区，手动选择，然后复制到txt就好啦！（而且不用担心空的表格，它会自动变成空格或者省略的）  
    ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70]  
    ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 1]

# 矩阵 #

B=inv(A);%1. 矩阵求逆

B=det(A);%2. 求行列式

C=A*B;%行乘列 3. 矩阵相乘

C=A.*B;%对应元素相乘4. 矩阵点乘

A(i,:)%提取第i行  5. 选择矩阵的一列或一行
     A(:,i)%提取第i列

sort(y,index); 6. 矩阵排序

A=[];  6. 矩阵清空
    A(1,:)=[];%把第一行清空，后面的行会自动往前补

1.  获得矩阵的列数

size(A,2)；%获得矩阵A的列数
    [r,w]=size(B);%获得矩阵的行和列数

1.  矩阵拼接

C=[A;B];
    C=[A,B];

1.  画图  
    画曲线

t=0:pi/60:2*pi;
    x=5*sin(t);
    y=5*cos(t);
    z=3*t;
    c=5*cos(asin(x/5));%asin是arcsin即sin的反函数
    plot3(x,y,z);%三维曲线
    plot(x,c);%二维曲线
    plot(x,z);%二维曲线

1.  函数求导  
    f=diff(y);%f为函数y的导数
2.  matlab经常出现的误差函数  
    [matlab 万能实用的非线性曲线拟合方法][matlab]  
    ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 2]  
    ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 3]  
    ![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 4]

# 结构体/cell? #

## cell ##

1.  赋值和访问  
    小技巧：不知道怎么访问时，去工作区把这个表格选中复制一下到代码行或者命令行，就自动生成访问代码啦！

for i=1：4
        str=strcat ('tree3_cluster_',int2str(i),'.ply');
        pc1{ i}=pcread(str);%数组索引必须为正整数或逻辑值，所以只能加1
        pc{ i}.Location=pc1{ i}.Location;
        pc{ i}.Color=pc1{ i}.Color;
        pc{ i}.Count=pc1{ i}.Count;
        nc=nc+1;
    end

这里生成了两种cell组，不知道具体怎么叫  
第一组pc1包含4个cell，每一个cell都是pointCloud类型的，这里面包含Location、Color、Normal、Intensity、Count|、XLimits、YLimits、ZLimits这些点云里面的东西  
第二组pc是我自己定义的，方便读写的，也是包含4个cell，但是我只定义了里面包含Location、Color、Count  
![在这里插入图片描述][20210602194449162.png]  
![在这里插入图片描述][2021060219450380.png]  
访问：

pc{ 1, 1}.Location=A;

1.  把cell的一个组清空，其他组前移一位

pc(:,1)=[];%如果是1行n列（即n个cell）的组，要清空第一个，并且后面的移上来，要这样写。pc(1,1)=[];这样写就出错，不太懂为啥
    pc{ :,1}=[];%这样只能把里面的结构每一个都清空，但是这个cell还存在，只不过里面没值

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 5]

# 图像 #

1.  保存图片  
    imwrite(‘图片.jpg’,g);%g是图像矩阵

# 点云 #

1.  读取点云  
    tips:如果有点云的txt文件那可以直接读取，方便快捷

pc = pcread('Road.ply'); %读入点云文件

注意：读取的pointCloud格式的数据（此处为pc），只读，不能写，要导出来，在别的地方来写操作

1.  显示点云

pcshow(pc0);%使用存储在点云对象中的位置和颜色显示点。

1.  定义新点云

ptCloud=pointCloud(A);%把矩阵A保存为pointCloud格式 ，自动把A保存到Location，自动计算count

1.  保存点云

pcwrite(ptCloud,'pointcloud.ply');

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20210602200630957.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 1]: https://img-blog.csdnimg.cn/20210602200708456.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70
[matlab]: https://blog.csdn.net/ljyljyok/article/details/81624496?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522159704792419725211949239%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall.%2522%257D&request_id=159704792419725211949239&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~all~first_rank_v2~rank_v28-6-81624496.pc_first_rank_v2_rank_v28&utm_term=matlab+%E6%8B%9F%E5%90%88%E6%9B%B2%E7%BA%BF%E8%BE%93%E5%87%BA%E6%96%B9%E7%A8%8B&spm=1018.2118.3001.4187
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 2]: https://img-blog.csdnimg.cn/20200810212848944.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 3]: https://img-blog.csdnimg.cn/20200810212940666.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 4]: https://img-blog.csdnimg.cn/20200810212955859.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70
[20210602194449162.png]: https://img-blog.csdnimg.cn/20210602194449162.png
[2021060219450380.png]: https://img-blog.csdnimg.cn/2021060219450380.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5_size_16_color_FFFFFF_t_70 5]: https://img-blog.csdnimg.cn/20210601212421366.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQwNzkxMDk5,size_16,color_FFFFFF,t_70