| 查看: 284 | 回复: 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号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
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
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
100