24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2165  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268分085602化学工程调剂 +4 月照花林。 2026-04-09 4/200 2026-04-09 12:48 by yutian743
[考研] 086004 求调剂 309 +6 Yin DY 2026-04-08 6/300 2026-04-09 12:27 by wang_dand
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 求调剂 +15 熊二想上岸 2026-04-04 15/750 2026-04-08 22:52 by 猪会飞
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +4 江淮北月 2026-04-05 4/200 2026-04-06 22:44 by chenzhimin
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 301求调剂 +6 细胞相关蛋白 2026-04-02 10/500 2026-04-06 08:34 by jp9609
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 283求调剂 +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
信息提示
请填处理意见