收藏几个经典的函数,用来当作数学优化方法求解的例子,再收藏几个有趣的函数图像。

一、平面心形图像

1. 笛卡尔心形函数

第一个当然是大名鼎鼎的笛卡尔心形函数,它的平面直角坐标系方程为:

\begin{equation}(x^2+y^2)^2+4ay(x^2+y^2)-4a^2x^2=0\end{equation}

其中 a 是一个可调参数,下面的图为 a=1 时的图像。

matlab 代码:

ezplot('(x^2+y^2)^2+4*2*x*(x^2+y^2)-4*2^2*y^2=0')

笛卡尔心形图比较像一个心脏的形状。

2 函数1

经过网上搜索,平面坐标系中,下面的函数最像普通的心形图像,而且函数形式简单,不需要分段。

\begin{equation}x^2+(y-(x^2)^{1/3})^2=9 \end{equation}

改变右面的数值可以变化心形图像的大小,MATLAB 代码:

> ezplot('x^2+(y-(x^2)^(1/3))^2=9')

3. 函数2

还有一个函数:

\begin{equation}-x^2y^3+(x^2+y^2-1)^3=0\end{equation}

代码:

ezplot('-x^2*y^3+(x^2+y^2-1)^3=0',[-1.5,1.5])

4. 函数3

\begin{equation}17x^2-16|x|y+17y^2=200\end{equation}

代码:

ezplot('17*x.^2-16*abs(x).*y+17*y.^2=200'

5. 函数4

\begin{align}f(x)&=\sqrt{2\sqrt{x^2}-x^2}\nonumber\\g(x)&=-2.14\sqrt{\sqrt{2}-\sqrt{|x|}}\nonumber\end{align}

代码:

x=linspace(-2,2,1000);
y1=sqrt(2*sqrt(x.^2)-x.^2);
y2=-2.14*sqrt(sqrt(2)-sqrt(abs(x)));
plot(x,y1,'b',x,y2,'b');
axis([-2.5,2.5,-3,1.5]);

6. 函数5

\begin{align}x&=16(sint)^3\nonumber\\y&=13cost-5cos(2t)-2cos(3t)-cos(4t)\nonumber\end{align}

t=linspace(-6,6,1000);
x=16*(sin(t)).^3;
y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t);
plot(x,y);

二 . 立体心形图像

立体心形貌似只能由笛卡尔方程得到:

\begin{equation}(x^2+ 9y^2/4 + z^2- 1)^3 - x^2z^3 - 9y^2z^3/80=0\end{equation}

画图时用 isosurface 和 patch  函数:

1. 图形1

f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;
[x,y,z]=meshgrid(linspace(-3,3));
val=f(x,y,z);
[p,v]=isosurface(x,y,z,val,0);
patch('faces',p,'vertices',v,'facevertexcdata',jet(size(v,1)),'facecolor','w','edgecolor','flat');
view(3);
grid on;
axis equal;

2. 图形2

f=@(x,y,z)(x.^2+ (9./4).*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9./80).*y.^2.*z.^3;
[x,y,z]=meshgrid(linspace(-1.5,1.5));
val=f(x,y,z);
isosurface(x,y,z,val,0);
axis equal;
view(3);
colormap([1 0.2 0.2])

3. 轴对称振荡器函数

\begin{equation}f(x)=xsin(x)\end{equation}

ezplot('x*sin(x)')

4.  刚好有一个极大点,一个极小点的二元函数

\begin{equation}f(x,y)=xe^{-x^{2}-y^{2}}\end{equation}

ezmesh(@(x,y) x.*exp(-x.^2-y.^2))

5. peaks 函数(二元高斯分布的概率密度函数)(三个极大点,三个极小点)

\begin{equation}f(x,y)=3(1-x)^{2}e^{-x^2-(y+1)^2}-10(\frac{1}{5}x-x^3-y^5)e^{-x^2-y^2}-\frac{1}{3}e^{-(x+1)^2-y^2}\end{equation}

 f=@(x,y)3*(1-x).^2.*exp(-(x.^2) - (y+1).^2)- 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2)- 1/3*exp(-(x+1).^2 - y.^2);ezmesh(f);

转载于个人公众号:Python 统计分析与数据科学

matlab心形图大全,几个经典的函数图像,有趣的函数图像,matlab画图相关推荐

  1. matlab心形图c语言,MATLAB心形图.doc

    MATLAB心形图 三维图 表达你对他的爱 [X,Y,Z] = meshgrid(linspace(-3,3,101)); F = -X.^2.*Z.^3-(9/80).*Y.^2.*Z.^3+(X. ...

  2. [MATLAB] 心形图

    效果图: 代码: f=@(x,y,z)(x.^2+(9.4).*y.^2+z.^2-1).^3-x.^2.*z.^3-(9./80).*y.^2.*z.^3; [X,Y,Z]=meshgrid(lin ...

  3. matlab怎么画一箭穿心,MATLAB学习与使用:如何绘制三维心形图 经验告诉你该这样...

    MATLAB是MATrix & LABoratory(矩阵实验室)的缩写,是一款强大的科学软件,具有编程.绘图.仿真等功能.利用MATLAB绘制一款三维的心形图,然后送给心爱的姑娘,也是理工男 ...

  4. 【matlab】 动态绘制三维心形图

    [转载]matlab动态绘制三维心形图 matlab动态绘制三维心形图程序如下,核心是isosurface函数,verts存储了心形曲面各点的坐标,而faces存储了各点的连接顺序. clear fi ...

  5. matlab绘图——根据方程组绘图和心形图绘制

    目录 1.如何绘图 1.1标签设置 1.2图像颜色参数 1.3图像线形参数 : 2.具体例子 2.1.根据方程组绘制图像 2.2.绘制心形图 3.mesh函数 3.1 meshgrid() 3.2 m ...

  6. 三维旋转心形图(matlab 含源码及注释)

    一.二维心形线 (1)笛卡尔心形线,分为开口向左.右.上.下的四种.可用极坐标方程画图,在此不做演示.读者感兴趣或需要可自行百度查找. (2)第二种心形线 方程: 实现代码: ezplot('(x^2 ...

  7. 【Python】函数图像绘制:二维图像、三维图像、散点图、心形图

    [Python]函数图像绘制:二维图像.三维图像.散点图.心形图 所有需要用的包 二维图像 三维图像 散点图绘制 心形图绘制 所有需要用的包 from mpl_toolkits.mplot3d imp ...

  8. MATLAB之心形图绘制

    一.静态心形图绘制 (1)效果展示 (2)静态心形原始代码 1 clc; 2 clear all; 3 const=0; 4 % 均布三位坐标 5 x=-5:0.05:5; 6 y=-5:0.05:5 ...

  9. 使用R语言绘制心形图

    今天七夕,正好看到高等数学的心形线,想到心形线的函数应该可以用R语言来绘制,就尝试了一下. 心形线的参数方程为: 其中r是半径,t为弧度. 有了参数方程之后,我们的作图就变得简单了,其基本思路是:首先 ...

最新文章

  1. 数据库死锁查询及处理
  2. abap--关于异常的处理
  3. 工具-eclipse-hibernate关于快速生成*.hbm.xml文件,与bean文件的操作,不使用hibernate tool解决办法之一,主要为了快速开发,写的小工具
  4. 通过福禄克LinkIQ网络电缆测试仪排除 VLAN 故障
  5. cvtcolor python opencv_二值分析 | OpenCV + skimage如何提取中心线
  6. Leetcode--287. 寻找重复数(Java)
  7. python接口自动化(三十二)--Python发送邮件(常见四种邮件内容)番外篇——上(详解)...
  8. 领扣(LeetCode)七进制数 个人题解
  9. 《软件工程导论第6版》--张海藩 牟永敏 课后答案及其详解 第4章 形式化说明技术
  10. AI笔记: 计算机视觉之图像边缘检测: Robert算子, Sobel算子, Laplace算子, LoG算子, Canny算子
  11. 园区网络IP地址规划
  12. html半透明遮罩,div半透明遮罩效果
  13. python项目:基于OpenCV的学生网课睡意检测系统
  14. MyEclipse2017CI安装包及图解破解步骤
  15. 安装谷歌插件,程序包无效,解决方法
  16. 微软dpm服务器,使用 DPM 的 ReFS 卷在服务器中Windows Server 2016
  17. Spring Mvc时间格式参数,前后时间传参,单个字段的时间格式化问题
  18. Python爬虫实战:手机APP数据抓取分析!谁说不能爬取app数据的?
  19. Oracle 查看和修改数据库时区
  20. USB PD v1.0快速充电通信原理

热门文章

  1. 海尔银悦扫地机器人怎样规划_海尔银悦智能扫地机器人TAB-T550WSC 陀螺仪导航定位规划家用全自动洗擦拖地机吸尘器一体机...
  2. python身份证号判断性别。1.判断身份证长度是否正确 2.判断性别
  3. R语言pheatmap包热图legend位置调整
  4. Word表格占满整页,怎样删除后面的回车和空白页?
  5. 物联网毕业设计 - 单片机宠物远程智能管理系统
  6. 黑龙江科技大学计算机推免,黑龙江科技大学计算机应用技术怎么样
  7. 计算机b类b级2020,2020年A级、B级驾驶证降级新规定
  8. Pandas基本操作教程
  9. 专业心理测试软件吗,您真的了解心理测评软件吗?
  10. Android Studio实现飞机大战