24小时热门版块排行榜    

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

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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 据悉今年马上要出结果了 +5 瞬息宇宙 2026-08-10 6/300 2026-08-10 10:38 by Tide man
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +9 Lanmanbaby 2026-08-09 17/850 2026-08-10 10:35 by Tide man
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +7 yufeiwaner 2026-08-09 7/350 2026-08-10 09:08 by zju2000
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[考博] 【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
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口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变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 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
[基金申请] 求各位大神看下 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 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 纯娱乐,不喜欢勿喷 +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
信息提示
请填处理意见