24小时热门版块排行榜    

查看: 782  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ANUsWpSnuWQD 2026-08-06 4/200 2026-08-07 04:46 by nSM3ys5fefM3
[基金申请] 听说今天filecode变了 +19 布布和一二 2026-08-06 34/1700 2026-08-07 03:13 by 虫友是什么虫
[基金申请] 基金中了 +6 laoda193707 2026-08-06 6/300 2026-08-06 23:34 by dragonxp
[基金申请] filecode +10 等待解的谜 2026-08-06 14/700 2026-08-06 21:16 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 7/350 2026-08-06 21:05 by archvillain
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见