MATLAB指令

最新推荐文章于 2024-07-09 15:39:11 发布

山人自有锦囊妙计

最新推荐文章于 2024-07-09 15:39:11 发布

阅读量2.4k

点赞数 3

分类专栏： MATLAB

本文链接：https://blog.csdn.net/qq_20386411/article/details/83542299

版权

MATLAB 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

matlab207的安装：

安装方法：https://blog.csdn.net/gisboygogogo/article/details/76793803

当时复制了此文件就可以打开安装的matlab，没有任何错误提示框弹出。

当一个工程文件中的文件过多，移植程序时，先看清楚程序中调用的路径是否能对上，对上后再运行程序

路径：

1、选择工程文件所在的目录（不同.m文件放在同一目录下）

2、关于路径的问题，好像只能将.m文件放在一个文件夹中，在文件夹中再放文件夹（调用函数）总是提示某个函数找不到

全部放在一个文件夹中，就没有问题了

数据获取方法：

1、load直接加载txt文本数据

用fopen函数不能读出数据，只能打开文件，取出数据必须用load函数。若不能打开指定文件，fopen函数返回-1

2、加速度数据

取第一列时间数据和x轴加速度数据

3、读取csv格式的数据，直接使用csvread函数

filename='/Users/Steven Xiao/Desktop/wifi感知/X4M03开发板/datum.csv';
aa = csvread(filename);
plot(aa(1,:))

4、保存变量

csvwrite('total_wave.csv',total_wave);，第一个参数是文件名和文件格式，第二个参数是变量矩阵，保存在运行路径下

下标索引必须为正整数

matlab下标索引必须为正整数类型或逻辑类型。取范围时只能从正数开始取

数据保存

1，将所需变量保到txt文本中

filename='F:\指夹、手环、Biopac实验数据（志愿者）\E4数据\李志成_1532571159_A01671\BVP.csv';%2018/7/26 10:12:39
bvp=load(filename);

bvp1=bvp(19200:23040);
fop=fopen('2_25.txt','w');
for i=1:(length(bvp1))
    a=bvp1(i,1);

    fprintf(fop,'%g ',bvp1(i,1));
    fprintf(fop,'\n');

end
fclose(fop);

①虽然结果中不是一行行显示，但是已经进行了换行处理

②fclose(fop);必须在end后

基本操作

1、close all % 可以一次关闭所有的figure

调试方法：

断点直接使用鼠标打在语句的最前面

(2) 执行：F5。

(3) 单步执行：F10。

(4) step in：F11。

(5) step out：Shift+F11

时间戳(单位为秒)转化为时间方法：

1、直接使用在线工具转化http://tool.chinaz.com/Tools/unixtime.aspx

同样把时间转化为时间戳

常见错误:

1、关于找不到.m文件的提示

在一个.m文件中调用另一个.m文件时，直接用.m的文件名作为函数，不用.m中的函数名function后面的内容。所以函数名最好和.m文件的名称一样

main中的函数名不能写成ideal_bandpassing，否则ctr+D时就会出现下面错误

画图：

https://blog.csdn.net/jiangnankid/article/details/22690165

1、

①画茎状图

M=4; %原离散信号有4点
n=[0:1:M-1]; %原信号是1行4列的矩阵
xn=[1 1 1 1]; %构建原始信号
subplot(3,1,1);
stem(n,xn); %画图
title('原始信号');

画柱状图

color_set=[0.8660 0.9236 0.8957 0.9116];
%bar的颜色索引
color_background=['c' 'm' 'y' 'k' 'r' 'g' 'b'];
%对figure的标题、横坐标、纵坐标的标注进行设置
hold on
%title('Tiltle Name');
ylabel('sVRI');
ax = gca;
ax.XTick = [1 2 3 4];
ax.XTickLabels = {'Rest','Hard Task','Relax','Esay Task'};
%对不同的数据所对应的bar进行不同颜色条的设计
color_bar=bar(1,0.8660);
set(color_bar,'FaceColor',color_background(1));%选择‘c’这个颜色
color_bar=bar(2,0.9236);
set(color_bar,'FaceColor',color_background(2));%选择‘m’这个颜色
color_bar=bar(3,0.8957);
set(color_bar,'FaceColor',color_background(3));%选择‘y’这个颜色
color_bar=bar(4,0.9116);
set(color_bar,'FaceColor',color_background(4));%选择‘y’这个颜色
axis([-inf,inf,0.7,1.0]) 
%画柱状图
data = [0.8660 ,0.9236 ,0.8957,0.9116];
b = bar(data);
color_background=['c' 'm' 'y' 'k' 'r' 'g' 'b'];
ch = get(b,'children');
set(gca,'XTickLabel',{'Rest','Hard Task','Relax','Esay Task'})
set(ch,'FaceVertexCData',[0 0 1;0 1 1;1 1 1;1 0 1;])
ylabel('sVRI');
axis([-inf,inf,0.7,1.0])

在柱状图上加高度

color_set11=[0.8660];
color_set1=num2cell(color_set11);

color_bar=bar(1,0.8660);
set(color_bar,'FaceColor',color_background(1));%选择‘c’这个颜色
text(1-0.1,0.8660+0.018,color_set1);

多条曲线加上标注

figure
plot(f02)

hold on
plot(s01)

hold on
plot(t12)
legend('0-2','0-1','1-2')

一个legend即可

②曲线加粗

plot(x,y,'Linewidth',3)，可以将x，y直接换成矩阵数据

③曲线颜色plot(横轴，纵轴，‘形状颜色’)，形状和颜色连在一起

如‘*r’表示离散的*点，所有点之间没有连接起来

‘-*r’表示将所有*用直线连接起来

*/○的大小：plot(breathfreq,'-ko','markersize',4)

b 蓝 .
c 青。圈
g 绿 × ×标记
k 黑－实线
m 紫红 * 星号
r 红：点线
w 白－. 点划线
y 黄－－虚线

④figure也可以传入参数，figure('visible','off')，只是不画出来，但是实际画图的时候会占用figure

2、横轴点数转化为时间轴

①将矩阵中的横轴（数字轴）转化为对应的时间轴

len=length(se);
t=0.005*(1:len);%很轴分成多分
plot(t,second_ant_csi.')%指定横轴为时间，数据中的横轴含义直接被覆盖了

注：指定横轴，纵坐标可以直接传入数据，这里的一行数据要直接转化成一列数据

②只指定某一个坐标轴的范围ylim([0 0.6])/xlim([0 0.6])

使用 axis([XMIN XMAX YMIN YMAX]) ,对于不需要设置范围的轴，使用参数-inf,inf，即正负无穷之间

>> x = 1:10;
>> y = x.^2;
>> subplot(3,1,1)
>> plot(x,y,'r')
>> subplot(3,1,2)
>> plot(x,y,'g')
>> xlim([0,15])
>> subplot(3,1,3)
>> plot(x,y,'b')
>> axis([-1,12,-inf,inf])

必须在plot后面

3、多图和曲线

①怎样实现一个界面画多个图

    subplot(2,2,1);% 新建窗口 
    plot(I)             % 画出原始的波形

    e=I;
    subplot(2,2,3);                                      % 新建窗口
    s1=wden(e,'minimaxi','s','one',5,'db3');   % 选用 db3 小波对信号进行5层分解，并对细节系数选用 minimaxi 阈值模式和尺度噪声 
    plot(s1);

直接使用subplot()函数,前两个参数确定将界面分成多少块，最后一个参数决定图画在那一块

②在一个坐标系中画多个曲线hold on...hold off

    subplot(2,2,1)
    plot(csi15_1)
    hold on 
    plot(csi15_2)
    plot(csi15_3)
    hold off

注：

hold on 和hold off，是相对使用的.通常是一个图上画两个曲线进行比较。

% 前者的意思是，你在当前图的轴（坐标系）中画了一幅图，再画另一幅图时，原来的图还在，与新图共存，都看得到

% 后者表达的是，你在当前图的轴（坐标系）中画了一幅图，此时，状态是hold off,则再画另一幅图时，原来的图就看不到了，在轴上绘制的是新图，原图被替换了

③给最大值点位置加上坐标线

plot(f_PSD,plot_Pxx1);
j = f_PSD>1 & f_PSD<2 ; %指定频率范围求最大值和最小值
max_psd=max(plot_Pxx1(j)); %求出指定横轴范围内的纵轴的最大值
ax=axis;%获得当前坐标的范围
hold on;%保持图像
plot([f_PSD(find(plot_Pxx1==max_psd)),f_PSD(find(plot_Pxx1==max_psd))],[ax(3),ax(4)],'r:');%这里只是绘制了竖线

注：

1、先画出图像后再画纵横线

2、plot([ax(1),ax(2)],[maxy,maxy],'r:',[maxx,maxx],[ax(3),ax(4)],'r:');%绘制纵横线

5、给每一条曲线加上注释