24小时热门版块排行榜    

查看: 2463  |  回复: 10
本帖产生 1 个 博学EPI ,点击这里进行查看

zhaoyanqiao

铁虫 (正式写手)

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

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

» 猜你喜欢

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

向大家学习
已阅   回复此楼   关注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的回帖

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

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见