【matlab学习】matlab自定义函数的编写

判断一个数是否为素数

function p=myprime(a)
for i=2:a-1
    if 0==rem(a,i)
        p=0;
        break;
    elseif i==a-1
        p=1;
    end
end

(测试代码)

clc;
n = input('请输入一个数：');
res = myprime(n);
if res == 1
    disp('是素数')
else
    disp('不是素数')
end  

求两个数的最大公约数和最小公倍数

function z=mycalulate(a,b)
c = 1;
while c~= 0
    c = mod(a,b);
    a = b;
    b = c;
end
z = a;
end

(测试代码)

x = input('输入第一个数:');
y = input('输入第二个数:');
m1 = mycalulate(x,y);
m2 = x*y/m1;
disp(['最大公约数:' num2str(m1)]);
disp(['最小公倍数:' num2str(m2)]);

求两个矩阵的乘积

function C=work(A,B)
[m1,n1]=size(A);
[m2,n2]=size(B);
if(n1~=m2)
    C=[];
    C;
end
C=zeros(m1,n2);
for i=1:m1
    for j=1:n2
        for k=1:n1
            C(i,j)=C(i,j)+A(i,k)*B(k,j);
        end
    end
end
C;
end

(测试代码)

clc;
clear;
a=[1 3 5;2 4 6]
b=[1 2;4 5;6 7]
c=work(a,b)