利用Matlab绘制各类特殊图形的实例代码
作者:MatlabFans_Mfun 发布时间:2021-01-05 21:10:31
Matlab绘图介绍
强大的绘图功能是Matlab的特点之一,Matlab提供了一系列的绘图函数,用户不需要过多的考虑绘图的细节,只需要给出一些基本参数就能得到所需图形,这类函数称为高层绘图函数。此外,Matlab还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素(如坐标轴、曲线、文字等)看做一个独立的对象,系统给每个对象分配一个句柄,可以通过句柄对该图形元素进行操作,而不影响其他部分。
Matlab绘制特殊图形
1. 绘制极坐标图
说明:使用polarplot函数绘制极坐标图,每组数据表示一条闭合曲线,共有20条曲线构成20条封闭同心曲线。
t = linspace(0,2*pi,500);
y = 1+0.3*sin(20*t)+0.1*sin(30*t)+0.1*sin(40*t);
polarplot(t,y,t,1.1*y,t,1.2*y,t,1.3*y,t,1.4*y,t,1.5*y,t,1.6*y,t,1.7*y,t,1.8*y,t,1.9*y,...
t,2.0*y,t,2.1*y,t,2.2*y,t,2.3*y,t,2.4*y,t,2.5*y,t,2.6*y,t,2.7*y,t,2.8*y,t,2.9*y,'linewidth',1.5);
2. 单条曲线绘制分段函数(反比例函数y=1/x)
说明:反比例函数在x接近于0时,趋近于奇异(1/0 趋近于无穷大),使用nan非数对图形进行镂空,可以实现一条曲线绘制y=1/x的整个定义域。
x = linspace(-3,3,500);
f = @(x) 1./x;
y = f(x);
for ii = 1:length(x)
if abs(x(ii))<0.03
x(ii) = nan;
y(ii) = nan;
end
end
plot(x,y)
legend('y=1/x')
3. 正方体内绘制随机分布的颜色片图
说明: 使用plot3绘制正方体12条边,使用fill3函数绘制颜色片,颜色片位置和填充颜色随机生成。
a = 20;
b = 20;
c = 20;
A1 = [0 0 0
a 0 0
a b 0
0 b 0
0 0 0];
A2 = [0 0 c
a 0 c
a b c
0 b c
0 0 c];
A3 = [0 0 0
0 0 c
0 b c
0 b 0
0 0 0];
A4 = [a 0 0
a 0 c
a b c
a b 0
a 0 0];
% 绘图
figure
hold on
plot3(A1(:,1),A1(:,2),A1(:,3),'k');
plot3(A2(:,1),A2(:,2),A2(:,3),'k');
plot3(A3(:,1),A3(:,2),A3(:,3),'k');
plot3(A4(:,1),A4(:,2),A4(:,3),'k');
view(3)
XYZ0 = [0 0 0
1 0 0
1 0 1
0 0 1];
for ii = 1:100
p = [rand*(a-2) rand*(b-2) rand*(c-2)]+1; % 球心
X2= p(1) + XYZ0(:,1);
Y2 = p(2) + XYZ0(:,2);
Z2 = p(3) + XYZ0(:,3);
ColorSpec = rand(1,3);
fill3( X2,Y2,Z2,ColorSpec )
end
4. 在大圆内随机生成若干互不相交的小圆
说明:先绘制大圆,在大圆内随机生成小圆圆心,判断小圆圆心和其它小圆的距离,如果小圆圆心距离小于直径,说明小圆相交,舍去该小圆。
tic
X = [];
Y = [];
n = 0;
while n < 1000
r = rand * (750-5);
theta = rand * 2*pi;
x0 = r*cos(theta);
y0 = r*sin(theta);
s = min( (x0-X).^2 + (y0-Y).^2 ); % 最小圆心距离
if s < 10^2
continue;
else
X = [X;x0]; % 圆心坐标集合
Y = [Y;y0];
n = n +1 % 点数
end
end
toc
alpha = linspace(0,2*pi,100);
x = 5*cos(alpha);
y = 5*sin(alpha);
figure
plot( 150*x,150*y,'r' )
hold on
for ii = 1:n
x0 = X(ii);
y0 = Y(ii);
plot(x0+x,y0+y,'k')
end
5. 在长方体内随机生成若干球体
说明:patch函数绘制长方体和球体表面,并可填充颜色。
x = 100;
y = 80;
z = 50;
theta = linspace(0,2*pi,50);
phi = linspace(0,2*pi,50);
[theta,phi] = meshgrid(theta,phi);
r = 2;
X0 = r*cos(phi).*cos(theta);
Y0 = r*cos(phi).*sin(theta);
Z0 = r*sin(phi);
% 绘图
X = [0 x x 0
0 x x 0
0 0 0 0
x x x x
0 x x 0]';
Y = [0 0 y y
0 0 y y
0 y y 0
0 y y 0
y y y y]';
Z = [0 0 0 0
z z z z
0 0 z z
0 0 z z
0 0 z z ]';
figure
patch(X,Y,Z,'r');
view(3)
hold on
for ii = 1:50
p = [rand*(100-4*r) rand*(80-4*r) rand*(50-4*r)]+2*r; % 球心
X2= p(1) + X0;
Y2 = p(2) + Y0;
Z2 = p(3) + Z0;
patch( X2,Y2,Z2,'y' )
end
6. 绘制圆柱体与球体曲面相交,并绘制相交曲线
说明: mesh函数绘制曲面图,求解方程得到交线参数方程,plot3函数绘制三维交线。
%% 曲面1
t = linspace(0,2*pi,200);
s = linspace(0,2*pi,200);
[t,s] = meshgrid(t,s);
x = 2*cos(t);
y = 2*sin(t).*cos(s);
z = 2*sin(t).*sin(s);
figure
mesh(x,y,z)
%% 曲面2
t2 = linspace(0,2*pi,200);
z2 = linspace(-3,3,200);
[t2,z2] = meshgrid(t2,z2);
x2 = 1 + cos(t2);
y2 = sin(t2);
hold on
mesh(x2,y2,z2)
%% 交线
t3 = linspace(0,2*pi,200);
y3 = sin(t3);
x3 = 1 + cos(t3);
z3 = sqrt(4-2*x3);
plot3(x3,y3,z3,'r','linewidth',5)
hold on
plot3(x3,y3,-z3,'r','linewidth',5)
8. 绘制三维抛物曲面
说明:使用nan非数对图形进行镂空
x = -20:0.1:20;
y = -20:0.1:20;
[X,Y] = meshgrid(x,y);
p = 0.2;
q = 0.1;
Z = X.^2/(2*p) + Y.^2/(2*q);
Z = (Z<=500) .* Z + ((Z>500)-1) ./ ((Z>500)-1); % 图形镂空
mesh(X,Y,Z)
9. 抛物曲面随参数变化形成动画
说明:每个步长内动态更新绘制三维曲面,形成动画效果。
x = linspace(-1,1,20);
y = linspace(-1,1,20);
[X,Y] = meshgrid(x,y);
figure
a = 1;
Z = a.*X.^2 + Y.^2;
h = surf(X,Y,Z);
zlim([0,15])
for a = 1:0.1:10
Z = a.*X.^2 + Y.^2;
set(h,'zdata',Z);
drawnow
pause(0.1)
end
10. 使用不同频率的正弦波合成方波
说明:傅里叶级数,利用不同频率的正弦波合成方波,三角函数项数越多,合成方波越精确。
t=0:0.000001:1;
f1=6*sin(10*pi*t)/pi;
f2=6*sin(10*pi*t)/pi+2*sin(30*pi*t)/pi;
f3=6*sin(10*pi*t)/pi+2*sin(30*pi*t)/pi+6*sin(50*pi*t)/(5*pi);
%% 循环段
N = 10; % 点数
f4 = 0; % 初始值
for ii = 1:N
f4 = f4 + 3*2*sin((2*ii-1)*10*pi*t)/pi/(2*ii-1);
end
figure
subplot(2,2,1),plot(t,f1)
subplot(2,2,2),plot(t,f2)
subplot(2,2,3),plot(t,f3)
subplot(2,2,4),plot(t,f4)
总结
来源:https://blog.csdn.net/MatlabFans_Mfun/article/details/118738103
猜你喜欢
- 前言又见面了,小伙伴儿们,发现最近大家喜欢看一些简单的小案例?!咳咳,下面进入正题。每个人的电脑里都会有不想让别人知道的隐私,或者是上班时间
- 由于笔者最近在做一个跨数据库操作的测试,开始做IBatisNet (IBatis.DataMapper.1.6.2/IBatis.DataA
- 1. 直接打印错误try:# your codeexcept KeyboardInterrupt:print("quit"
- 废话不多说,直接上代码:# coding:utf-8from urllib import requestfrom urllib import
- 好久没有更新博客了,今天看到论坛上有位朋友问起全屏布局,有点像vc的界面。来了兴趣,就写了一个。运用IE6的怪异模式,通过绝对定位来实现的。
- Python装饰器(decorator)在实现的时候,被装饰后的函数其实已经是另外一个函数了(函数名等函数属性会发生改变),为了不影响,Py
- 一 前言温习python 多进程语法的时候,对 join的理解不是很透彻,本文通过代码实践来加深对 join()的认识。multiproce
- Internet Explorer 8 Beta 测试了一年多之后,今天,IE8 终于发布了。它绝对好过 IE7,还有一些不错的新功能,如
- pygame.mixer是一个用来处理声音的模块,其含义为“混音器”。游戏中对声音的处理一般包括制造声音和播放声音两部分,这里仅学习了播放声
- 在大多数语音识别任务中,我们都缺少文本和音频特征的alignment,Connectionist Temporal Classificati
- server:#coding=utf-8from BaseHTTPServer import BaseHTTPRequestHandleri
- 实际上,在web开发中,cookie仅仅是一个文本文件,当用户访问站点时,它就被存储在用户使用的计算机上,其中,保存了一些信息,当用户日后再
- 实现过程最近在尝试openai的相关函数,所以今天就来盘点盘点下 createChatCompletion 函数,听说是专门为了聊天创立的a
- 四年前写的一个内容管理系统,应用在公司内部网上,昨天DBA说其中的SQL语句未使用参数化的调用,导致服务器负担加重,资源占用大。并列出了几个
- --利用T-SQL语句,实现数据库的备份与还原的功能 ----体现了SQL Server中的四个知识点: ----1. 获取SQL Serv
- 1,file_put_contents()函数 2,使用PHP内置缓存机制实现页面静态化:output_bufferingphp中outpu
- 各大云计算提供商(亚马逊、谷歌和微软)目前都使用了键/值存储方式。然而,在San Francisco召开的MSDN开发者大会上,微软宣布他们
- 1.TCP是一种面向连接的可靠地协议,在一方发送数据之前,必须在双方之间建立一个连接,建立的过程需要经过三次握手,通信完成后要拆除连接,需要
- 代码如下所示:scole = input("input your scole:")if scole>90: &nb
- 前言相信对于python开发人员来说,机器上有不同的python版本是很正常的,因为开发的项目有的用2.6或2.7,有的就要用3.0+版本,