版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

[求助] 如何迭代

如何应用Matlab来计算下式？

$[latex]2u_n=\sum\limits_{j=1}^{n}a_ju_{n-j}+f_n.$ [/latex]

$其中：[latex]a_1=2,a_j=j+2,(j>1),u_1=1,f_n=n^2.$ [/latex]

回复此楼

» 猜你喜欢

新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有194人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复

1楼 2012-10-24 08:03:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

pengyizhai

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 1.5
红花: 1
帖子: 2
在线: 3.8小时
虫号: 2080154
注册: 2012-10-23

【答案】应助回帖

★
感谢参与，应助指数 +1
fegg7502: 金币+1, 3ks 2012-10-24 10:43:38

似乎缺一个初值u0,否则可能初始赋值为0,

赞一下

回复此楼

2楼2012-10-24 09:10:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

引用回帖:

2楼: Originally posted by pengyizhai at 2012-10-24 09:10:09
似乎缺一个初值u0,否则可能初始赋值为0,

谢谢，写错了，u0=1,而不是u1=1.

赞一下

回复此楼

3楼2012-10-24 10:54:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

clc

clear

m=10;   %n的上限

a(1)=2;

u(1)=1;  %这里u从1开始,u(1)对应原公式u0,u(n+1)对应原公式中n

f(1)=1;

u(2)=3/2; %u(2)对应n=1,此时a1*u0=2,fn=1

for n=2:m

    f(n)=n^2;

    a(n)=n+2;

    u(n+1)=(a*flipud(u')+f(n))/2;

end

赞一下

回复此楼

showmethemoney

4楼2012-10-24 12:08:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

引用回帖:

4楼: Originally posted by csgt0 at 2012-10-24 12:08:50
clc
clear
m=10; %n的上限
a(1)=2;
u(1)=1; %这里u从1开始,u(1)对应原公式u0,u(n+1)对应原公式中n
f(1)=1;
u(2)=3/2; %u(2)对应n=1,此时a1*u0=2,fn=1
for n=2:m
f(n)=n^2;
a(n)=n+2;
u( ...

谢谢，您这个程序不太合适呀

赞一下

回复此楼

5楼2012-10-24 15:20:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

引用回帖:

5楼: Originally posted by 我要飞 at 2012-10-24 15:20:35
谢谢，您这个程序不太合适呀...

你要什么样的啊，符号形的？

赞一下

回复此楼

showmethemoney

6楼2012-10-24 15:39:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

引用回帖:

6楼: Originally posted by csgt0 at 2012-10-24 15:39:57
你要什么样的啊，符号形的？...

谢谢，我用手算了前几项得到：u0=1,u1=1.5,u2=5.5,u3=16.5,u4=42.25,....
二用您的程序得到的结果为： 1.0e+004 *

Columns 1 through 8

0.0001 0.0001 0.0006 0.0015 0.0041 0.0107 0.0272 0.0688

Columns 9 through 11

0.1741 0.4401 1.1118

赞一下

回复此楼

7楼2012-10-24 16:04:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
计算强帖: 2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

显示的问题，比如用得到结果
vpa(u,6)
ans =
[ 1.0, 1.5, 5.5, 15.5, 41.25, 106.5, 271.5, 688.375, 1741.38, 4401.0, 11118.4]
另外你确定你算得没问题？
u4时n=3,应该有u4=(a1*u3+a2*u2+a3*u1+3^2)/2=(11+6+5+9)/2=15.5

赞一下

回复此楼

showmethemoney

8楼2012-10-24 16:14:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

我要飞

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 1515.3
散金: 144
红花: 5
帖子: 458
在线: 591.2小时
虫号: 671362
注册: 2008-12-09
专业: 计算数学与科学工程计算

引用回帖:

8楼: Originally posted by csgt0 at 2012-10-24 16:14:29
显示的问题，比如用得到结果
vpa(u,6)
ans =

另外你确定你算得没问题？
u4时n=3,应该有u4=(a1*u3+a2*u2+a3*u1+3^2)/2=(11+6+5+9)/2=15.5...

非常感谢了，如何设置可以显示出最后的结果为：[ 1.0, 1.5, 5.5, 15.5, 41.25, 106.5, 271.5, 688.375, 1741.38, 4401.0, 11118.4]，另外，u4=15.5,我算错了，谢谢你的提醒！

赞一下

回复此楼

9楼2012-10-24 16:27:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖