Matlab小练

Matlab小练

题目1

已知蝴蝶曲线方程为：

请给出代码实现下图的绘制：

三个蝴蝶的中心为(1,1), (15,15), (30,30)，s=1

代码

clear all;
clc;
syms t f x y x0 y0 s
for i=1:3
t=0:0.1:40*pi;
s=1;
f=exp(cos(t))-2*cos(4*t)-(sin(t/12)).^(5);
if i==1
    x0=1;y0=1;
elseif i==2
    x0=15;y0=15;
else 
    x0=30;y0=30;
end
x=x0+s.*f.*sin(t);
y=y0+s.*f.*cos(t);
plot(x,y);
hold on
end

图像展示

题目2

使用for...end循环编写程序，求出0, 10, 20, ..., 200的正平方根，依次放入列向量X中。

代码

clear all;
clc;
X=[];
for i=0:10:200
    X=[X; i^(0.5)];
end
X

结果展示

题目3

已知方程组

判断方程组是否有唯一解，若有给出唯一解。要求：使用矩阵除法方法。

代码

clear all;
clc;

%系数矩阵
a=[2 1 -1; 1 -4 3; 5 2 -5];
b=[3; 2; 1];
%增广矩阵
c=[a b];

ra=rank(a);
rc=rank(c);
d=size(a);

if (ra==rc && ra==d(2))
    fprintf('方程有唯一解\n');
end;
x=inv(a)*b

结果展示

题目4

使用rand()函数生成4行4列的随机矩阵a。提取a中第1和3列元素构成a1，提取a1的第2和4行元素构成a2。最后计算a2的最小值。

代码

clear all;
clc;
a=rand(4)
a1=[a(:,1) a(:,3)]
a2=[a1(2,:);a1(4,:)]
min(min(a2))

结果展示