首页 > 其他 > 详细

MATLAB数据处理快速学习教程

时间:2014-05-13 19:40:27      阅读:546      评论:0      收藏:0      [点我收藏+]

转自:http://blog.csdn.net/abcjennifer/article/details/7706581

本篇内容集合了MATLAB中的基本操作、数据存储与计算、数据的直线与曲线拟合与画图、for-if-while语句的使用方法对一般matlab的操作进行了基本和详细的应用举例,非常适合初学者进行matlab学习快速掌握。下面分四个模块分别进行讲述:


========================BasicOperations========================

>> A=rand(3,2)

 

A =

 

   0.8147    0.9134

   0.9058    0.6324

   0.1270    0.0975

 

>> A=[1 2;3 4;5 6]

 

A =

 

    1     2

    3     4

    5     6

 

>> A(3,2)

 

ans =

 

    6

 

>> A(2,:)

 

ans =

 

    3     4

 

>> A([1 3],:) %select the 1st and 3rd row

 

ans =

 

    1     2

    5     6

 

>> save test.mat A

>> save testtext.txt A –ascii

 

bubuko.com,布布扣


>> A(:,2)=[10,11,12]

 

A =

 

    1    10

    3    11

    5    12

 

>> A=[A,[101;102;103]]%append another column vector

 

A =

 

    1    10   101

    3    11   102

    5    12   103

 

>> A(:)%put all elements of A into a single vector

 

ans =

 

    1

    3

    5

   10

   11

   12

  101

  102

  103

 

B=[11,12;13,14;15,16]

 

B =

 

   11    12

   13    14

   15    16

 

>> C=[A B]

 

C =

 

    1    10   101   11    12

    3    11   102   13    14

    5    12   103   15    16

 

A=[1 2;3 4;5 6]

 

A =

 

    1     2

    3     4

    5     6

 

>> C=[A;B]

 

C =

 

    1     2

    3     4

    5     6

   11    12

   13    14

15    16

 

========================Computingon Data========================

 

>>A.*B%对应位置的数据相乘,element product

 

ans =

 

   11    24

   39    56

75    96

 

>>A‘%转置

 

ans =

 

    1     3     5

    2     4     6

 

>> A<3

 

ans =

 

    1     1

    0     0

    0     0

 

>> find(A<3)

 

ans =

 

    1

    4

 

A =

 

    8     1     6

    3     5     7

    4     9     2

 

>> [r,c]=find(A>=6)

 

r =

 

    1

    3

    1

    2

 

 

c =

 

    1

    2

    3

    3

 

a=[1 15 2 0.5]

 

a =

 

   1.0000   15.0000    2.0000   0.5000

 

>> sum(a)

 

ans =

 

  18.5000

 

>> prod(a)

 

ans =

 

   15

 

>> floor(a)%取下界

 

ans =

 

    1    15     2    0

 

>> ceil(a)%取上界

 

ans =

 

    1    15     2    1

 

rand(3)%创建3*3random矩阵,每个值在[0,1]之间

 

ans =

 

   0.6463    0.2760    0.1626

   0.7094    0.6797    0.1190

   0.7547    0.6551    0.4984

 

>> max(rand(3),rand(3)) %在两个random3*3矩阵中找对应位置的max

 

ans =

 

   0.9597    0.2238    0.5060

   0.5472    0.7513    0.8143

0.5853    0.8407   0.8909

 

A=magic(3)

 

A =

 

    8     1     6

    3     5     7

    4     9     2

 

>> max(A,[],1) %找每列最大值,1表示第一维,即列

 

ans =

 

    8     9     7

 

>> max(A,[],2) %找每行最大值,2表示第二维,即行

 

ans =

 

    8

    7

    9

 

max(A) %defaultis column max

 

ans =

 

    8     9     7

 

>> max(max(A))

 

ans =

 

    9

 

>> A(:)

 

ans =

 

    8

    3

    4

    1

    5

    9

    6

    7

    2

 

>> max(A(:))

 

ans =

 

    9

 

>>

>>

>>

>> A=magic(9)

 

A =

 

   47    58    69   80     1    12   23    34    45

   57    68    79    9    11    22   33    44    46

   67    78     8   10    21    32   43    54    56

   77     7    18   20    31    42   53    55    66

    6    17    19   30    41    52   63    65    76

   16    27    29   40    51    62   64    75     5

   26    28    39   50    61    72   74     4    15

   36    38    49   60    71    73    3    14    25

   37    48    59   70    81     2   13    24    35

 

>> sum(A,1)%columnsum

 

ans =

 

  369   369   369  369   369   369  369   369   369

 

>> sum(A,2)%sumeach row

 

ans =

 

  369

  369

  369

  369

  369

  369

  369

  369

  369

 

>> eye(9)

 

ans =

 

    1     0     0    0     0     0    0     0     0

    0     1     0    0     0     0    0     0     0

    0     0     1    0     0     0    0     0     0

    0     0     0    1     0     0    0     0     0

    0     0     0    0     1     0    0     0     0

    0     0     0    0     0     1    0     0     0

    0     0     0    0     0     0    1     0     0

    0     0     0    0     0     0    0     1     0

    0     0     0    0     0     0    0     0     1

 

>> A.*eye(9)%takethe element product of the 2 matrix

 

ans =

 

   47     0     0    0     0     0    0     0     0

    0    68     0    0     0     0    0     0     0

    0     0     8    0     0     0    0     0     0

    0     0    0    20     0    0     0     0    0

    0     0     0    0    41     0    0     0     0

    0     0     0    0     0    62    0     0     0

    0     0     0    0     0     0   74     0     0

    0     0     0    0     0     0    0    14     0

    0     0     0    0     0     0    0     0    35

 

>> sum(sum(A.*eye(9)))

 

ans =

 

  369

 

>> A=magic(3)

 

A =

 

    8     1     6

    3     5     7

    4     9     2

 

>> temp=pinv(A) %矩阵求逆

 

temp =

 

   0.1472   -0.1444    0.0639

  -0.0611    0.0222    0.1056

  -0.0194    0.1889   -0.1028

 

>> temp*A

 

ans =

 

   1.0000   -0.0000   -0.0000

  -0.0000    1.0000    0.0000

   0.0000    0.0000    1.0000

 

========================PlottingData========================

 

>>t=[0:0.01:0.98];

>> y1=sin(2*pi*4*t);

>> plot(t,y1)

bubuko.com,布布扣

>> hold on;%plotnew figure on the old ones

>> y2=cos(2*pi*4*t);

>> plot(t,y2,‘r‘)

 bubuko.com,布布扣

>> xlabel(‘time‘)
>> ylabel(‘value‘)
>> legend(‘sin‘,‘cos‘)
>> title(‘my plot‘)
>> print -dpng ‘myplot.png‘%save as a file in default catalog

>> close

bubuko.com,布布扣


%分别显示两幅图像

>> figure(1);plot(t,y1);

>> figure(2);plot(t,y2)
>> 

%一幅图中显示两个subplot

 figure
subplot(1,2,1);plot(t,y1)
subplot(1,2,2);plot(t,y2)
axis([0.5 1 -1 1])

bubuko.com,布布扣

 

  A=magic(5)

imagesc(A)

bubuko.com,布布扣


>>imagesc(A),colorbar,colormap gray;

 bubuko.com,布布扣

 

 

 ========================for if while statements========================
  1. > v=zeros(10,1)  
  2. for i=1:10  
  3. v(i)=2^i;  
  4. end  
  5. v  

v =


     0
     0
     0
     0
     0
     0
     0
     0
     0
     0




v =


           2
           4
           8
          16
          32
          64
         128
         256
         512
        1024

  1. indices=1:10;  
  2. indices  
  3. for i=indices, disp(i);  
  4. end  
  5. i=1;  
  6. while i<=5,  
  7. v(i)=100;  
  8. i=i+1;  
  9. end;  
  10. v  
  11. i=1;  
  12. while true,  
  13. v(i)=999;  
  14. i=i+1;  
  15. if i==6,  
  16. break;  
  17. end  
  18. end  
  19. v  
  20. v(1)=2;  
  21. if v(1)==1,  
  22. disp(‘The value is one‘);  
  23. else if v(1)==2,  
  24. disp(‘The value is two‘);  
  25. else  
  26. disp(‘The value is neither one or two‘);  
  27. end  
  28. end  


另外还有一些matlab高级应用:

http://emuch.net/bbs/viewthread.php?tid=4340045

MATLAB数据处理快速学习教程,布布扣,bubuko.com

MATLAB数据处理快速学习教程

原文:http://www.cnblogs.com/zhizhan/p/3724229.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!