m文件
也就是matlab腳本
曹阔,很簡單,我們先來寫一個(gè):
進(jìn)入matlab,在command window
輸入edit
熄求,進(jìn)入editor:
寫一個(gè)m文件如下
function f = myfact(n)
f = prod(1:n);
end
點(diǎn)擊save
,命名為jiecheng.m
逗概,存一下弟晚。
這個(gè)程序是用來計(jì)算階乘
的
回到 command window,輸入jiecheng(5)
,你會(huì)看到返回了5!
這里,m文件中卿城,prod
函數(shù)用于求數(shù)組元素的乘積枚钓。意思是product.
排列數(shù),A(n,m)表示從n個(gè)東西內(nèi)選擇m個(gè)做排列瑟押,計(jì)算方法是n! ÷ (n-m)!
,我們寫成matlab函數(shù)如下:
function f = myfunc(n,m)
f = prod(1:n)/ prod(1:n-m);
end
存成pailie.m
搀捷,請(qǐng)自己測試以下這個(gè)函數(shù),pailie(5,1)
多望。
組合數(shù)的計(jì)算方式是C(n,m) = A(n,m) ÷ m!
請(qǐng)自己編寫函數(shù)測試一下嫩舟。
如果你看到了我寫的這篇文章,但是不打算去測試一下后面的兩個(gè)函數(shù)怀偷。
那么你的性格不適合學(xué)習(xí)計(jì)算機(jī)家厌,請(qǐng)轉(zhuǎn)系吧。去一個(gè)不需要學(xué)習(xí)計(jì)算機(jī)的專業(yè)椎工,對(duì)你自己會(huì)更好饭于。這是我誠懇的建議。
還有一個(gè)官方的視頻维蒙,你可以看看:
Writing a MATLAB Program (4 min, 57 sec)
下面有一些別的函數(shù)掰吕,你自己試試吧:
function out=c(a, b, theta)
out = sqrt(a.^2+b.^2-2*a.*b.*cos(theta));
%余弦定理計(jì)算第三邊
你計(jì)算:SideThree = c(2, 3, pi/6)
會(huì)得到:
SideThree =
1.6148
這里是計(jì)算最大公約數(shù)的m文件:
function y = gcd(a,b)
% GCD Greatest common divisor
% gcd(a,b) is the greatest common divisor
% of the integers a and b, not both zero.
a = round(abs(a)); b = round(abs(b));
if a == 0 & b == 0
error('The gcd is not defined when both numbers are zero')
else
while b ~= 0
r = rem(a,b);
a = b; b = r;
end
end
你要記茁:這里只是一個(gè)例子,MATLAB很強(qiáng)大八千,實(shí)際上:你直接gcd(a,b)
即可。
2015年5月4日17:03:32
@kaikai