Matlab小练

最新推荐文章于 2023-08-19 21:43:54 发布

没没没没没超疯还是疯了

最新推荐文章于 2023-08-19 21:43:54 发布

阅读量349

点赞数

分类专栏： Matlab 文章标签： matlab 线性代数矩阵

本文链接：https://blog.youkuaiyun.com/m0_63033104/article/details/127600726

版权

Matlab 专栏收录该内容

2 篇文章

订阅专栏

博客围绕Matlab进行小练习，包含四个题目。一是实现蝴蝶曲线方程绘图；二是用循环求正平方根并放入列向量；三是用矩阵除法判断方程组是否有唯一解并求解；四是生成随机矩阵，提取元素构成新矩阵并计算最小值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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))

结果展示

在这里插入图片描述