行业知识
MATLAB | 儿童节一起来玩转MATLAB叭
2026-10-27 13:18:43
来源：互联网转载
hey, 儿童节要到啦，提前祝各位大朋友，小朋友儿童节快乐，本篇将会在前半部分给出一些有趣且简短的代码，在后半部分展示一下我准备的游戏大礼包里面的游戏效果，完整代码获取方式请见文末，让我们一起玩转MATLAB叭(请尽量使用R2018b及之后版本)：
有趣代码
该部分有一部分是我以前的代码进行的重写，有一部分是粉丝问我要的代码，还有一部分是对大佬代码的改写：
平面大钻石
% author : slandarerfigure('render','painters')hold on; axis equal; axis([-30,30,-30,30])CList = [21,97,169; 68,107,178; 169,217,246; 76,177,227; 92,134,195; 55,124,199]./255;set(gca, 'XColor','none', 'YColor','none', 'Color',[0,64,115]./510)t = 0:20:360-1;T = [cosd(t).*7; sind(t).*7; t];plot([1; 0]*T(1,:), [1; 0]*T(2,:), 'LineWidth',2)for i = 1:5 t = T(3,:); L = T + [cosd(t - 20)*(7 - i); sind(t - 20)*(7 - i); t*0 - 20]; R = T + [cosd(t + 20)*(7 - i); sind(t + 20)*(7 - i); t*0 + 20]; plot([T(1,:); L(1,:)], [T(2,:); L(2,:)], 'LineWidth',2) plot([T(1,:); R(1,:)], [T(2,:); R(2,:)], 'LineWidth',2) T = [L, R];endLHdl = findobj(gca, 'Type','line');for i = 1:length(LHdl) LHdl(i).Color = CList(randi([1,6],[1,1]),:);endP = T(1:2,:).';[k,~] = convhull(P);plot(P(k,1), P(k,2), 'LineWidth',2, 'Color',CList(1,:)) 羽毛
% author : slandarer% remixed from Juan Villacrés / Flower R-Dfigure('Color','k', 'Renderer','painters', 'InvertHardCopy', 'off')hold ont = linspace(-3, 3, 1000);x = 7*sin(7.32*t)./(1 + cos(1.42*t).^2);y = 7*cos(1.42*t).*sin(7.32*t).^4;for i = 0:pi/3:2*pi m = [cos(i) -sin(i); sin(i) cos(i)]*([x; y]); patch([m(1,:), NaN], [m(2,:), NaN], - [sqrt(m(1,:).^2+m(2,:).^2), NaN],... 'EdgeColor','interp', 'Marker','none', 'MarkerFaceColor','flat', 'LineWidth',1)endaxis equal offcolormap(bone) 樱花树
一个用循环一个用的递归：
% author : slandarerhold on; axis([-1,5,0,5])set(gca, 'XColor','none', 'YColor','none', 'Color',[.5,.5,.5])T=[1.2; 0; pi/2]; a = pi/10;for i = 1:16 L = .6*.9^i; I = randi(25, [1,size(T,2)]) > 9; if i == 1,I = ~1; end L1 = T(:,I); t = L1(3,:); R1 = L1 + [cos(t - a)*L; sin(t - a)*L; t*0 - a]; R2 = L1 + [cos(t + a)*L; sin(t + a)*L; t*0 + a]; L2 = T(:,~I); t = L2(3,:); R3 = L2 + [cos(t)*L; sin(t)*L; t*0]; T = [R1, R2, R3]; X = [L1(1,:), L1(1,:), L2(1,:); R1(1,:), R2(1,:), R3(1,:)]; X(end+1,:) = nan; Y = [L1(2,:), L1(2,:), L2(2,:); R1(2,:), R2(2,:), R3(2,:)]; Y(end+1,:) = nan; plot(X(:), Y(:), 'Color',[0 0 0] + i*.3/16, 'LineWidth',5*0.8^i) if i > 14 scatter(T(1,:),T(2,:), i*2-20, 'CData',[.86,.68,.68]/(1 - .13*(i - 15))); endend
% 原代码出自：Sebastian Kraemer% 注释改编自：slandarerhold on;% 绘制背景C = [123,159,190;228,214,200]./255;CN = size(C,1);CM = linspace(1,0,100).'*ones(1,100);M(:,:,1) = interp1(linspace(0,1,CN),C(:,1),CM);M(:,:,2) = interp1(linspace(0,1,CN),C(:,2),CM);M(:,:,3) = interp1(linspace(0,1,CN),C(:,3),CM);image([.5,14],[3.2,13],M)E = F([0,1]);% 绘制山fill([-6,7,10,12,15,28],[.2,5.5,7.8,7.8,5.5,.2],[.2,5.5,7.8,7.8,5.5,.2], 'EdgeColor','none');% 绘制树fill(real(E),imag(E),2, 'EdgeColor','none');scatter(real(E),imag(E),30,'CData',[1,.7,.7],'Marker','h','MarkerEdgeAlpha',.05);colormap bone;axis off;axis([.5,14,0,13])function V = F(p)V = [];if abs(diff(p)) > 8e-3 % ^ C____D % | /\ / % |/ \ / % A----B % | | % | | % 0----1 % 假设0，1节点为原节点，则乘虚数i是为了得到与01向量方向垂直的向量 % 原始节点再加上该向量便得到A,B节点， % 后面的1i^cos(...)^.5是为了得到C,D节点 % cos(...)的范围为[-1, 1], ^cos(...)是为了将1i旋转至[-pi, pi]范围， % 再^.5将其旋转范围调整至[-pi/2, pi/2]范围，-1i再在前面乘i*diff(p) % 就能接下来生成C,D，并且在向量的右半边方向 % 这样就能看似随机，实则依靠norm(p)对左右两个树枝的宽度和角度进行调整 Z = p + .5i*diff(p)*[4; 4 - 1i + 1i^cos(283*norm(p))^.5]; % 使用平行四边形的AC边和CB边作为原节点再进行递归，使树进行分叉 V = [p(1); F(Z(1:2)); F(Z(2:3)); p(2)]; endend
毕达哥拉斯树
% author : slandarerclc; clearX = [0, 1, 1, 0];Y = [0, 0, 1, 1];MF = @(T) [cos(T), sin(T); -sin(T), cos(T)];M1 = MF(pi/2); M2 = MF(-pi/2); M3 = MF(pi/2.5);N = 10; CL = turbo(N);ax = gca; ax.DataAspectRatio = [1,1,1];ax.NextPlot = 'add';ax.XColor = 'none';ax.YColor = 'none';fill(X,Y, CL(1,:))for i = 2:N XM = X(:,3)./2 + X(:,4)./2; XV = X(:,3)./2 - X(:,4)./2; YM = Y(:,3)./2 + Y(:,4)./2; YV = Y(:,3)./2 - Y(:,4)./2; XYM = [XV, YV]*M3 + [XM, YM]; XR = [XYM(:,1), X(:,3)]; XL = [X(:,4), XYM(:,1)]; YR = [XYM(:,2), Y(:,3)]; YL = [Y(:,4), XYM(:,2)]; XYR3 = [XR(:,1) - XR(:,2), YR(:,1) - YR(:,2)]*M2 + [XR(:,2), YR(:,2)]; XYR4 = [XR(:,2) - XR(:,1), YR(:,2) - YR(:,1)]*M1 + [XR(:,1), YR(:,1)]; XYL3 = [XL(:,1) - XL(:,2), YL(:,1) - YL(:,2)]*M2 + [XL(:,2), YL(:,2)]; XYL4 = [XL(:,2) - XL(:,1), YL(:,2) - YL(:,1)]*M1 + [XL(:,1), YL(:,1)]; X = [XR, XYR3(:,1), XYR4(:,1); XL, XYL3(:,1), XYL4(:,1)]; Y = [YR, XYR3(:,2), XYR4(:,2); YL, XYL3(:,2), XYL4(:,2)]; fill(X.',Y.', CL(i,:))end
% author : slandarerE = F([0,1], 13);fill(real(E),imag(E),imag(E), 'EdgeColor','none');hold on; axis equal offcolormap(cool)function V = F(p, n)T = 36;V = [];if n > 0 Z = p + 1i*diff(p)*[1; 1 + cosd(T)*cosd(T-90) + cosd(T)*sind(T-90)*1i]; V = [p(1); F(Z(1:2), n - 1); F(Z(2:3), n - 1); p(2)];endend
彩色涟漪
% author : slandarerhold on; axis tight equalXY = rand(200,2).*[2,1.2];n = 100; s = 2^n;r = randi(size(XY,1), n, 1);for j = 1:n scatter(XY(r,1),XY(r,2), s/2^(j-1), jet(n), 'filled', 'MarkerFaceAlpha',.2);endcamva(3)
随机油画
% author : slandarer% 复刻自 : Martin Enderclc; clearr=@(n) randi([0,n-1], [1,1]);P(:,:,1) = randMat();P(:,:,2) = randMat();P(:,:,3) = randMat();imshow(uint8(P))function M = randMat()r = @(n) randi([0,n-1],[1,1]);M = zeros(512,512);for i = 1:512 for j = 1:512 M(i,j) = getC(i,j); endend function C = getC(i,j) if M(i,j) == 0 if r(99) ~= 0 M(i,j) = getC(mod(i+r(2),512)+1,mod(j+r(2),512)+1); else M(i,j) = r(256); end end C = M(i,j); endend
绿水青山
% 原代码出自：Tim% 注释改编自：slandarerclc; clear; close all;a = 200;X = (.5:a)'./a;CL = (- cos(X.*2.*pi) + 1).^.2;r = repmat((X - .5)'.^2, a, 1) + repmat((X - .5).^2, 1, a);Z = abs(ifftn(exp(7i.*rand([a,a]))./r.^.9)).*(CL*CL').*30;% 画山surf(repmat(X, 1, a), repmat(X', a, 1), Z);m = 50;l = (m:-1:1)./m;% 绘制云彩hold onfor n = 1:m surf(repmat(X, 1, a), repmat(X', a, 1), ones([a,a]).*n, ones([a,a,3]), 'EdgeAlpha',0, 'FaceAlpha',max(.2, l(n))./2);endzlim([-a/2,a]);shading flat;CList = [140,116,78; 171,135,87; 187,170,101; 162,181,117; 95,175,146; 79,183,160; 72,162,163; 44,132,152; 36,116,141; 65,64,123; 40,26,44]./255;CList = interp1(1:size(CList, 1), CList, linspace(1, size(CList, 1), 100));colormap(CList);camva(5);axis off
星球
% 原代码出自：Tim% 注释改编自：slandarer% 先画一个bone渐变色的球[a,b,c] = sphere(99);surf(a,b,c);colormap bonehold on% 在球面外生成一些随机点% 进行三角剖分后% 设置成半透明冷色% 一些透明三角形交错叠加形成炫酷星球x = randn(3,999); x = 1.01*x./vecnorm(x);p = delaunay(x');h = patch('faces',p, 'vertices',x', 'FaceVertexCData',pink(size(p,1)), 'FaceAlpha',.25);% 设置坐标区域比例axis equal off% 设置背景色set(gcf, 'color','k')set(gcf, 'InvertHardCopy','off')% 平滑星球表面配色shading flat% 在星球外生成一些随机点当作星星r=@()rand(1,3e2);scatter(r()*10-5, r()*10-5, r().^2*200, '.w');camva(2) 画圈圈
% bird 1K = 1:9830;t = linspace(0,2*pi,200);X = @(k) (sin(pi.*k./2e4)).^12.*(cos(31.*pi.*k./1e4).^16.*sin(6.*pi.*k./1e4)./2+sin(31.*pi.*k./1e4).^20./6)... +3.*k./2e4+cos(31.*pi.*k./1e4).^6.*sin(pi./2.*(k./1e4-1).^7-pi./5);Y = @(k) -9./4.*cos(31.*pi.*k./1e4).^6.*cos(pi./2.*(k./1e4-1).^7-pi./5).*(2./3+(sin(pi.*k./2e4).*sin(3.*pi.*k./2e4)).^6)... +3./4.*cos(3.*pi.*(k-1e4)./1e5).^10.*cos(9.*pi.*(k-1e4)./1e5).^10.*cos(36.*pi.*(k-1e4)./1e5).^14+7./10.*((k-1e4)./1e4).^2;R = @(k) sin(pi.*k./2e4).^10.*(1./4.*cos(31.*pi.*k./1e4+25.*pi./32).^20+1./20.*cos(31.*pi.*k./1e4).^2)+1./30.*(3./2-cos(62.*pi.*k./1e4).^2);CX = [X(K') + cos(t).*R(K'), K'.*nan]';CY = [Y(K') + sin(t).*R(K'), K'.*nan]';plot(CX(:),CY(:), 'Color',[0,0,0,.2]);set(gca, 'DataAspectRatio',[1,1,1], 'XColor','none', 'YColor','none');
% bird 2K = -2e4:2e4;t = linspace(0,2*pi,200);X = @(k) k./15e3+sin(17.*pi./20.*(k./2e4).^5).*cos(41.*pi.*k./2e4).^6+... (1./3.*cos(41.*pi.*k./2e4).^16+1./3.*cos(41.*pi.*k./2e4).^80).*cos(pi.*k./4e4).^12.*sin(6.*pi.*k./2e4);Y = @(k) 1./2.*(k./2e4).^4-cos(17.*pi./20.*(k./2e4).^5).*(11./10+45./20.*cos(pi.*k./4e4).^8.*cos(3.*pi.*k./4e4).^6).*cos(41.*pi.*k./2e4).^6+... 12./20.*cos(3.*pi.*k./2e5).^10.*cos(9.*pi.*k./2e5).^10.*cos(8.*pi.*k./2e5).^10;R = @(k) 1./50+1./40.*sin(41.*pi.*k./2e4).^2.*sin(9.*pi.*k./2e5).^2+1./17.*cos(41.*pi.*k./2e4).^2.*cos(pi.*k./4e4).^10;CX = [X(K') + cos(t).*R(K'), K'.*nan]';CY = [Y(K') + sin(t).*R(K'), K'.*nan]';plot(CX(:),CY(:), 'Color',[0,0,0,.2]);set(gca, 'DataAspectRatio',[1,1,1], 'XColor','none', 'YColor','none');
% butterfly 1K = 1:4e4;t = linspace(0, 2*pi, 200);X = @(k) 3./2.*cos(141.*pi.*k./4e4).^9.*(1-1./2.*sin(pi.*k./4e4)).*(1-1./4.*cos(2.*pi.*k./4e4).^30.*(1+cos(32.*pi.*k./4e4).^20)).*... (1-1./2.*sin(2.*pi.*k./4e4).^30.*sin(6.*pi.*k./4e4).^10.*(1./2+1./2.*sin(18.*pi.*k./4e4).^20));Y = @(k) cos(2.*pi.*k./4e4).*cos(141.*pi.*k./4e4).^2.*(1+1./4.*cos(pi.*k./4e4).^24.*cos(3.*pi.*k./4e4).^24.*cos(21.*pi.*k./4e4).^24);R = @(k) 1./100+1./40.*(cos(141.*pi.*k./4e4).^14+sin(141.*pi.*k./4e4).^6).*(1-cos(pi.*k./4e4).^16.*cos(3.*pi.*k./4e4).^16.*cos(12.*pi.*k./4e4).^16);CX = [X(K') + cos(t).*R(K'), K'.*nan]';CY = [Y(K') + sin(t).*R(K'), K'.*nan]';plot(CX(:),CY(:), 'Color',[0,0,0,.2]);set(gca, 'DataAspectRatio',[1,1,1], 'XColor','none', 'YColor','none');
% butterfly 2K = 1:4e4;t = linspace(0, 2*pi, 200);X = @(k) 6./5.*cos(141.*pi.*k./4e4).^9.*(1-1./2.*sin(pi.*k./4e4).^3).*(1-1./4.*cos(2.*pi.*k./4e4).^30.*(1+2./3.*cos(30.*pi.*k./4e4).^20)-... sin(2.*pi.*k./4e4).^10.*sin(6.*pi.*k./4e4).^10.*(1./5+4./5.*cos(24.*pi.*k./4e4).^20));Y = @(k) cos(2.*pi.*k./4e4).*cos(141.*pi.*k./4e4).^2.*(1+1./4.*cos(pi.*k./4e4).^24.*cos(3.*pi.*k./4e4).^24.*cos(19.*pi.*k./4e4).^24);R = @(k) 1./100+1./40.*(cos(2820.*pi.*k./4e4).^6+sin(141.*pi.*k./4e4).^2).*(1-cos(pi.*k./4e4).^16.*cos(3.*pi.*k./4e4).^16.*cos(12.*pi.*k./4e4).^16);CX = [X(K') + cos(t).*R(K'), K'.*nan]';CY = [Y(K') + sin(t).*R(K'), K'.*nan]';plot(CX(:),CY(:), 'Color',[0,0,0,.2]);set(gca, 'DataAspectRatio',[1,1,1], 'XColor','none', 'YColor','none');
柔线球
function Iris% 灵感来自 Oliver Brotherhood 的 Processing 作品 Iris% 由 slandarer 使用 MATLAB 进行复刻% 生成连续随机角度randList = (rand(30,3) - .5).*2.*3.15;thetaList = interp1(linspace(1, 3000, 30), randList, 1:3000, 'makima').';% 渐变颜色% #046380, #16193B, #35478C, #4E7AC7, #7FB2F0, #ADD5F7CList = [4 99 128; 22 25 59; 53 71 140; 78 122 199; 127 178 240; 173 213 247]./255;CList = interp1(linspace(1, 3000, size(CList, 1)), CList, 1:3000, 'linear');% 坐标区域修饰hold on; axis tight equal off% 随机椭圆半径A = rand(2,1); B = rand(2,1);for i = 1:3000 X = [0; A.*cos(thetaList([1,2],i)); cos(thetaList(3,i))]; Y = [0; B.*sin(thetaList([1,2],i)); sin(thetaList(3,i))]; coe2 = ((linspace(0, 1, 100)).^((0:3)')).*((1 - linspace(0, 1, 100)).^((3:-1:0)')); bezierPnts = ([1,3,3,1].*coe2.') * [X,Y]; plot(bezierPnts(:,1),bezierPnts(:,2), 'Color',[CList(i,:), .2], 'LineWidth',.5) % pause(.01)endend
旋转多边形
% author : slandarer% remixed from Daniel Pereira / Ghost Pentagon FlowerT = linspace(0, 360, 6).';I = 20:-1:1;fill(I.*sind(T + 18.*I), I.*cosd(T + 18.*I), I, 'edgecolor','none')colormap(bone); axis equal off; set(gcf, 'color','w');
炫酷像素画
% @author:slandarer[X,Y]=meshgrid(0:1023);P(:,:,1)=bitand(bitand(mod(X,Y),mod(Y,X)),255);P(:,:,2)=bitand(bitxor(mod(X,Y),mod(Y,X)),255);P(:,:,3)=bitand(bitor(mod(X,Y),mod(Y,X)),255);imshow(uint8(P))% [X,Y]=meshgrid(0:1023);% P(:,:,1)=mod(bitand(mod(X,Y),mod(Y,X)),255);% P(:,:,2)=mod(bitxor(mod(X,Y),mod(Y,X)),255);% P(:,:,3)=mod(bitor(mod(X,Y),mod(Y,X)),255);% imshow(uint8(P))
水晶爱心
% author : slandarerhold on; grid onaxis([-22,22, -20,20, -10,10])axis offs = pi/8; t = [0:.2:s, s:.02:pi-s, pi-s:.2:pi+s, pi+s:.02:2*pi-s, 2*pi-s:.2:2*pi];x = 16*sin(t).^3;y = 13*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);lW = 'LineWidth';plot3(x,y,x.*0, 'Color',[186,110,64]./255, lW,1)set(gca, 'DataAspectRatio',[1,1,1])for i = 1:length(t) for j = 1:6 L = rand(1)*2.5 + 1.5; A = rand(1,3) - .5; A = A./norm(A); B = [x(i), y(i), 0]; C = [A(3), A(3), -A(1)-A(2)]; P0 = [.8.*A.*L + B + C./norm(C).*.14.*L, 1]'; a = B(1); b = B(2); c = 0; u = A(1); v = A(2); w = A(3); for k = 1:4 o = k*pi/2; co = cos(o); so = sin(o); cn = 1 - co; RM = [u^2 + (v^2 + w^2)*co, u*v*cn - w*so, u*w*cn + v*so, (a*(v^2 + w^2) - u*(b*v + c*w))*cn + (b*w - c*v)*so; u*v*cn + w*so, v^2 + (u^2 + w^2)*co, v*w*cn - u*so, (b*(u^2 + w^2) - v*(a*u + c*w))*cn + (c*u - a*w)*so; u*w*cn - v*so, v*w*cn + u*so, w^2 + (u^2 + v^2)*co, (c*(u^2 + v^2) - w*(a*u + b*v))*cn + (a*v - b*u)*so]; P(:,k) = RM*P0; end F = [1,3,4; 1,4,5; 1,5,6; 1,6,3; 2,3,4; 2,4,5; 2,5,6; 2,6,3]; patch('Faces',F, 'Vertices',[B;B + A.*L; P'], 'FaceColor',[0,71,177]./255,... 'FaceAlpha',.2, 'EdgeColor',[0,63,159]./255,... 'EdgeAlpha',.25, lW,.01) endend 山水画
% python出处：童晶|《Python趣味创意编程》% MATLAB代码改写 ：slandarer% 计算基础数据a = 800; b = 600;d = (.5:a)'/a;s = (-cos(d*2*pi) + 1).^.2;f = d-.5; r = f'.^2 + f.^2;% 坐标区域修饰hold onaxis([0,a,0,b]);set(gca, 'XTick',[], 'YTick',[], 'DataAspectRatio',[1,1,1], 'Color',[.67,.7,.9])% 绘制渐变背景[X,Y] = meshgrid(1:a, 301:b);V = repmat(linspace(1, 0, 300)', [1,a]);t = [0,1]; I = @interp1;V = cat(3, I(t, [.25,.68], V), I(t, [1/3,.7], V),I(t, [.5,.9], V));surf(X,Y,X.*0, 'CData',V, 'EdgeColor','none');% 绘制云彩P = abs(ifftn(exp(3i*rand(a))./r.^.8)).*(s*s');C = zeros([a,a,3]);C(:,:,1) = 1; C(:,:,2) = .75; C(:,:,3) = .88;y = (1:a)./a.*.8 + .2;image([0,a], [400,b], C, 'AlphaData',P.*(y'));% 绘制8座山c1 = [230,25,90]; c2 = [210,70,10];for i = 1:8 H = abs(ifftn(exp(5i*rand(a))./r.^1.05)).*(s*s').*10; nh = (8-i)*30 + H(400,:); hm = ceil(max(nh)); C = zeros([hm,a,3]); tc = c1 + (c2 - c1)./8.*i; tc = hsv2rgb(tc./[360,100,100]); C(:,:,1) = tc(1); C(:,:,2) = tc(2); C(:,:,3) = tc(3); P = ones(hm,a); P((1:hm)' > nh) = nan; image([-50,850], [0,hm], C, 'AlphaData',P.*.98);end
玫瑰花球
function roseball% author : slandarers = @sin; c = @cos; f = @surf;[x,t] = meshgrid((0:24)./24, (0:.5:575)./575.*20.*pi+4*pi);p = (pi/2)*exp(-t./(8*pi));u = 1 - (1 - mod(3.6*t, 2*pi)./pi).^4./2 + s(15*t)/150;y = 2*(x.^2 - x).^2.*s(p);r = u.*(x.*s(p) + y.*c(p));h = u.*(x.*c(p) - y.*s(p)) + .35;L = [.02 .04 .39; .02 .06 .69; .01 .26 .99; .17 .69 1];p={'EdgeAlpha',0.05, 'EdgeColor','none', 'FaceColor','interp', 'CData',sH(h,L)};hold onx = r.*c(t); y = r.*s(t);f(x,y,h, p{:})f(x,y,- h, p{:})rx = pi - acos(- 1/sqrt(5));Rx = [1, 0, 0; 0, c(rx), -s(rx); 0, s(rx), c(rx)];yz = 72*pi/180;Rz = @(n) [c(yz/n), - s(yz/n), 0; s(yz/n), c(yz/n), 0; 0, 0, 1];Rz1 = Rz(1); Rz2 = Rz(2);[U,V,W] = rT(x, y, h, Rx);for k = 1:5, [U,V,W] = rT(U, V, W, Rz1); f(U, V, W, p{:}), end[U,V,W] = rT(U, V, W, Rz2);for k = 1:5, [U,V,W] = rT(U, V, W, Rz1); f(U, V, - W, p{:}), endaxis equal offview(11, -.07) function c = sH(H, cL) X = rescale(H, 0, 1); c = interp1(rescale(1:size(cL,1), 0, 1), cL, X); end function [U,V,W] = rT(X, Y, Z, R) U = X; V = Y; W = Z; for i = 1:numel(X) v = [X(i); Y(i); Z(i)]; n = R*v; U(i) = n(1); V(i) = n(2); W(i) = n(3); end endend MATLAB小游戏截图
以下游戏全部是由本人编写的MATLAB版本的小游戏的截图，不过其中有部分代码为6年前写的代码，可能会极其不优雅。。大家见谅，主要没时间去重构了，以下是小游戏效果截图：
完
链接：https://pan.baidu.com/s/1a1Xd4OQrrYPIuevRNLeNNA?pwd=slan
提取码：slan
下一篇：怎么查找阿里邮箱服务器?
