送小公主一个蝴蝶结（matlab代码实现）

目录

1 那一低头的温柔

2 蝴蝶结代码

---

1 那一低头的温柔

每一次给她系鞋带，她低头看着我，让我想起徐志摩的那一低头的温柔与娇羞，好美的花朵啊，甜到心里；我每次给她系一个蝴蝶结，她就像一个花仙子，她身边萦绕着飞舞的蝴蝶，这个春天有你，空气都是甜的。爱一旦入心，就是一辈子的真心相眷，哪怕你的世界我永远无法抵达，我的心，也只会为你一人倾情绝恋。

走过许多路，遇见许多人，可这世俗红尘中纵是有再多相遇，他人再好，也再难走进我心里，因为我的心早已被你满满占据。

你是爱的唯一，是我思念的全部，红尘万里路，牵着你之手，看花，看月，听风，听雨，但在我的灵魂世界里，却是风也是你，雨也是你，晴天是你，飞絮是你。

---

2 蝴蝶结代码

%=====欢迎关注公众号：电力系统与算法之美========
dx=pi/60;
col=0:dx:pi;
az=0:dx:2*pi;
[phi,theta]=meshgrid(az,col);
%% 计算 l=3 的网格上的 P_{l}^{m}(\cos \theta)
l=3;
Plm=legendre(l,cos(theta));
% 由于 legendre 为 m 的所有值计算答案，因此 Plm 会包含一些额外的函数值。
% 提取 m=2 的值并丢弃其余值。
% 使用 reshape 函数将结果定向为与 phi 和 theta 具有相同大小的矩阵。
m=2;
if l~=0
    Plm=reshape(Plm(m+1,:,:),size(phi));
end
%% 
a=(2*l+1)*factorial(l-m);
b=4*pi*factorial(l+m);
C=sqrt(a/b);
Ylm=C.*Plm.*exp(1i*m*phi);
%% 球面坐标转换为笛卡尔坐标并绘图
[Xm,Ym,Zm]=sph2cart(phi,pi/2-theta,abs(real(Ylm)));
surf(Xm,Ym,Zm)
title('wlz','interpreter','latex')