| 查看: 257 | 回复: 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。 |
回复此楼» 猜你喜欢
面上本子正文33页，违规吗？会被低分嘛？
已经有12人回复
-
Sci. Bull. 悲剧经验
已经有7人回复
-
提交了我也来说说感想
已经有5人回复
-
今年审到国自然15份，谈谈感受
已经有20人回复
-
国自然评分
已经有4人回复
-
评审有感
已经有22人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 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
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 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
5