matlab中speye,Matlab相关操作

最新推荐文章于 2023-08-12 21:16:55 发布

HR-Henry陈恒

最新推荐文章于 2023-08-12 21:16:55 发布

阅读量1.1k

点赞数

文章标签： matlab中speye

>> a=3

a =

>> a=3; //加分号就暂时不输出a

>> a

a =

>> disp(a) //disp函数就是display a的意思，显示出a的数值

>> disp(sprintf('2 decimals:%0.2f',a)) //sprintf就是可以将参数转换成字符串

2 decimals:3.00

>> a=pi;

>> a

a =

3.141592653589793

>> format long //可以调整格式

>> a

a =

3.141592653589793

>> format short

>> a

a =

3.1416

>> A=[1,2,3] //逗号分隔即为一列

A =

1 2 3

>> [1,2,3;2,3,4] //分号分隔为一行

ans =

1 2 3

2 3 4

>> v=1:0.2:2 //意思为从1开始，0.2为单位，跳到2

v =

1 至 4 列

1.0000 1.2000 1.4000 1.6000

5 至 6 列

1.8000 2.0000

>> v=1:6 //从1跳到6，默认1为单位

v =

1 2 3 4 5 6

>> ones(2,3) //生成2行3列的全为1的矩阵

ans =

1 1 1

>> ans=2*ones(2,3)

ans =

2 2 2

>> w=zeros(1,3) //1行3列全为0的矩阵

w =

0 0 0

>> w=rand(1,3) //1行3列的随机，估计是【0，1)

w =

0.8147 0.9058 0.1270

>> rand(3,3)

ans =

0.9134 0.2785 0.9649

0.6324 0.5469 0.1576

0.0975 0.9575 0.9706

>> rand(3,3)

ans =

0.9572 0.1419 0.7922

0.4854 0.4218 0.9595

0.8003 0.9157 0.6557

>> randn(1,3) //更加高级的随机函数

ans =

-1.2075 0.7172 1.6302

>> hist(ans) //生成直方图

>> eye(4) //生成单位矩阵

ans =

1 0 0 0

0 1 0 0

0 0 1 0

0 0 0 1

>> help eye // help加函数名可以了解更多

eye - 单位矩阵

此 MATLAB 函数返回标量 1。

I = eye

I = eye(n)

I = eye(n,m)

I = eye(sz)

I = eye(classname)

I = eye(n,classname)

I = eye(n,m,classname)

I = eye(sz,classname)

I = eye('like',p)

I = eye(n,'like',p)

I = eye(n,m,'like',p)

I = eye(sz,'like',p)

另请参阅 ones, speye, zeros

eye 的参考页

名为 eye 的其他函数

>> eye(4,2)

ans =

1 0

0 1

0 0

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

A =

1 2

3 4

5 6

>> size(A)

ans =

3 2

>> sz=size(A)

sz =

3 2

>> size(sz)

ans =

1 2

>> size(A,1) //A矩阵的行数

ans =

>> size(A,2) //A矩阵的列数

ans =

>> v=[1 2 3 4]

v =

1 2 3 4

>> length(v) //对只有一行的矩阵给出长度

ans =

>> length(A) //对矩阵给出行数

ans =

>> length([1;2;3;4;5]) //对向量给出长度

ans =

>> pwd //现在所在位置

ans =

'/Users/chenzeyuan/Documents/MATLAB'

>> who //查询当前的变量

您的变量为:

A ans sz v

>> whos //详细查询变量

Name Size Bytes Class Attributes

A 3x2 48 double

ans 1x34 68 char

sz 1x2 16 double

v 1x4 32 double

>> clear ans //删去ans这个变量

>> A(3,2) //A3行2列的数

ans =

>> A(2,:) %means every element along that row/column

ans =

3 4

>> A(:,2)

ans =

>> A([1 3],:) //给出第1行和第3行所有列的矩阵

ans =

1 2

5 6

>> A(:,2)=[10;11;12] //把10，11，12赋给A的第2列

A =

1 10

3 11

5 12

>> A=[A,[100;101;102]] //加上100，101，102一列

A =

1 10 100

3 11 101

5 12 102

>> A(:) //把矩阵变向量

ans =

100

101

102

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

B =

11 12

13 14

15 16

>> C=[A B] //矩阵A和矩阵B合起来

C =

1 10 100 11 12

3 11 101 13 14

5 12 102 15 16

>> B=[B,[13;14;15]]

B =

11 12 13

13 14 14

15 16 15

>> C=[A;B] //把矩阵B加在矩阵A下

C =

1 10 100

3 11 101

5 12 102

11 12 13

13 14 14

15 16 15

>> size(C)

ans =

6 3

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

A =

1 2

3 4

5 6

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

B =

11 12

13 14

15 16

>> C=[1 1;2 2]

C =

1 1

2 2

>> A*C //进行矩阵的运算

ans =

5 5

11 11

17 17

>> A .^ 2 //A矩阵所有数乘以2

ans =

1 4

9 16

25 36

>> log(A) //矩阵所有数求log

ans =

0 0.6931

1.0986 1.3863

1.6094 1.7918

>> exp(A) //矩阵所有数进行e的幂运算，次数为矩阵对应的值

ans =

2.7183 7.3891

20.0855 54.5982

148.4132 403.4288

>> abs(A) //求绝对值

ans =

1 2

3 4

5 6

>> ones(3,1) //3行1列全为1

ans =

>> ones(3,2)

ans =

1 1

>> A+ones(3,2)

ans =

2 3

4 5

6 7

>> A+1 //A+1==A+ones(3,2)，都是对所有数加一

ans =

2 3

4 5

6 7

>> A' //转置矩阵

ans =

1 3 5

2 4 6

>> (A')'

ans =

1 2

3 4

5 6

>> a=[1 15 2 0.5]

a =

1.0000 15.0000 2.0000 0.5000

>> val=max(a) //求最大值

val =

>> [val,ind]=max(a) //前面为最大值，后为其下标

val =

ind =

>> max(A) //对于一个矩阵来说，就是对每列求最大值

ans =

5 6

>> a<3 //比3小的为1，大的为0

ans =

1×4 logical 数组

1 0 1 1

>> find(a<3) //找到比3小的数的下标

ans =

1 3 4

>> A=magic(3) //生成3X3的矩阵，行与列加起来都一样

A =

8 1 6

3 5 7

4 9 2

>> [r,c]=find(A>=7) //对于矩阵的find，r表示行，c表示列

r =

c =

>> a

a =

1.0000 15.0000 2.0000 0.5000

>> sum(a) //计算总和，都加起来

ans =

18.5000

>> prod(a) //全部相乘

ans =

>> sum(A) //对于矩阵，对每列求和

ans =

15 15 15

>> prod(A) //对于矩阵，每列求乘积

ans =

96 45 84

>> floor(a) //向下求整

ans =

1 15 2 0

>> ceil(a) //向上求整

ans =

1 15 2 1

>> rand(3) //生成随机3X3矩阵

ans =

0.8147 0.9134 0.2785

0.9058 0.6324 0.5469

0.1270 0.0975 0.9575

>> max(A)

ans =

8 9 7

>> b=rand(3)

b =

0.9649 0.9572 0.1419

0.1576 0.4854 0.4218

0.9706 0.8003 0.9157

>> c=rand(3)

c =

0.7922 0.0357 0.6787

0.9595 0.8491 0.7577

0.6557 0.9340 0.7431

>> max(b,c) //求b，c矩阵对应位置的最大值

ans =

0.9649 0.9572 0.6787

0.9595 0.8491 0.7577

0.9706 0.9340 0.9157

>> max(A,[],1) //对于列里找最大

ans =

8 9 7

>> max(A,[],2) //对于行里找最大

ans =

>> max(max(A))

ans =

>> A(:)

ans =

>> max(A(:))

ans =

>> A=magic(9)

A =

1 至 7 列

47 58 69 80 1 12 23

57 68 79 9 11 22 33

67 78 8 10 21 32 43

77 7 18 20 31 42 53

6 17 19 30 41 52 63

16 27 29 40 51 62 64

26 28 39 50 61 72 74

36 38 49 60 71 73 3

37 48 59 70 81 2 13

8 至 9 列

34 45

44 46

54 56

55 66

65 76

75 5

4 15

14 25

24 35

>> sum(A,1) //对每列求和

ans =

1 至 7 列

369 369 369 369 369 369 369

8 至 9 列

369 369

>> sum(A,2) //对每行求和

ans =

369

>> eye(9) //生成单位矩阵

ans =

1 至 7 列

1 0 0 0 0 0 0

0 1 0 0 0 0 0

0 0 1 0 0 0 0

0 0 0 1 0 0 0

0 0 0 0 1 0 0

0 0 0 0 0 1 0

0 0 0 0 0 0 1

0 0 0 0 0 0 0

8 至 9 列

0 0

1 0

0 1

>> A .* eye(9)

ans =

1 至 7 列

47 0 0 0 0 0 0

0 68 0 0 0 0 0

0 0 8 0 0 0 0

0 0 0 20 0 0 0

0 0 0 0 41 0 0

0 0 0 0 0 62 0

0 0 0 0 0 0 74

0 0 0 0 0 0 0

8 至 9 列

0 0

14 0

0 35

>> flipud(eye(9)) //换个对角线

ans =

1 至 7 列

0 0 0 0 0 0 0

0 0 0 0 0 0 1

0 0 0 0 0 1 0

0 0 0 0 1 0 0

0 0 0 1 0 0 0

0 0 1 0 0 0 0

0 1 0 0 0 0 0

1 0 0 0 0 0 0

8 至 9 列

0 1

1 0

0 0

>> pinv(A) //自反矩阵

ans =

1 至 4 列

0.0005 -0.0012 0.0017 0.0126

0.0003 0.0003 0.0128 -0.0122

-0.0011 0.0140 -0.0120 0.0003

0.0126 -0.0120 0.0003 0.0003

-0.0108 0.0003 0.0003 0.0003

0.0003 0.0005 0.0001 -0.0011

0.0003 0.0003 0.0003 0.0005

0.0003 0.0003 -0.0011 0.0017

0.0003 0.0003 0.0003 0.0003

5 至 8 列

-0.0121 0.0003 0.0003 0.0003

0.0003 -0.0011 0.0017 0.0003

0.0005 0.0001 0.0003 0.0003

-0.0011 0.0017 0.0005 0.0001

0.0003 0.0003 0.0003 0.0003

0.0017 0.0003 0.0003 0.0126

0.0001 0.0003 0.0126 -0.0134

0.0003 0.0128 -0.0122 0.0003

0.0127 -0.0120 -0.0011 0.0018

9 列

0.0003

0.0114

-0.0120

0.0017

0.0003

0.0001

>> A=magic(3)

A =

8 1 6

3 5 7

4 9 2

>> pinv(A)

ans =

0.1472 -0.1444 0.0639

-0.0611 0.0222 0.1056

-0.0194 0.1889 -0.1028

>> pinv(A)*A

ans =

1.0000 0.0000 -0.0000

-0.0000 1.0000 0.0000

0.0000 -0.0000 1.0000

>> t=[0:0.01:0.98]

t =

1 至 4 列

0 0.0100 0.0200 0.0300

5 至 8 列

0.0400 0.0500 0.0600 0.0700

9 至 12 列

0.0800 0.0900 0.1000 0.1100

13 至 16 列

0.1200 0.1300 0.1400 0.1500

17 至 20 列

0.1600 0.1700 0.1800 0.1900

21 至 24 列

0.2000 0.2100 0.2200 0.2300

25 至 28 列

0.2400 0.2500 0.2600 0.2700

29 至 32 列

0.2800 0.2900 0.3000 0.3100

33 至 36 列

0.3200 0.3300 0.3400 0.3500

37 至 40 列

0.3600 0.3700 0.3800 0.3900

41 至 44 列

0.4000 0.4100 0.4200 0.4300

45 至 48 列

0.4400 0.4500 0.4600 0.4700

49 至 52 列

0.4800 0.4900 0.5000 0.5100

53 至 56 列

0.5200 0.5300 0.5400 0.5500

57 至 60 列

0.5600 0.5700 0.5800 0.5900

61 至 64 列

0.6000 0.6100 0.6200 0.6300

65 至 68 列

0.6400 0.6500 0.6600 0.6700

69 至 72 列

0.6800 0.6900 0.7000 0.7100

73 至 76 列

0.7200 0.7300 0.7400 0.7500

77 至 80 列

0.7600 0.7700 0.7800 0.7900

81 至 84 列

0.8000 0.8100 0.8200 0.8300

85 至 88 列

0.8400 0.8500 0.8600 0.8700

89 至 92 列

0.8800 0.8900 0.9000 0.9100

93 至 96 列

0.9200 0.9300 0.9400 0.9500

97 至 99 列

0.9600 0.9700 0.9800

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

y1 =

1 至 4 列

0 0.2487 0.4818 0.6845

5 至 8 列

0.8443 0.9511 0.9980 0.9823

9 至 12 列

0.9048 0.7705 0.5878 0.3681

13 至 16 列

0.1253 -0.1253 -0.3681 -0.5878

17 至 20 列

-0.7705 -0.9048 -0.9823 -0.9980

21 至 24 列

-0.9511 -0.8443 -0.6845 -0.4818

25 至 28 列

-0.2487 -0.0000 0.2487 0.4818

29 至 32 列

0.6845 0.8443 0.9511 0.9980

33 至 36 列

0.9823 0.9048 0.7705 0.5878

37 至 40 列

0.3681 0.1253 -0.1253 -0.3681

41 至 44 列

-0.5878 -0.7705 -0.9048 -0.9823

45 至 48 列

-0.9980 -0.9511 -0.8443 -0.6845

49 至 52 列

-0.4818 -0.2487 -0.0000 0.2487

53 至 56 列

0.4818 0.6845 0.8443 0.9511

57 至 60 列

0.9980 0.9823 0.9048 0.7705

61 至 64 列

0.5878 0.3681 0.1253 -0.1253

65 至 68 列

-0.3681 -0.5878 -0.7705 -0.9048

69 至 72 列

-0.9823 -0.9980 -0.9511 -0.8443

73 至 76 列

-0.6845 -0.4818 -0.2487 -0.0000

77 至 80 列

0.2487 0.4818 0.6845 0.8443

81 至 84 列

0.9511 0.9980 0.9823 0.9048

85 至 88 列

0.7705 0.5878 0.3681 0.1253

89 至 92 列

-0.1253 -0.3681 -0.5878 -0.7705

93 至 96 列

-0.9048 -0.9823 -0.9980 -0.9511

97 至 99 列

-0.8443 -0.6845 -0.4818

>> plot(t,y1) //生成一个图

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

y2 =

1 至 4 列

1.0000 0.9686 0.8763 0.7290

5 至 8 列

0.5358 0.3090 0.0628 -0.1874

9 至 12 列

-0.4258 -0.6374 -0.8090 -0.9298

13 至 16 列

-0.9921 -0.9921 -0.9298 -0.8090

17 至 20 列

-0.6374 -0.4258 -0.1874 0.0628

21 至 24 列

0.3090 0.5358 0.7290 0.8763

25 至 28 列

0.9686 1.0000 0.9686 0.8763

29 至 32 列

0.7290 0.5358 0.3090 0.0628

33 至 36 列

-0.1874 -0.4258 -0.6374 -0.8090

37 至 40 列

-0.9298 -0.9921 -0.9921 -0.9298

41 至 44 列

-0.8090 -0.6374 -0.4258 -0.1874

45 至 48 列

0.0628 0.3090 0.5358 0.7290

49 至 52 列

0.8763 0.9686 1.0000 0.9686

53 至 56 列

0.8763 0.7290 0.5358 0.3090

57 至 60 列

0.0628 -0.1874 -0.4258 -0.6374

61 至 64 列

-0.8090 -0.9298 -0.9921 -0.9921

65 至 68 列

-0.9298 -0.8090 -0.6374 -0.4258

69 至 72 列

-0.1874 0.0628 0.3090 0.5358

73 至 76 列

0.7290 0.8763 0.9686 1.0000

77 至 80 列

0.9686 0.8763 0.7290 0.5358

81 至 84 列

0.3090 0.0628 -0.1874 -0.4258

85 至 88 列

-0.6374 -0.8090 -0.9298 -0.9921

89 至 92 列

-0.9921 -0.9298 -0.8090 -0.6374

93 至 96 列

-0.4258 -0.1874 0.0628 0.3090

97 至 99 列

0.5358 0.7290 0.8763

>> plot(t,y1);

>> hold on; //让上面生成的图不消失

>> plot(t, y2,'r') //这里生成的图加到上面去

>> xlabel('time') //x轴名字

>> ylabel('value') //y轴名字

>> legend('sin','cos') //提示sin是什么颜色，cos是什么颜色

>> title('my plot') //图的标题

>> print -dpng 'myplot.png' //生成图

>> close //关闭图

>> figure(1);plot(t,y1); //给每个图命名，同时不会被下面生成的图取代

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

>> close

>> subplot(1,2,1) //生成1X2的图，可进入第一个

>> plot(t,y1) //在第一个图里生成

>> subplot(1,2,2) //生成上述图，可进入第二个

>> plot(t,y2) //在第二个图里生成

>> close

>> subplot(1,2,1)

plot(t,y1)

subplot(1,2,2)

plot(t,y2)

>> axis([0.5 1,-1 1]) //前两个参数控制x的范围，后两个控制y的范围

>> clf; //清除图

>> A=magic(5)

A =

17 24 1 8 15

23 5 7 14 16

4 6 13 20 22

10 12 19 21 3

11 18 25 2 9

>> imagesc(A) //有眼色标记的矩阵的图

>> imagesc(A),colorbar,colormap gray //调出colorbar同时将颜色调为灰色

>> v=zeros(10,1)

v =

>> for i=1:10,

v(i)=2^i;

end;

>> v

v =

128

256

512

1024

>> indices=1:10;

>> indices

indices =

1 至 7 列

1 2 3 4 5 6 7

8 至 10 列

8 9 10

>> for i=indices,

disp(i);

end;

>> v

v =

128

256

512

1024

>> i=1;

>> while i<=5,

v(i)=100;

i=i+1;

end;

>> v

v =

100

128

256

512

1024

>> v(1)=2;

>> if v(1)==1,

disp(v(1));

elseif v(1)==2,

disp('ufck pff');

else,

disp('co');

end;

ufck pff

>> hh(3) //在搜索路径里放置.m后缀的文件，里面生成函数，即可调用

ans =

>> A=magic(3)

A =

8 1 6

3 5 7

4 9 2

>> x=[1;2;3]

x =

>> A .* x

ans =

8 1 6

6 10 14

12 27 6

>> A*x

ans =

>> sum(A*x)

ans =

>> A(:,1:2)

ans =

8 1

3 5

4 9

>> A(1:3,1:2)

ans =

8 1

3 5

4 9

>> A/4

ans =

2.0000 0.2500 1.5000

0.7500 1.2500 1.7500

1.0000 2.2500 0.5000

>> A ./ 4

ans =

2.0000 0.2500 1.5000

0.7500 1.2500 1.7500

1.0000 2.2500 0.5000

>> log(A)

ans =

2.0794 0 1.7918

1.0986 1.6094 1.9459

1.3863 2.1972 0.6931

>> A .^ 2

ans =

64 1 36

9 25 49

16 81 4

HR-Henry陈恒

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫