| 查看: 233 | 回复: 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。 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
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
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★
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
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