24小时热门版块排行榜    

查看: 785  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 14/700 2026-08-09 23:33 by Tide man
[基金申请] 2026国自然放榜时间 +4 布布和一二 2026-08-08 4/200 2026-08-09 22:24 by otj2008
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见