网页
资讯
视频
图片
知道
文库
贴吧
地图
采购
进入贴吧
全吧搜索
吧内搜索
搜贴
搜人
进吧
搜标签
日
一
二
三
四
五
六
签到排名:今日本吧第
个签到,
本吧因你更精彩,明天继续来努力!
本吧签到人数:0
一键签到
成为超级会员,使用一键签到
一键签到
本月漏签
0
次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行
补签
。
连续签到:
天 累计签到:
天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
09月19日
漏签
0
天
mathematica吧
关注:
19,580
贴子:
73,267
看贴
图片
吧主推荐
游戏
11
回复贴,共
1
页
<<返回mathematica吧
>0< 加载中...
求助
求助为何函数无法正常运行
只看楼主
收藏
回复
不会追女人
还分不清
4
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
图一是我想运行的代码,后几张是论文的公式,我的最终目的是想通过这个代码绘制三维图像,就像最后一张图那样
不会追女人
还分不清
4
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
我没用plot3d 就算i从1到100都是null。我不知道哪里有问题
不会追女人
还分不清
4
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
r[i] = h + (2 i - 1)/(2 n)*Subscript[L, 2]
Sin[\[Alpha]] Cot[\[Theta]]; \[Alpha] =
10 Degree; Subscript[\[CurlyPhi], 1] = 0; Subscript[\[CurlyPhi], 2] \
= 0; Subscript[L, 2] = 65; Subscript[L, 1] = 65; Subscript[W, max] = \
-1.6; h = 1500; n = 100
W[x_, y_] :=
Sum[ (Subscript[W,
max]/(2 (1 + Tan[\[Alpha]] Cot[\[Theta]]) Erf[
Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - ((i - 1)/n) Subscript[L, 2]
Cos[\[Alpha]] Cos[Subscript[\[CurlyPhi],
1]] (1 + Tan[\[Alpha]] Cot[\[Theta]]))]) -
Erf[Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - (i/n)*Subscript[L, 2]
Cos[\[Alpha]] Cos[Subscript[\[CurlyPhi],
1]] (1 + Tan[ \[Alpha]] Cot[\[Theta]]))]) ((Erf[
Sqrt[Pi] (x + (i - 1)/n*Subscript[L, 2]
Sin [Subscript[\[CurlyPhi], 1]])/r[i]] -
Erf[Sqrt[
Pi] (x + (i - 1)/n*Subscript[L, 2]
Sin[Subscript[\[CurlyPhi],
1]] - (Subscript[L, 1] +
Subscript[L, 2] (2 i - 1)/(2 n)) (Sin [
Subscript[\[CurlyPhi], 1]] +
Cos[Subscript[\[CurlyPhi], 1]] +
Tan[Subscript[\[CurlyPhi], 2]]))/r[i]]))/2, {i, 1,
100}];
Plot3D[W[x_, y_], {x, -100, 100}, {y, -100, 100}, PlotRange -> All,
AxesLabel -> {"x", "y", "w"}]
qflwq
还未摆脱
8
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
\[Alpha] = 10 Degree;
\[Theta] = 85 Degree;
\[CurlyPhi]1 = 0;
\[CurlyPhi]2 = 0;
L2 = 65;
L1 = 65;
Wmax = -1.6;
h = 1500;
n = 100;
r[i_] := h + (2 i - 1)/(2 n) * L2 Sin[\[Alpha]] Cot[\[Theta]];
W[x_, y_] =
Total@N@Table[(Wmax/(2 (1 + Tan[\[Alpha]] Cot[\[Theta]]) Erf[
Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - ((i - 1)/
n) L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))]) -
Erf[Sqrt[Pi]/
r[i] (y -
h Cot[\[Theta]] - (i/n)*
L2 Cos[\[Alpha]] Cos[\[CurlyPhi]1] (1 +
Tan[\[Alpha]] Cot[\[Theta]]))]) ((Erf[
Sqrt[Pi] (x + (i - 1)/n*L2 Sin[\[CurlyPhi]1])/r[i]] -
Erf[Sqrt[
Pi] (x + (i - 1)/n*
L2 Sin[\[CurlyPhi]1] - (L1 +
L2 (2 i - 1)/(2 n)) (Sin[\[CurlyPhi]1] +
Cos[\[CurlyPhi]1] + Tan[\[CurlyPhi]2]))/r[i]]))/
2, {i, 1, 100}];
Plot3D[W[x, y], {x, -100, 100}, {y, -100, 100}, PlotRange -> All,
AxesLabel -> {"x", "y", "w"}]
qflwq
还未摆脱
8
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
按照论文公式重写了一遍,把参数赋值和画图放在了一块,这个应该ok。
r[i_] := (h + (2 i - 1)/(2 n) L2 Sin[\[Alpha]] Cos[\[Phi]1])/
Tan[\[Beta]];
a[i_] := i/n L1 Cos[\[Alpha]] Cos[\[Phi]1];
W[i_, y_] :=
Wmax/(2 (1 +
Tan[\[Alpha]] Cot[\[Theta]]) (Erf[
Sqrt[\[Pi]]/
r[i] (y - h Cot[\[Theta]] -
a[i - 1] (1 + Tan[\[Alpha] ] Cot[\[Theta]]))] -
Erf[Sqrt[\[Pi]]/
r[i] (y - h Cot[\[Theta]] -
a[i] (1 + Tan[\[Alpha]] Cot[\[Theta]]))]))
W[i_, x_, y_] :=
W[i, y]/2 (-Erf[
Sqrt[\[Pi]]*(
x + (i - 1)/
n L2 Sin[\[Phi]1] - (L + (L2 (2 i - 1))/(
2 n) (Sin[\[Phi]1] + Cos[\[Phi]1] Tan[\[Phi]2])))/r[i]] +
Erf[Sqrt[\[Pi]]*(x + (i - 1)/n L2 Sin[\[Phi]1])/r[i]])
\[Alpha] = 10 Degree; \[Phi]1 = 10 Degree; \[Phi]2 =
10 Degree; L2 = 65; L1 = 65; Wmax = -1.6; h = 900; n = 20; \[Theta] =
85 Degree; L = 65;
\[Beta] = 20 Degree;
Wdata[x_, y_] = Total@Table[N@W[i, x, y], {i, 100}];
Plot3D[Wdata[x, y], {x, -1000, 2000}, {y, -2000, 1000},
PlotRange -> All, AxesLabel -> {"x", "y", "W[x, y]"},
PlotLabel -> "W[x, y]的三维图"]
登录百度账号
扫二维码下载贴吧客户端
下载贴吧APP
看高清直播、视频!
贴吧页面意见反馈
违规贴吧举报反馈通道
贴吧违规信息处理公示