[求助] δ函数或者comb怎么编程

» 猜你喜欢

三无产品还有机会吗已经有4人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
压汞仪和BET测气凝胶孔隙率已经有4人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

1楼2013-04-16 20:46:21

ifelseend

捐助贵宾 (著名写手)

小蜜蜂

应助: 32 (小学生)
金币: 280.2
散金: 1839
红花: 17
帖子: 1589
在线: 212.9小时
虫号: 1792988
注册: 2012-05-03
专业: 原子和分子物理

【答案】应助回帖

感谢参与，应助指数 +1

你考虑一下，能不能把 $\delta$ 函数化简掉。就是继续推导推导。

always_move_forward,_step_by_step.

2楼2013-04-17 08:39:43

jyd3611

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.8
帖子: 33
在线: 10小时
虫号: 2059571
注册: 2012-10-14
性别: GG
专业: 光学

引用回帖:

2楼: Originally posted by ifelseend at 2013-04-17 08:39:43
你考虑一下，能不能把\delta函数化简掉。就是继续推导推导。

变成e指数的傅里叶变换形式吗？

3楼2013-04-17 15:45:38

ifelseend

捐助贵宾 (著名写手)

小蜜蜂

应助: 32 (小学生)
金币: 280.2
散金: 1839
红花: 17
帖子: 1589
在线: 212.9小时
虫号: 1792988
注册: 2012-05-03
专业: 原子和分子物理

引用回帖:

3楼: Originally posted by jyd3611 at 2013-04-17 15:45:38
变成e指数的傅里叶变换形式吗？...

具体的吧，说不好，只是说，你看看这样是不是可以看看可以这样做。有的时候 $\delta$ 函数在中间部分。也许你变化一下就没有了。比如时域变到频域什么的。

always_move_forward,_step_by_step.

4楼2013-04-17 21:20:42

Hades_h

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 962.8
红花: 1
帖子: 235
在线: 68.5小时
虫号: 2212832
注册: 2012-12-29
性别: GG
专业: 光学

Fortran编过，matlab不会。

坚持到底！

5楼2013-04-18 17:03:44

jyd3611

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.8
帖子: 33
在线: 10小时
虫号: 2059571
注册: 2012-10-14
性别: GG
专业: 光学

引用回帖:

5楼: Originally posted by Hades_h at 2013-04-18 17:03:44
Fortran编过，matlab不会。

能把你变过的代码贴上吗，我参考下，谢谢了！！！

6楼2013-04-18 23:15:25

Hades_h

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 962.8
红花: 1
帖子: 235
在线: 68.5小时
虫号: 2212832
注册: 2012-12-29
性别: GG
专业: 光学

★ ★ ★ ★ ★ ★
华丽的飘过: 金币+6, 3q 2013-04-26 03:52:51

引用回帖:

6楼: Originally posted by jyd3611 at 2013-04-18 23:15:25
能把你变过的代码贴上吗，我参考下，谢谢了！！！...

Program main

implicit none

Real,External :: Gamma

Write(*,*) "Gamma=",Gamma(X)    !输出Gamma积分结果

End

!子函数,计算Gamma积分

REAL Function Gamma(X)

implicit none

integer i

Real X,Y,S,T,U

Real :: A(11)=(/0.0000677106,-0.0003442342,0.0015397681,-0.0024467480,0.0109736958,-0.0002109075,0.0742379071,0.0815782188,0.4118402578,0.4227843370,1.0/)

IF (X<=0) THEN

Write(*,*) "ERR**X<0"

Gamma=-1.0

End IF

Y=X

IF (Y<=1.0) THEN

  T=1.0/(Y*(Y+1.0))

  Y=Y+2.0

  ELSE IF (Y<=2.0) THEN

  T=1.0/Y

  Y=Y+1.0

  ELSE IF (Y<=3.0) THEN

  T=1.0

  ELSE

  T=1.0

10 IF (Y>3.0) THEN

  Y=Y-1.0

  T=T*Y

  GOTO 10

  END IF

ENd IF

S=A(1)

U=Y-2.0

Do i=1,10

S=S*U+A(i+1)

END do

S=S*T

Gamma=S

Return

END

!子函数结束

!----------------------------------------------------------------

赞一下(1人)

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

jyd3611

坚持到底！

7楼2013-04-20 16:10:30

jyd3611

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.8
帖子: 33
在线: 10小时
虫号: 2059571
注册: 2012-10-14
性别: GG
专业: 光学

送红花一朵

引用回帖:

7楼: Originally posted by Hades_h at 2013-04-20 16:10:30
Program main

implicit none

Real,External :: Gamma

Write(*,*) "Gamma=",Gamma(X) !输出Gamma积分结果

End

!子函数,计算Gamma积分

REAL Function Gamma(X)

implicit none

...

谢谢！！！

8楼2013-04-20 19:52:04