24小时热门版块排行榜    

查看: 793  |  回复: 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的回帖

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的回帖

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的回帖

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 ...

直接在theta角加上90度即可。。。
除了努力,别无它法
6楼2012-09-13 13:10:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yht_ll 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +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有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见