24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2135  |  回复: 10
本帖产生 1 个 博学EPI ,点击这里进行查看

zhaoyanqiao

铁虫 (正式写手)

[求助] 我想请问一个matlab编程的问题,sin(v/2)=3*v/14/pi

我想求解“sin(v/2)=3*v/14/pi;”,不知道matlab程序如何编写,恳请高人解答,并麻烦附上程序代码,谢谢了

» 猜你喜欢

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

向大家学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

luomia

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhaoyanqiao: 金币+10, 博学EPI+1, ★★★★★最佳答案, 非常感谢,我怎么就没想到迭代呢 2012-07-07 13:43:11
引用回帖:
7楼: Originally posted by zhaoyanqiao at 2012-07-07 09:40:24
v=0:0.1:10;
y1=sin(v/2);
y2=3*v/14/pi;
plot(v,y1)
hold on
plot(v,y2)
grid on

可以看到,在v5.5左右时有一个交点,也就是我所需要的值...

迭代法程序如下:
function v=diedai(tol,a,n)
v=a;
f=0;
for i=1:n
    f=14*pi*sin(v/2)/3;
    if abs(f-v)<=tol
        break;
    else
        v=f;
    end
end
end
不过迭代次数很惊人
diedai(0.0001,5,20000)

ans =

   -0.3535

>> diedai(0.0001,5,200000)

ans =

  -13.6383

>> diedai(0.0001,5,2000000)

ans =

   -5.5123

可以看到一共应有3个解,0,+/-5.5123
一次次萍水相逢
9楼2012-07-07 11:32:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

不是不想

新虫 (小有名气)

【答案】应助回帖

s=solve('sin(v/2)-3*v/14/pi=0 ')      试一下
2楼2012-07-06 14:13:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

luomia

木虫 (著名写手)

【答案】应助回帖

>> s=solve('sin(v/2)=3*v/14/pi');
>> s

s =

0
一次次萍水相逢
3楼2012-07-06 15:02:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyanqiao

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by 不是不想 at 2012-07-06 14:13:07
s=solve('sin(v/2)-3*v/14/pi=0 ')      试一下

>> s=solve('sin(v/2)-3*v/14/pi=0 ')

s =

14/3*sin(0)*pi

这个不是我想要的结果,谢谢回复
向大家学习
4楼2012-07-06 21:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyanqiao

铁虫 (正式写手)

引用回帖:
3楼: Originally posted by luomia at 2012-07-06 15:02:39
>> s=solve('sin(v/2)=3*v/14/pi');
>> s

s =

0

对不起,程序运行下来是这样的:
>> s=solve('sin(v/2)=3*v/14/pi');
>> s

s =

14/3*sin(0)*pi

>>
仍然不是我想要的结果
我需要的是除了0之外的那个值。
向大家学习
5楼2012-07-06 21:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luomia

木虫 (著名写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by zhaoyanqiao at 2012-07-06 21:20:36
对不起,程序运行下来是这样的:
>> s=solve('sin(v/2)=3*v/14/pi');
>> s

s =

14/3*sin(0)*pi

>>
仍然不是我想要的结果
我需要的是除了0之外的那个值。...

大哥,你这个问题只有0解。sin(v/2)=3*v/14/pi等价于v=14*sin(v/2)/(3*pi)。然后你可以到matlab里输入
>> v=0:50;
>> f=14*sin(v/2)/(3*pi);
>> plot(v,f)
>> grid on
>> axis equal
你就会看到不可能存在其它解。
一次次萍水相逢
6楼2012-07-07 07:52:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyanqiao

铁虫 (正式写手)

引用回帖:
6楼: Originally posted by luomia at 2012-07-07 07:52:46
大哥,你这个问题只有0解。sin(v/2)=3*v/14/pi等价于v=14*sin(v/2)/(3*pi)。然后你可以到matlab里输入
>> v=0:50;
>> f=14*sin(v/2)/(3*pi);
>> plot(v,f)
>> grid on
>> axis ...

v=0:0.1:10;
y1=sin(v/2);
y2=3*v/14/pi;
plot(v,y1)
hold on
plot(v,y2)
grid on

可以看到,在v5.5左右时有一个交点,也就是我所需要的值
向大家学习
7楼2012-07-07 09:40:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luomia

木虫 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by zhaoyanqiao at 2012-07-07 09:40:24
v=0:0.1:10;
y1=sin(v/2);
y2=3*v/14/pi;
plot(v,y1)
hold on
plot(v,y2)
grid on

可以看到,在v5.5左右时有一个交点,也就是我所需要的值...

你这个问题由于涉及到PI,需要指定计算精度。用常规迭代法很难收敛,以下程序希望能有所帮助。
>> for v=5:0.00000001:6;
f=14*pi*sin(v/2)/3;
if abs(f-v)<=0.000001
break;
end
end
>> vpa(v)

ans =

5.5122612

>> vpa(f)

ans =

5.5122622

>>

» 本帖已获得的红花(最新10朵)

一次次萍水相逢
8楼2012-07-07 11:09:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyanqiao

铁虫 (正式写手)

送鲜花一朵
引用回帖:
8楼: Originally posted by luomia at 2012-07-07 11:09:53
你这个问题由于涉及到PI,需要指定计算精度。用常规迭代法很难收敛,以下程序希望能有所帮助。
>> for v=5:0.00000001:6;
f=14*pi*sin(v/2)/3;
if abs(f-v)<=0.000001
break;
end
end
>> vp ...

非常感谢,这种编程方法我学习了
向大家学习
10楼2012-07-07 13:45:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhaoyanqiao 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 求调剂:085600材料与化工,考材科基,总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 一志愿双非085502,267分,过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666!
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
信息提示
请填处理意见