在MATLAB和数学面前遇到它，我们是该理性还是感性呢（转）

715211229

在MATLAB和数学面前，如果我们遇到了它，此时该理性还是感性呢？


真的有这样的函数吗，网上有人这样说：
第一个常量用于控制杯罩
后面的用于控制突起
常数用于控制突起大小
次方正负控制凸凹

这里提供MATLAB代码，给大家刷一刷下
x=0:0.01:1;
y=3.*x.*log10(x)-(1/30).*exp(-  ( (30.*x-30./exp(1)).^4));
plot(y,x)
grid on
axis equal
复制代码
没有MATLAB的兄弟可以到Google绘图测试下，什么不知道怎么使用Google绘图，哦，直接在Google搜索框中输入函数表达式回车就行了，比如y=sin(x)。
  

也要有些人说二维的不算什么，求3D效果图！Okay，现在就提供三维的Mathematica代码
Plot3D[11 x*Log10[x] y (y - 1) +
  Exp[-((25 x - 25/E)^2 + (25 y - 25/2)^2)^3]/25, {x, 0, 1}, {y, 0,
  1}, Mesh -> None, BoxRatios -> Automatic, ViewPoint -> {1, 1, 1},
ViewVertical -> {1, 0, 0}]

Plot3D[Piecewise[{{12.5 (x)*Log10[x]*y*(y - 1) +
     Exp[-((25*(x) - 25/E)^2 + (25*y - 25/2)^2)^3]/25,
    x <= 0.8556}, {12.5 (x - 0.8)*Log10[x - 0.8]*y*(y - 1) +
     Exp[-((25*(x - 0.8) - 25/E)^2 + (25*y - 25/2)^2)^3]/25,
    x > 0.8556}}], {x, 0, 1.8}, {y, 0, 1}, Mesh -> All,
BoxRatios -> Automatic, ViewPoint -> {1, 1, 1},
ViewVertical -> {1, 1, 1}, Filling -> Axis, Filling -> Bottom,
ColorFunction -> "SunsetColors"]


MATLAB的代码
[X, Y] = meshgrid(0.01:0.01:1, 0.01:0.01:1);
Z = 11*X.*log10(X).*Y.*(Y-1)+exp(-((25 ...
            *X - 25/exp(1)).^2+(25*Y-25/2).^2).^3)./25;
surfl(X, Y, Z);
shading flat
colormap(pink)
view([20 30])

效果图片


有了三维函数，学习有限元的朋友可以进行静力和振动模态分析了，不妨共享下共享您的分析结果，比如什么力度比较合适哦！接下来提供罩杯型号，胸围与胸下围的差距
AA 约7.5cm
A 约10cm
B 约12.5cm
C 约15cm
D 约17.5cm
E 约20cm
F 约22.5cm
G 约25cm
H 约27.5cm
I 约30cm

请各位请发挥您的聪明才智测试并公布各杯参数，以后网友至此有对照依据！

cyrilsy

学以致用，真是太好了，谢谢