1599阅读
9回复

[求助]刚学有限差分1小时，弱弱的问一下

离线xbotao180

别来烦我

UID ：26430

注册：2009-03-01
登录：2022-01-17
发帖：32
等级：仿真新人

0楼发表于: 2009-06-09 00:43:36

d²y/dx_²=siny 其中 y(0)=0 y‘(1)=0 用有限差分刚开始不能自启动，就是还缺一个值，有个未知量，而采用符号计算，把未知量带到最后面利用y’(1)=

0来求，这样的化运算量非常大，大家都怎么搞的？？？？万分感谢

for i=1:1000

    if i==1
        y(i)=0;
    else if i==2
            y(i)=0.001^2*(-(cos(y(i-1)+n*sin(y(i-1)))))+2*y(i-1)-y(i-2);   %%%%在这个地方y(i-2)未知，只有通过末端位置来求，
        else
            y(i)=0.001^2*(-(cos(y(i-1)+n*sin(y(i-1)))))+2*y(i-1)-y(i ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

别来烦我别来烦我别来烦我别来烦我别来

离线vincentwl

UID ：15476

注册：2008-07-16
登录：2015-06-24
发帖：189
等级：八级仿真大师

1楼发表于: 2009-06-09 01:25:30

your question confused me……

http://www.meta-materials.com/

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

2楼发表于: 2009-06-09 09:21:15

可以贴出源代码，进行分析

数值计算最好不采用符号计算功能

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线xbotao180

别来烦我

UID ：26430

注册：2009-03-01
登录：2022-01-17
发帖：32
等级：仿真新人

3楼发表于: 2009-06-09 09:27:43

大家快来看看啊，马上要交

别来烦我别来烦我别来烦我别来烦我别来

离线cem-uestc

UID ：9061

注册：2008-03-07
登录：2019-01-05
发帖：2575
等级：荣誉管理员

4楼发表于: 2009-06-09 09:36:57

看不懂你要做什么？

欢迎光临
http://www.mwtee.com/home.php?mod=space&uid=13535

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

5楼发表于: 2009-06-09 11:00:10

回 4楼(cem-uestc) 的帖子

哈，我看明白了，不过要去吃饭了，:)

逆流而上

离线xbotao180

别来烦我

UID ：26430

注册：2009-03-01
登录：2022-01-17
发帖：32
等级：仿真新人

6楼发表于: 2009-06-09 13:03:41

看明白了，说说啊，我知道应该改成矩阵的格式，但右边不显含x,怎么办

别来烦我别来烦我别来烦我别来烦我别来

离线gwzhao

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

7楼发表于: 2009-06-09 15:26:47

回 6楼(xbotao180) 的帖子

好像我也没理解对呀。

逆流而上

离线mars982133

UID ：14427

注册：2008-06-28
登录：2012-07-30
发帖：134
等级：禁止发言

8楼发表于: 2009-06-16 09:30:48

用户被禁言,该主题自动屏蔽!

离线huadodo

UID ：3119

注册：2007-05-30
登录：2020-06-24
发帖：574
等级：积极交流六级

9楼发表于: 2009-06-16 09:48:50

一到学期末，论坛上这种作业贴就越来越多了

发帖回复