1701007216
1701007217
运行程序输出图形如图5-53所示。
1701007218
1701007219
1701007220
1701007221
1701007222
图5-53 静态的五星红旗
1701007223
1701007224
该红旗是一面静态的图像，那么当五星红旗随风飘扬，红旗是一种波浪状的不断飘扬。采用MATLAB程序进行飞扬的五星红旗模拟，直观的一睹五星红旗飘扬的神态，具体程序如下：
1701007225
1701007226
clc,clear,close all %清屏和清除变量 warning off %消除警告 flag_size = [400 600]; %窗口大小 flag_texture = imread(‘FiveStarRedFlag.png’); %加载图像 flag_texture = double(flag_texture); %双精度 %%Constants cycle = 2; cycle2 = 3; frames = 10; %动画，10帧 wave_length = flag_size(1)./cycle; wave_length2 = flag_size(2)./cycle2; rot_angle = 15.*pi/180; rot_angle2 = 5.*pi/180; amp1 = 1; %幅度 amp2 = 0.5; %%Main Program flag_figure_handle = gcf; [flag_x flag_y] = meshgrid(1:flag_size(2), 1:flag_size(1)); %栅格化 for i = 1:frames cla; flag_z = amp1.*sin((cos(rot_angle).*flag_x - sin(rot_angle).*flag_ y).*2*pi./wave_length … - 2*pi*(i-1)/frames)+… amp2.*sin((cos(rot_angle2).*flag_x-sin(rot_ angle2).*flag_y).*2*pi./wave_length2 … - 2*pi*(i-1)/frames)+… sin((cosd(20).*flag_x-sind(20).*flag_y).*2*pi./600 - 2*pi*(i-1)/frames); end flag_surf_handle = warp(flag_x,flag_y,flag_z,flag_texture); view(2); %视图 set(flag_surf_handle, ‘FaceLighting’, ‘phong’,‘AmbientStrength’, 0.7); set(gca, ‘Visible’, ‘off’, ‘XLim’, [1 600], ‘YLim’, [1 400], ‘ZLim’, [-3 3]); camproj(‘perspective’) material dull
1701007227
1701007228
运行程序输出图形如图5-54所示。
1701007229
1701007230
1701007231
1701007232
1701007233
图5-54 飞扬的五星红旗
1701007234
1701007235
五星红旗随风飘扬，其波动情况和实际运动情况相符合，从程序角度来看实际，我们可以坐观风云，这就是程序下的几何图形之美，让我们开动思维，用我们的智慧筑起智能化设计的桥梁。
1701007236
1701007237
1701007238
1701007239
1701007240
我和数学有约：趣味数学及算法解析 [:1701004218]
1701007241
我和数学有约：趣味数学及算法解析 5.9 一颗红心永流传
1701007242
1701007243
一颗红心永流传，顾名思义，我肯定将描述一颗红心。红心就是我们眼中看到的像“心”一样的“类模型”，我们当中的很多人都喜欢它，一方面它精致，另一方面更是勿忘心安。一颗红心既然这么惹人喜欢，拥有一颗我们自己的红心，就显得那么理所当然了。
1701007244
1701007245
采用程序设计，我们可以设计一颗完美的红心，为我们所欣赏、所拥有，要想绘制三维的红心，首先来看看二维的心形图如何绘制。
1701007246
1701007247
二维红心的平面方程如下：
1701007248
1701007249
1701007250
1701007251
1701007252
采用程序进行数值解模拟如下：
1701007253
1701007254
figure(1) ezplot(‘(x^2-1+y^2)^3-x^2*y^3’) axis([-2.5 2.5 -2 2])
1701007255
1701007256
运行程序输出图形如图5-55所示。
1701007257
1701007258
1701007259
1701007260
1701007261
图5-55 二维心形图
1701007262
1701007263
二维心形图很逼近我们平时所看到的三维心形图的截面图，由二维图形，我们进行三维设计，就显得比较容易。具体的红心设计MATLAB程序如下：
1701007264
1701007265
clc,clear,close all %清屏和清除变量 warning off %消除警告 figure(‘color’,[1,1,1]) %设置背景为白色 x = linspace(-2,2,100); %均匀化划分点数 [X,Y,Z] = meshgrid(x,x,x); %平面栅格化 I1 = (X.^2+9/4*Y.^2+Z.^2-1).^3-X.^2.*Z.^3-9/80*Y.^2.*Z.^3; p = patch(isosurface(X,Y,Z,I1,0)); set(p, ‘FaceColor’, ‘red’, ‘EdgeColor’, ‘none’); view(3); %视图 axis equal ; %坐标轴属性，x和y等长 axis off; %不显示坐标轴 light(‘Posi’,[0 -2 3]); %在(0,-2,3)点处建立一个光源 lighting phong set(gca,‘nextplot’,‘replacechildren’); %增强显示颜色亮度 brighten(0.6) %亮度 campos([2,13,10]); %照相机位置 axis off vis3d; %隐藏坐标轴 view([9.93997681865585 29.6119315085364]); %视图