24小时热门版块排行榜    

查看: 795  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yht_ll

金虫 (正式写手)

[求助] 大家帮我看看,为什么这个图多转了90度

这是我的代码。。
N=8;

theta=0:0.01:2*pi;
a0=cos((pi/2)*cos(theta));
a1=sin(N*(pi/2)*cos(theta));
a2=cos((pi/4)*(1-sin(theta)));
a3=N*sin((pi/2)*cos(theta));
a4=sin(theta);
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%
rho=abs(((a0./a4).*(a1./a3)).*a2);
figure(1);
polar(theta,rho);

画出的图如下:


理论上,我需要的最大最长的椭圆朝向应该在0度的位置,现在在90度,是哪里出现问题了,还请高手多少指教啊。。小女子不甚感激。。。谢谢先!

GIZLEHB6@{2Z~V5GWL)VE41.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

除了努力,别无它法
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yht_ll

金虫 (正式写手)

引用回帖:
4楼: Originally posted by simonyangq at 2012-09-13 10:50:18
我根据你的公式,将所有计算theta的正弦和余弦互换,得到了你期望的结果。
程序如下:
a0=cos((pi/2)*sin(theta));
a1=sin(N*(pi/2)*sin(theta));
a2=cos((pi/4)*(1-cos(theta)));
a3=N*sin((pi/2)*sin(the ...

你把SIN全都改成了COS是吧。。。结果是对了。。
那还是这个函数自己本身的问题。。。非常感谢您的帮助,谢谢!
除了努力,别无它法
5楼2012-09-13 13:06:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

luomia

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2012-09-13 18:46:23
rho=abs(((a0./a4).*(a1./a3)).*a2);
你可以查看一下,这个表达式在thta=0时的值为INF,正无穷。而在这一点周围的其它值都非常小,近似为0,所以在图上不可能出现你想要的椭圆。
一次次萍水相逢
2楼2012-09-11 15:01:48
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

simonyangq

新虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2012-09-13 18:46:34
虽然不知道你要模拟什么过程,但从函数上看,
你的a0/(a4*a3)就已经是90度最大,0度最小,而且a1*a2在0度并不能获得最大,
结局就是你画的那个图了。
建议重新考察公式中的正弦或余弦角度定义是否正确,然后在绘图
3楼2012-09-13 10:47:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

simonyangq

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yht_ll: 金币+10 2012-09-13 13:06:18
xiegangmai: 金币+2, 谢谢参与! 2012-09-13 18:46:41
我根据你的公式,将所有计算theta的正弦和余弦互换,得到了你期望的结果。
程序如下:
a0=cos((pi/2)*sin(theta));
a1=sin(N*(pi/2)*sin(theta));
a2=cos((pi/4)*(1-cos(theta)));
a3=N*sin((pi/2)*sin(theta));
a4=cos(theta);
rho=abs(((a0./a4).*(a1./a3)).*a2);
polar(theta,rho);
4楼2012-09-13 10:50:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +6 Tide man 2026-08-14 7/350 2026-08-14 18:36 by 正能量1斤
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 好奇怪的filecode +6 布布和一二 2026-08-08 7/350 2026-08-14 13:38 by Equinoxhua
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见