| 查看: 302 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】如何编写该算法
|
|||
|
已知函数f(n),f(N)为0到n上第N点函数值，其关系为: f(1)*sin[(N-1)*pi/n]+f(2)*sin[(N-2)*pi/n]+...f(N)*sin[(N-N)*pi/n]+f(N+1)*sin[(N-(N+1))*pi/n]+...+f(n)*sin[(N-n)*pi/n]=0,求N的值 希望各位给小弟帮忙用matlab如何编写算法求N值，小弟以前没有学习过matlab，最近才开始学习，由于很急，望各位高手给予解决，小弟我刚到这里不久，没有M币答谢给各位，在此先谢谢各位的帮助！～小弟我万分感激。希望高手给出算法。 注：f(n)是一个关于n的函数，为了编程方便在这可以假设为f(n)=2*n+1。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★
lvzhu2007(金币+4,VIP+0):thanks
cffang(金币+2,VIP+0):谢谢，希望和你交流一下！
lvzhu2007(金币+4,VIP+0):thanks
cffang(金币+2,VIP+0):谢谢，希望和你交流一下！
|
>> x = fsolve(@(x)3*sin((x-1)*pi)+5*sin((x-2)*pi/2)+7*sin((x-3)*pi/3),0) Optimization terminated: first-order optimality is less than options.TolFun. x = 0 >> x = fsolve(@(x)3*sin((x-1)*pi)+5*sin((x-2)*pi/2)+7*sin((x-3)*pi/3),1) Optimization terminated: first-order optimality is less than options.TolFun. x = 4.6319 >> ezplot('3*sin((x-1)*pi)+5*sin((x-2)*pi/2)+7*sin((x-3)*pi/3)-y') 见附件 |
2楼2009-01-06 11:04:41
3楼2009-01-06 22:52:32
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★
cffang(金币+5,VIP+0):谢谢，希望和你交流一下！
cffang(金币+5,VIP+0):谢谢，希望和你交流一下！
|
function f=myfunsolve(N) n=1000; I=1:n; J=2*I+1; f=sum(J.*sin((N-I)*pi/n)); >> fsolve(@myfunsolve,0) Optimization terminated: first-order optimality is less than options.TolFun. ans = -319.3346 %----------------------------------------------------------------------- >> n=1000; I=1:n; J=2*I+1; N=[-400:.1:0]'; f=sum(repmat(J,size(N,1),1).*sin((repmat(N,1,n)-repmat(I,size(N,1),1))*pi/n),2); plot(N,f); >> hold on; >> scatter(-319.3,0,'O','filled') >> text(-319.3,0,'root of the system of nonlinear equation') >> |
4楼2009-01-07 12:54:04
50