【答案】应助回帖

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有9人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

换能器的ANSYS仿真已经有9人回复
关于MATLAB蒙特卡洛模拟问题！希望能帮个忙！在此先行表示感谢！已经有7人回复
Fortran数组的调用问题已经有10人回复
数组p[5][0]与数组p[5]的区别已经有3人回复
matlab中数组不能与矩阵相乘吗？！已经有5人回复
ANSYS四维TABLE建立移动三维热源，程序总是出错，求原因解释！已经有16人回复
fortran程序在SSH计算中心上运行的问题！已经有13人回复
fortran中的数组越界问题，已经有3人回复
IDL可视化工具入门与提高闫殿武著机械工业出版社已经有74人回复
matlab中返回数组中某个值的下标已经有4人回复
matlab如何得到图中点对应元素的的位置，请具体点。谢谢。已经有9人回复
求助udf的输入输出已经有16人回复
Fortran中数组的使用已经有4人回复
请教 fortran 运行错误的原因已经有13人回复
程序无法运行怎么回事已经有13人回复
【讨论】关于谭浩强的书已经有6人回复

1楼2015-06-28 22:54:24

回帖支持 ( 显示支持度最高的前 50 名 )

guo_zy_123

木虫 (正式写手)

应助: 32 (小学生)
金币: 4141
红花: 11
帖子: 442
在线: 197.5小时
虫号: 2601957
注册: 2013-08-17
性别: MM
专业: 动力学与控制

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
cyberspecter: 金币+70, ★★★很有帮助 2015-07-01 22:58:46

function [y] = y(k0)

if k0==0
y=1;
elseif k0==1
y=0;
else
S = [];
for i = [0:1:k0]
s = Y(i,S);
S = [S s];
end
y = S(k0-1);

end

end

function  [s] = Y(k0,S)
s1=0;
if k0==0
s=1;
elseif k0==1
s=0;
elseif k0==2
s=-0.65;
else
   s1=0;
   k=k0-2;
   for k2=1:k
      for k1=1:k2
         s1=-0.3/(k+1)/(k+2)*S(k1+1)*S(k2-k1+1)*S(k-k2+1)+s1;
      end
   end
s=s1-1/(k+1)/(k+2)*S(k+1);
end
end

将这两个函数贴到1个脚本里，用matlab计算就行了，解决了那个反复递归的再计算问题，将中间变量直接存储处理的

赞一下(1人)

老古董

7楼2015-07-01 11:47:11

普通回帖

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名

专家经验: +68
应助: 119 (高中生)
金币: 15500.8
散金: 186
红花: 29
帖子: 2631
在线: 589.7小时
虫号: 2666806
注册: 2013-09-20
性别: GG
专业: 信号理论与信号处理
管辖: 程序语言

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 谢谢应助 2015-07-04 06:17:06
xzhdty: 金币+1 2015-07-04 06:18:19

实际中Y(0)=1;Y(1)=0;程序里Y(1)=1;Y(2)=0;
实际中k1=0,k2=100；程序中k1=1,k2=101
这样就可以了，程序中下标也就是指的第几个数而已

学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？

2楼2015-06-30 13:25:05

cyberspecter

至尊木虫 (职业作家)

应助: 0 (幼儿园)
金币: 21399.6
散金: 500
帖子: 4420
在线: 1028.6小时
虫号: 1077529
注册: 2010-08-19
专业: 动力学与控制

引用回帖:

2楼: Originally posted by qqinly at 2015-06-30 13:25:05
实际中Y(0)=1;Y(1)=0;程序里Y(1)=1;Y(2)=0;
实际中k1=0,k2=100；程序中k1=1,k2=101
这样就可以了，程序中下标也就是指的第几个数而已

这个是知道的，但是里面有个k2-k1可能出现k(0),这就不好玩了。所以编程还是有些糊涂，能不能具体实现一下哩？

[ 发自小木虫客户端 ]

3楼2015-06-30 14:18:37

1134684571

木虫 (小有名气)

应助: 34 (小学生)
金币: 3450.5
红花: 10
帖子: 288
在线: 162.4小时
虫号: 2631662
注册: 2013-09-03
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cyberspecter: 金币+50 2015-07-01 22:59:18

function  [s] = Y(k0)
s1=0;
if k0==0
s=1;
elseif k0==1
s=0;
elseif k0==2
s=-0.65;
else
   s1=0;
   k=k0-2;
   for k1=1:k
      for k2=1:k1
         s1=-0.3/(k+1)/(k+2)*Y(k1)*Y(k2-k1)*Y(k-k2)+s1;
      end
   end
s=s1-1/(k+1)/(k+2)*Y(k);
end
试一下递归调用，此程序还有问题，你自己可再修改一下。要保存维数组的的话，重新定义一个就可以。。。。

4楼2015-06-30 21:15:38

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名

专家经验: +68
应助: 119 (高中生)
金币: 15500.8
散金: 186
红花: 29
帖子: 2631
在线: 589.7小时
虫号: 2666806
注册: 2013-09-20
性别: GG
专业: 信号理论与信号处理
管辖: 程序语言

【答案】应助回帖

引用回帖:

3楼: Originally posted by cyberspecter at 2015-06-30 14:18:37
这个是知道的，但是里面有个k2-k1可能出现k(0),这就不好玩了。所以编程还是有些糊涂，能不能具体实现一下哩？
...

呃，特例的话，可以在for循环中用if语句实现

学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？

5楼2015-06-30 23:36:30

guo_zy_123

木虫 (正式写手)

应助: 32 (小学生)
金币: 4141
红花: 11
帖子: 442
在线: 197.5小时
虫号: 2601957
注册: 2013-08-17
性别: MM
专业: 动力学与控制

【答案】应助回帖

感谢参与，应助指数 +1

老古董

6楼2015-07-01 11:45:20

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

【答案】应助回帖

感谢参与，应助指数 +1

可以先直接定义为递归函数,然后用递归函数赋值给数组

我用符号计算方式得到结果依次是这样的:

$\left\{1,0,-\frac{13}{20},0,\frac{247}{2400},0,-\frac{13819}{720000},0,\frac{246103}{57600000},0,-\frac{45642571}{51840000000},0,\frac{12297341329}{6842880<br /> 0000000}\right\}$

如果求精确解,容易看出计算量增加很快; 但如果不要精确结果,只是数值解,可以快很多;

» 本帖已获得的红花（最新10朵）

cyberspecter

ส็็็็็็็็็็็็็็็็็็็็

8楼2015-07-01 14:38:58

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

★ ★
xzhdty: 金币+2, 谢谢应助 2015-07-04 06:18:03

引用回帖:

8楼: Originally posted by cooooldog at 2015-07-01 14:38:58
可以先直接定义为递归函数,然后用递归函数赋值给数组

我用符号计算方式得到结果依次是这样的:

\left\{1,0,-\frac{13}{20},0,\frac{247}{2400},0,-\frac{13819}{720000},0,\frac{246103}{57600000},0,-\frac{4 ...

$\left\{1,0,-\frac{13}{20},0,\frac{247}{2400},0,-\frac{13819}{720000},0,\frac{246103}{57600000},0,-\frac{45642571}{51840000000},0,\frac{12297341329}{68428800000000}\right\}$

ส็็็็็็็็็็็็็็็็็็็็

9楼2015-07-01 14:40:17

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

引用回帖:

9楼: Originally posted by cooooldog at 2015-07-01 14:40:17
\left\{1,0,-\frac{13}{20},0,\frac{247}{2400},0,-\frac{13819}{720000},0,\frac{246103}{57600000},0,-\frac{45642571}{51840000000},0,\frac{12297341329}{68428800000000}\right\}...

看上去50个以后就衰减到近似于数值0了; 定义这种数列有什么意义呢?

$\left\{1,0,-\frac{13}{20},0,\frac{247}{2400},0,-\frac{13819}{720000},0,\frac{246103}{57600000},0,-\frac{45642571}{51840000000},0,\frac{12297341329}{68428800000000},0,-\frac{50813173303}{1368576000000000},0,\frac{25171458575317}{3284582400000000000},0,-\frac{1443087089010899}{913711104000000000000},0,\frac{12436741879271537731}{38193124147200000000000000},0,-\frac{89082918821543310379}{1326708523008000000000000000},0,\frac{1926405551376121266766819}{139145189893079040000000000000000},0,-\frac{2581785558773946078444548263}{904443734305013760000000000000000000},0,\frac{52265584131416030501849132921}{88799930277219532800000000000000000000},0,-\frac{54279758207944515170226725056093}{447271227764732067840000000000000000000000},0,\frac{68046316678125846378035627488107643}{2719409064809570972467200000000000000000000000},0,-\frac{157420340614314662890612736394073721791}{30511769707163386311081984000000000000000000000000},0,\frac{58424737000113334556542796492487419902507}{54921185472894095359947571200000000000000000000000000},0,-\frac{15397590807939717243622504829551977258357629}{70199260704444634614623895552000000000000000000000000000},0,\frac{28673442067015724731762786705914682261022379679}{634010165099089436835234761932800000000000000000000000000000},0,-\frac{276332469221813335801493121801075715273000959007471}{29633635116731440277678872772739072000000000000000000000000000000},0,\frac{1077994209212315826795846232448917004404185658377449149}{560668376408558850053684272860223242240000000000000000000000000000000},0,-\frac{124350755108958379365748025941838479074905770413237332629}{313671226801545086381385525627205976064000000000000000000000000000000000},0,\frac{102401504615292173007468054528688287560050554168964411577723}{1252763858544429256001292468667662068613120000000000000000000000000000000000},0,-\frac{15444862496714134143901647584296702576052169821944152986628344857}{916396762525250000764945440830394803190497280000000000000000000000000000000000000}\right\}$

ส็็็็็็็็็็็็็็็็็็็็

10楼2015-07-01 14:58:38