數值計算Matlab的一些簡單使用

v=[1;2;3;4]

a=dot(v,v)=30 %對應列相乘，然後相加：1*1+2*2+3*3+4*4=30

a=dot(v,v,1)

b=dot(v,v,2) %對應行相乘，然後相加

e=eig(v) %返回乙個列向量，其中包括方正v的特徵值

(b-a)x=c x=(b-a)\c %ax=b 則x=a'b; xa=b 則x=ba'

c(:,2)=v 第二列換成v

c(2,:)=v 第二行換成v

c(1:10,2)=v 第2列的前10個數換成v

v=c(:,2) 把c的第二列儲存到c

size(c,2) c的列數 size(c,1) 行數

a=rand(1,100) %rand生成乙個0到1之間的隨機數組 1x100

a(1,1) %a矩陣的第一行，第一列的數

求a矩陣裡面大於0.1的數的個數：

for i=1:100

if a(1,i)>0.1
a(1,i)=1
else
a(1,i)=0
endendn=sum(a)

%可以直接用matlab給的函式 factorial(n)

n=input 
f=nwhile n>1
n=n-1
f=f*n
enddisp(['n!=' num2str(f)])
%或者function [y]=fac(x)
if x==0
y=1else
y=1for i=x:-1:1
y=y*i
endend

j=1;

while  j<101
if 0.1<=a(1,j)&&a(1,j)<=0.2
break
else
j=j+1;
endendif j==101
disp('kein element 0.499 <= a_i<= 0.501')
else 
x=['i=',num2str(j),',ai=',num2str(a(1,j))];
disp(x)
end

x=-2*pi:pi/100:2*pi;

y=sin(x);
plot(x,y);
xlabel('x軸名字')
ylabel('y軸名字')
title('a')

function 3d(nodes,elements,sol)

zeile=size(elements,1);x=nodes(:,1);y=nodes(:,2);z=sol;t=zeros(2*zeile,3);
for i=1:size(elements,1)
t(2*i-1,:)=[elements(i,1:2),elements(i,4)];
t(2*i,:)=[elements(i,2:4)];
end% t=[1 2 5
%    5 4 1
%    2 3 6
%    6 5 2
%    4 5 8
%    4 8 7
%    5 6 9
%    5 8 9]
trisurf(t,x,y,z)
xlabel('x'),ylabel('y'),zlabel('f(x,y)')
end% [x,y] = meshgrid(-1:1,-1:1);
% z = x.*x+y.*y;
% t = delaunay(x,y);
% subplot(1,2,1),trisurf(t,x,y,z)
% subplot(1,2,2),mesh(x,y,z)

MATLAB的一些簡單命令 1

關鍵字 if a y prctile input linspace disp 注釋 本部落格講解了一些matlab常用簡單命令的含義，這些命令簡單卻容易忘記，所以便整合了一下。至於較為複雜的命令，則將會單獨寫部落格記錄。為什麼 段設定選項沒有matlab，理學院學生emmm，1 否定用 表示，比如 ...

MATLAB讀取txt的一些方法簡單記錄

文件內容是數字和字元混合的時 1 比如文件testreaddata.txt 為 1 a w 2 c 3 7 b d 8 e 9 3 f s 6 r 5 m 3 n 6 data cell m,n 定義cell矩陣，儲存檔案內容 fid fopen d projects matlabproject d...

matlab的一些語法知識

matlab作為 工具，如今已經成為工科生必會的工具，其語法和c語言很相似，那麼就來複習一下吧。round 函式和fix 函式的作用都是取整，區別是，round是四捨五入，而fix是都捨掉了，即向0取整。還有floor 是不大於該數的最接近的整數，ceil 為不小於該數的最接近的整數。除是整個相除，...