九层妖塔 起于垒土


Matlab 李萨如图模拟_动态模拟李萨如图形形成过程保存为.GIF动图

  • 一、数学分析
  • 二、程序设计及使用说明
  • 三、动图
  • 四、代码

一、数学分析

二、程序设计及使用说明




三、动图


四、代码


%此文件中的代码为M文件编程
close all;
clc;
clf;
axis([-1.5,1.5,-1.5,1.5])
axis manual
grid on
grid minor
title('李萨如图模拟')
xlabel('X轴')
ylabel('Y轴')A=1;
B=1;
w=1;
Fx=0;                               %对定值提前赋值
t=linspace(0,4*pi,500);             %定义t为0到4π有500个元素的行矩阵
mkdir 'C:\Users\Lenovo\Desktop\'/ 李萨如图形  %在桌面的绝对路径下新建文件夹
for c=1:54                          %循环54次,生成54个图片close all;pause(3);switch c                        case 1q=1;p=1;Fy=0;                                                 %选择变量值mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率1比1         %新建文件夹   oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1'); %改变当前路径,保存旧路径case 2q=1;p=1;Fy=pi/4;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 3q=1;p=1;Fy=pi/2;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 4q=1;p=1;Fy=3*pi/4;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 5q=1;p=1;Fy=pi;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 6q=1;p=1;Fy=5*pi/4;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 7q=1;p=1;Fy=3*pi/2;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 8q=1;p=1;Fy=7*pi/4;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 9q=1;p=1;Fy=2*pi;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率1比1');case 10q=2;p=1;Fy=0; mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率2比1            oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 11q=2;p=1;Fy=pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 12q=2;p=1;Fy=pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 13q=2;p=1;Fy=3*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 14q=2;p=1;Fy=pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 15q=2;p=1;Fy=5*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 16q=2;p=1;Fy=3*pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 17q=2;p=1;Fy=7*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 18q=2;p=1;Fy=2*pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率2比1');case 19q=3;p=1;Fy=0; mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率3比1            oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 20q=3;p=1;Fy=pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 21q=3;p=1;Fy=pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 22q=3;p=1;Fy=3*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 23q=3;p=1;Fy=pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 24q=3;p=1;Fy=5*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 25q=3;p=1;Fy=3*pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 26q=3;p=1;Fy=7*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');case 27q=3;p=1;Fy=2*pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比1');  case 28q=3;p=2;Fy=0/q; mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率3比2            oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 29q=3;p=2;Fy=pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 30q=3;p=2;Fy=pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 31q=3;p=2;Fy=3*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 32q=3;p=2;Fy=pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 33q=3;p=2;Fy=5*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 34q=3;p=2;Fy=3*pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 35q=3;p=2;Fy=7*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');case 36q=3;p=2;Fy=2*pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率3比2');  case 37q=4;p=3;Fy=0/q; mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率4比3            oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 38q=4;p=3;Fy=pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 39q=4;p=3;Fy=pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 40q=4;p=3;Fy=3*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 41q=4;p=3;Fy=pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 42q=4;p=3;Fy=5*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 43q=4;p=3;Fy=3*pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 44q=4;p=3;Fy=7*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');case 45q=4;p=3;Fy=2*pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率4比3');         case 46q=5;p=3;Fy=0/q; mkdir 'C:\Users\Lenovo\Desktop\李萨如图形\'/  频率5比3            oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 47q=5;p=3;Fy=pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 48q=5;p=3;Fy=pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 49q=5;p=3;Fy=3*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 50q=5;p=3;Fy=pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 51q=5;p=3;Fy=5*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 52q=5;p=3;Fy=3*pi/2/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 53q=5;p=3;Fy=7*pi/4/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');case 54q=5;p=3;Fy=2*pi/q;oldFolder = cd('C:\Users\Lenovo\Desktop\李萨如图形\频率5比3');  endfor n=1:length(t)plot(A*cos(Fx),B*cos(Fy),...                    %绘制起始点'Color','green',...   'Marker','.',...'linesty','-',...'LineWidth',25,...'MarkerFaceColor','green',...'MarkerEdgeColor','green',...'MarkerSize',25);hold onplot(A*cos(q*w*t(n)+Fx),B*cos(p*w*t(n)+Fy),...   %绘制当前点'Color','black',...   'Marker','^',...'linesty','-',...'LineWidth',1,...'MarkerFaceColor','black',...'MarkerEdgeColor','black',...'MarkerSize',10);axis([-1.5,1.5,-1.5,1.5])     axis manualhold ongrid ongrid minortitle('李萨如图模拟')xlabel('X轴')ylabel('Y轴')plot(A*cos(q*w*t(1:n)+Fx),B*cos(p*w*t(1:n)+Fy),... %绘制行进的曲线'Color','red',...   'Marker','.',...'linesty','-',...'LineWidth',3.0,...'MarkerFaceColor','red',...'MarkerEdgeColor','red',...'MarkerSize',1);axis([-1.5,1.5,-1.5,1.5])     axis manualgrid ongrid minortitle('李萨如图模拟')xlabel('X轴')ylabel('Y轴')hold offaxis manualgrid ongrid minorframe=getframe(gcf);                         %捕获坐标区或图窗作为影片帧imind=frame2im(frame);[imind,cm] = rgb2ind(imind,256);             %将 RGB 图像转换为索引图像if n==1                                          %将图像写入图形文件imwrite(imind,cm,'name.gif','gif', 'Loopcount',inf,'DelayTime',1e-4);elseimwrite(imind,cm,'name.gif','gif','WriteMode','append','DelayTime',1e-4);endendswitch c                                      %依次选择,将图片改名case 1  movefile ('name.gif', '频率1比1 相位差0.gif') ;case 2movefile ('name.gif', '频率1比1 相位差45.gif') ;case 3  movefile ('name.gif', '频率1比1 相位差90.gif') ;case 4  movefile ('name.gif', '频率1比1 相位差135.gif') ;case 5movefile ('name.gif', '频率1比1 相位差180.gif') ;case 6 movefile ('name.gif', '频率1比1 相位差225.gif') ;case 7  movefile ('name.gif', '频率1比1 相位差270.gif') ;case 8movefile ('name.gif', '频率1比1 相位差315.gif') ;case 9 movefile ('name.gif', '频率1比1 相位差360.gif') ;case 10  movefile ('name.gif', '频率2比1 相位差0.gif') ;case 11movefile ('name.gif', '频率2比1 相位差45.gif') ;case 12  movefile ('name.gif', '频率2比1 相位差90.gif') ;case 13movefile ('name.gif', '频率2比1 相位差135.gif') ;case 14movefile ('name.gif', '频率2比1 相位差180.gif') ;case 15 movefile ('name.gif', '频率2比1 相位差225.gif') ;case 16  movefile ('name.gif', '频率2比1 相位差270.gif') ;case 17movefile ('name.gif', '频率2比1 相位差315.gif') ;case 18 movefile ('name.gif', '频率2比1 相位差360.gif') ;    case 19  movefile ('name.gif', '频率3比1 相位差0.gif') ;case 20movefile ('name.gif', '频率3比1 相位差45.gif') ;case 21  movefile ('name.gif', '频率3比1 相位差90.gif') ;case 22  movefile ('name.gif', '频率3比1 相位差135.gif') ;case 23movefile ('name.gif', '频率3比1 相位差180.gif') ;case 24 movefile ('name.gif', '频率3比1 相位差225.gif') ;case 25 movefile ('name.gif', '频率3比1 相位差270.gif') ;case 26movefile ('name.gif', '频率3比1 相位差315.gif') ;case 27 movefile ('name.gif', '频率3比1 相位差360.gif') ;    case 28 movefile ('name.gif', '频率3比2 相位差0.gif') ;case 29movefile ('name.gif', '频率3比2 相位差45.gif') ;case 30  movefile ('name.gif', '频率3比2 相位差90.gif') ;case 31  movefile ('name.gif', '频率3比2 相位差135.gif') ;case 32movefile ('name.gif', '频率3比2 相位差180.gif') ;case 33 movefile ('name.gif', '频率3比2 相位差225.gif') ;case 34 movefile ('name.gif', '频率3比2 相位差270.gif') ;case 35movefile ('name.gif', '频率3比2 相位差315.gif') ;case 36movefile ('name.gif', '频率3比2 相位差360.gif') ;        case 37movefile ('name.gif', '频率4比3 相位差0.gif') ;case 38movefile ('name.gif', '频率4比3 相位差45.gif') ;case 39  movefile ('name.gif', '频率4比3 相位差90.gif') ;case 40  movefile ('name.gif', '频率4比3 相位差135.gif') ;case 41movefile ('name.gif', '频率4比3 相位差180.gif') ;case 42 movefile ('name.gif', '频率4比3 相位差225.gif') ;case 43movefile ('name.gif', '频率4比3 相位差270.gif') ;case 44movefile ('name.gif', '频率4比3 相位差315.gif') ;case 45movefile ('name.gif', '频率4比3 相位差360.gif') ;            case 46movefile ('name.gif', '频率5比3 相位差0.gif') ;case 47movefile ('name.gif', '频率5比3 相位差45.gif') ;case 48  movefile ('name.gif', '频率5比3 相位差90.gif') ;case 49  movefile ('name.gif', '频率5比3 相位差135.gif') ;case 50movefile ('name.gif', '频率5比3 相位差180.gif') ;case 51movefile ('name.gif', '频率5比3 相位差225.gif') ;case 52movefile ('name.gif', '频率5比3 相位差270.gif') ;case 53movefile ('name.gif', '频率5比3 相位差315.gif') ;case 54movefile ('name.gif', '频率5比3 相位差360.gif') ;          endend

Matlab 李萨如图模拟_动态模拟李萨如图形形成过程保存为.GIF动图相关推荐

  1. 【1.6万字】连续抓屏保存为Gif动图 【keyboard库、PIL库、imageio库和pygifsicle库 探索】

    目录 一.抓屏保存为Gif 二.keyboard库 三.PIL 库 3.0 安装Pillow 3.1 打开本地图片 3.2 创建一张新图片 3.3 Image模块的常用属性 3.4 图片的模式和模式转 ...

  2. R语言gganimate包创建可视化gif动图、并使用anim_save函数保存可视化gif动图(gganimate export to gif)

    R语言gganimate包创建可视化gif动图.并使用anim_save函数保存可视化gif动图(gganimate export to gif) 目录

  3. 随机模拟_随机模拟可帮助您掌握统计概念

    随机模拟 模拟有助于提炼概念 (Simulation helps distilling concepts) 掌握与统计相关的概念可能很困难 (Grasping statistics-related c ...

  4. python有趣的简单代码_简单代码一学就会,Python生成GIF动图

    文/IT可达鸭 图/IT可达鸭.网络 前言 最近在写技术文档的时候,发现一个问题.对于每个技术步骤,都需要一个截图,这样下来整篇文档都是截图,显得不是特别的专业. 为了解决这个问题,我想到一个方法,就 ...

  5. python的gif为正方形_简单代码一学就会,Python生成GIF动图

    文/IT可达鸭 图/IT可达鸭.网络 前言 最近在写技术文档的时候,发现一个问题.对于每个技术步骤,都需要一个截图,这样下来整篇文档都是截图,显得不是特别的专业. 为了解决这个问题,我想到一个方法,就 ...

  6. python制作gif动画_实用的Python(2)利用Python制作gif动图

    一.简介 moviepy是一个专门用于视频剪辑制作的模块,可以自动化完成很多繁琐的视频剪辑处理工作,除了处理视频数据之外,moviepy中还内置了可以制作gif动图的功能,通过使用moviepy.ed ...

  7. 动态图片怎么做?教你一键合成gif动图

    在浏览网页的时候看到一些搞笑有趣的图片,想要制作成gif动图的效果,但是不会使用PS的时候,应该如何制作呢?下面,小编为大家推荐一款简单实用的在线图片制作(https://www.yasuotu.co ...

  8. lisp调用布置天正的图块_操作太慢?那是你没有正确使用AutoCAD图块!

    不容置疑,会用CAD的朋友,都会用CAD中的图块功能,但是,对CAD图块的认识深度可能就不尽相同.熟练掌握图块特性和使用图块绘图,是每一个渴望成为AutoCAD高手必备的利器.虽然组成图块的各对象都有 ...

  9. ImageMagick将多张图片拼接成一张图片_如何快速将多张图片或视频制作成GIF动图表情...

    相信大家看到朋友手上很多搞笑和经典的动图表情都会收藏来自己使用,今天软件砖家就给大家介绍一下如何把多张图片或视频制作成GIF动图表情! 开启万兴优转,然后切换至工具箱的标签页然后点击GIF动图制作. ...

最新文章

  1. 基于Spring Boot+Cloud构建微云架构
  2. Android-环境问题
  3. MIPS指令以及数据通路
  4. 光端机和光纤交换机的区别?
  5. 八、计数排序及其应用分析
  6. [django]django模型中auto_now和auto_now_add
  7. 使用numpy实现神经网络模块
  8. 【java基础知识】Spring Boot启动报错com.sun.jersey.api.client.ClientHandlerException: java.net.ConnectExcepti
  9. flutter英语怎么说_美国人天天说的英语:“你搞反了”英语怎么说?
  10. 据说集齐这些特点的智能硬件百分百“作死”
  11. 抖音自动清关注服务器,各位有保留抖音自动关注的代码吗
  12. 3S基础知识:用MapX快速开发
  13. IM 即时通讯开发如何设计图片文件的服务端存储架构
  14. 【板栗糖GIS】DOS—如何在当前文件夹内部批量建子文件夹
  15. 微信小程序开发工具第一次使用,网络连接不上怎么办
  16. VC2010 无法启动程序 系统找不到指定文件
  17. Pushed master to new branch origin/master
  18. 【2022年的STM32】 04-GPIO特性、使用及与NXP GPIO比较
  19. java 素数求和,100-200之间所有素数求和程序代码(二个版本)
  20. Unity质量设置 Quality Settings

热门文章

  1. SVS打造智能指挥中心成功案例
  2. 什么是剩余电流动作继电器?
  3. js敏感信息泄露检测工具
  4. Latex中文配置中修复Package fontspec Error
  5. 阿里云部署redis
  6. 数字图像处理基础知识
  7. 南京理工大学计算机科学与工程学院复试名单,南京理工大学计算机科学与工程学院2015考研复试及录取办法...
  8. python画八角星_Python画八角星的程序
  9. 河南发文明确!这类企业数字化转型有了“时间表”
  10. php里drear_精通Dreamweaver CS3网站建设PHP篇