• /
• # Spiral galaxy

on 7 Oct 2021
• 10
• 76
• 1
• 0
• 277
R=(randn(1,1e5))+2
R = 1×100000
1.1772 1.8791 1.8538 2.8144 1.7708 3.5220 0.8706 3.2366 2.1968 1.0210 1.0951 1.9913 3.7272 2.4903 1.9489 2.3572 2.3391 3.5966 1.7794 1.5240 3.2847 2.1064 2.3015 2.7089 0.7204 2.4916 0.5949 2.8161 2.5975 0.9157
e=find(abs(R)<6)
e = 1×99996
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
R=R(e)
R = 1×99996
1.1772 1.8791 1.8538 2.8144 1.7708 3.5220 0.8706 3.2366 2.1968 1.0210 1.0951 1.9913 3.7272 2.4903 1.9489 2.3572 2.3391 3.5966 1.7794 1.5240 3.2847 2.1064 2.3015 2.7089 0.7204 2.4916 0.5949 2.8161 2.5975 0.9157
Theta=2*pi*randn(1,1e5)
Theta = 1×100000
1.2563 -5.1143 -3.8810 -2.1797 -4.6443 1.2953 -6.1351 -9.0861 -6.2247 0.1733 3.4588 -0.2528 10.9207 11.5668 0.9006 -1.3313 -0.8972 -2.6600 6.5525 -16.7713 -10.0739 17.8400 6.7094 10.5526 -5.3957 -1.9533 0.8799 7.4855 -13.6546 -5.1073
X=R.*cos(Theta(e))
X = 1×99996
0.3641 0.7351 -1.3697 -1.6097 -0.1205 0.9581 0.8611 -3.0527 2.1931 1.0057 -1.0405 1.9280 -0.2789 1.3464 1.2105 0.5592 1.4592 -3.1875 1.7153 -0.7407 -2.6166 1.1212 2.0955 -1.1610 0.4548 -0.9301 0.3791 1.0144 1.2053 0.3523
Y=R.*sin(Theta(e))
Y = 1×99996
1.1195 1.7293 1.2493 -2.3086 1.7667 3.3892 0.1284 -1.0754 0.1285 0.1761 -0.3415 -0.4980 -3.7167 -2.0949 1.5273 -2.2899 -1.8282 -1.6660 0.4734 1.3319 1.9857 -1.7832 0.9516 -2.4474 0.5586 -2.3115 0.4585 2.6270 -2.3010 0.8452
Z=exp(-((R.^2))/128).*sin(-9*log(R+3)+2*angle(X+Y*1i))/sqrt(8*pi)
Z = 1×99996
0.1581 0.1158 -0.0022 -0.1831 0.1903 -0.1790 0.1250 0.0162 -0.1609 0.0751 0.0969 -0.1294 -0.1789 0.1898 -0.0048 0.1685 0.1754 0.0546 -0.1610 0.0006 0.1551 0.1606 -0.1913 -0.1423 0.1162 -0.0715 0.0647 -0.1439 0.1744 0.0965
n=.04*randn(1,numel(e))
n = 1×99996
0.0075 0.0211 -0.0197 -0.0089 -0.0032 0.0419 -0.0062 0.0170 -0.0708 -0.0470 -0.0247 -0.0195 -0.0025 0.0526 -0.0140 -0.0043 0.0683 0.0125 -0.0386 0.0996 0.0074 -0.0793 0.0691 0.0450 0.0147 -0.0035 -0.1055 -0.0207 -0.0893 -0.0186
scatter3(X+n,Y+n,Z+n,1,abs(Z+n))
axis('equal','off')
colormap jet
set(gcf,'Color','k')
view(2)