h******g
发帖数: 11250 | 1 matlab怎么做积分变换？
比如说
f(x)={u(x)v(y)dy,对dy从0到正无穷积分
怎么求f(x0)? |
g********r
发帖数: 8017 | |
h******g
发帖数: 11250 | 3 牛爱
【在 g********r 的大作中提到】
: 獭都成灵童啦？哪位獭活佛转世的？
|
t****g
发帖数: 35582 | 4 我每次看见这个，加上獭老师住山里就想起了天山童獭。
【在 g********r 的大作中提到】
: 獭都成灵童啦？哪位獭活佛转世的？
|
h******g
发帖数: 11250 | 5 我看行
【在 t****g 的大作中提到】
: 我每次看见这个，加上獭老师住山里就想起了天山童獭。
|
d*****0
发帖数: 68029 | 6 长白山还是昆仑山？
【在 t****g 的大作中提到】
: 我每次看见这个，加上獭老师住山里就想起了天山童獭。
|
g********r
发帖数: 8017 | 7 天山童獭那是老妖怪。
【在 t****g 的大作中提到】
: 我每次看见这个，加上獭老师住山里就想起了天山童獭。
|
g********r
发帖数: 8017 | 8 我一直以为天山是另一座山。。。
【在 d*****0 的大作中提到】
: 长白山还是昆仑山？
|
c********y
发帖数: 30813 | 9 落基山啊
【在 d*****0 的大作中提到】
: 长白山还是昆仑山？
|
h******g
发帖数: 11250 | 10 天阴昆仑秦南岭
大太巫雪长武台
还有祁连横断山
【在 g********r 的大作中提到】
: 我一直以为天山是另一座山。。。
|
|
|
b*****e
发帖数: 14299 | 11
楼上都在打酱油，我开宣布答案是u(x0){v(y)dy，哈哈哈。
【在 h******g 的大作中提到】
: matlab怎么做积分变换？
: 比如说
: f(x)={u(x)v(y)dy,对dy从0到正无穷积分
: 怎么求f(x0)?
|
b*******t
发帖数: 33714 | 12 numerical还是symbolic?
【在 h******g 的大作中提到】
: matlab怎么做积分变换？
: 比如说
: f(x)={u(x)v(y)dy,对dy从0到正无穷积分
: 怎么求f(x0)?
|
b*******t
发帖数: 33714 | 13 。。。
我还想好心回答问题呢
楼上各位发包子吧 |
g********r
发帖数: 8017 | 14 chi
【在 b*******t 的大作中提到】
: 。。。
: 我还想好心回答问题呢
: 楼上各位发包子吧
|
c********y
发帖数: 30813 | 15 灵童很生气，你应该写成
u(x_0) \int v(y)dy
【在 b*****e 的大作中提到】
:
: 楼上都在打酱油，我开宣布答案是u(x0){v(y)dy，哈哈哈。
|
h******g
发帖数: 11250 | 16 numerical
【在 b*******t 的大作中提到】
: numerical还是symbolic?
|
h******g
发帖数: 11250 | 17 看来我问题没问好
应该是
f(x)={F(x*y)dy
【在 c********y 的大作中提到】
: 灵童很生气，你应该写成
: u(x_0) \int v(y)dy
|
b*****e
发帖数: 14299 | 18
那还是把x0代入，该咋积咋积。
【在 h******g 的大作中提到】
: 看来我问题没问好
: 应该是
: f(x)={F(x*y)dy
|
h******g
发帖数: 11250 | 19 问题是x0要根据前面的结果得到，不是个常量
本质问题就是matlab怎么对一个二元函数中的某一元积分
【在 b*****e 的大作中提到】
:
: 那还是把x0代入，该咋积咋积。
|
b*****e
发帖数: 14299 | 20
您老还是求mathmatica吧。
【在 h******g 的大作中提到】
: 问题是x0要根据前面的结果得到，不是个常量
: 本质问题就是matlab怎么对一个二元函数中的某一元积分
|
|
|
h******g
发帖数: 11250 | 21 不行啊，都有1000行code了，还zt到mathmatica里去？
或者您老来个matlab调用mathmatica?
顺便问一句，您老那俱乐部没通过审批，原因是没有详细描述
该怎么描述您老发句话？
【在 b*****e 的大作中提到】
:
: 您老还是求mathmatica吧。
|
c********y
发帖数: 30813 | 22 numerical的话，不就是和一元积分一样嘛。
除非你想一下子拿到一堆x0下的值，一个个算太慢，可以考虑一下优化的方法
【在 h******g 的大作中提到】
: 问题是x0要根据前面的结果得到，不是个常量
: 本质问题就是matlab怎么对一个二元函数中的某一元积分
|
x5
发帖数: 27871 | 23 把你的函数写出来把，说不定能直接积出来呢
【在 h******g 的大作中提到】
: 不行啊，都有1000行code了，还zt到mathmatica里去？
: 或者您老来个matlab调用mathmatica?
: 顺便问一句，您老那俱乐部没通过审批，原因是没有详细描述
: 该怎么描述您老发句话？
|
h******g
发帖数: 11250 | 24 怎么积？
你的意思是
y=1:1:10000000
求F(x0*y)
然后再积分？
用dblquad不行
【在 c********y 的大作中提到】
: numerical的话，不就是和一元积分一样嘛。
: 除非你想一下子拿到一堆x0下的值，一个个算太慢，可以考虑一下优化的方法
|
t****g
发帖数: 35582 | 25 我觉得洪老师是要问
if
f(x) = \int u(x,y)dy
f(x0) ?= \int u(x0, y)dy
【在 c********y 的大作中提到】
: numerical的话，不就是和一元积分一样嘛。
: 除非你想一下子拿到一堆x0下的值，一个个算太慢，可以考虑一下优化的方法
|
x5
发帖数: 27871 | 26 这得看x和y是不是independent的了
【在 t****g 的大作中提到】
: 我觉得洪老师是要问
: if
: f(x) = \int u(x,y)dy
: f(x0) ?= \int u(x0, y)dy
|
t****g
发帖数: 35582 | 27 估计积不出来，否则洪老师也不会郁闷了。
【在 x5 的大作中提到】
: 把你的函数写出来把，说不定能直接积出来呢
|
h******g
发帖数: 11250 | 28 f(w)={exp(-iwt)*exp(-t^2/2)*t dt
f(w)={exp(-iwt)*exp(-t^2/2)*t^2 dt
【在 x5 的大作中提到】
: 把你的函数写出来把，说不定能直接积出来呢
|
h******g
发帖数: 11250 | 29 跟傅立叶变换差不多，但不是
【在 t****g 的大作中提到】
: 估计积不出来，否则洪老师也不会郁闷了。
|
t****g
发帖数: 35582 | 30 你是说x不是y的函数就行了？
我觉得要seperatable才行吧。 u(x,y）= v(x)w(y)才行。
【在 x5 的大作中提到】
: 这得看x和y是不是independent的了
|
|
|
t****g
发帖数: 35582 | 31 原来是在搞复利也变缓？
【在 h******g 的大作中提到】
: f(w)={exp(-iwt)*exp(-t^2/2)*t dt
: f(w)={exp(-iwt)*exp(-t^2/2)*t^2 dt
|
v***a
发帖数: 23651 | 32 hahahah
我正想说
【在 b*****e 的大作中提到】
:
: 您老还是求mathmatica吧。
|
h******g
发帖数: 11250 | 33 差不多，但不是啊
fourier transform matlab里有函数的
【在 t****g 的大作中提到】
: 原来是在搞复利也变缓？
|
c********y
发帖数: 30813 | 34 自己写一个呗，simpson就挺够用了啊
【在 h******g 的大作中提到】
: 怎么积？
: 你的意思是
: y=1:1:10000000
: 求F(x0*y)
: 然后再积分？
: 用dblquad不行
|
t****g
发帖数: 35582 | 35 你这个不就是后面那个t或者t^2乘Gaussian的复利也变换么？
【在 h******g 的大作中提到】
: 差不多，但不是啊
: fourier transform matlab里有函数的
|
b*******t
发帖数: 33714 | 36 syms x y
u=2*x*y;
v=4+x;
F=int(u*v,y);
subs(F,3)
这种效果？
【在 h******g 的大作中提到】
: numerical
|
x5
发帖数: 27871 | 37 这个可积吧
第一个，用 u=exp(-t^2/2)变换一下？
【在 h******g 的大作中提到】
: f(w)={exp(-iwt)*exp(-t^2/2)*t dt
: f(w)={exp(-iwt)*exp(-t^2/2)*t^2 dt
|
c********y
发帖数: 30813 | 38 凹凹
【在 t****g 的大作中提到】
: 你是说x不是y的函数就行了？
: 我觉得要seperatable才行吧。 u(x,y）= v(x)w(y)才行。
|
v***a
发帖数: 23651 | 39 还是獭獭简单明了
【在 b*******t 的大作中提到】
: syms x y
: u=2*x*y;
: v=4+x;
: F=int(u*v,y);
: subs(F,3)
: 这种效果？
|
t****g
发帖数: 35582 | 40 数学大牛教一下我，都快忘光了
f(x) = \int u(x,y) dy
求f(x0)可以把x0代进去再做一元积分？不行吧。
【在 c********y 的大作中提到】
: 自己写一个呗，simpson就挺够用了啊
|
|
|
t****g
发帖数: 35582 | 41 这是啥意思？
【在 c********y 的大作中提到】
: 凹凹
|
b*******t
发帖数: 33714 | 42 因为我不会手动积。。
【在 v***a 的大作中提到】
: 还是獭獭简单明了
|
b*******t
发帖数: 33714 | 43 贪贪卖萌
【在 t****g 的大作中提到】
: 这是啥意思？
|
x5
发帖数: 27871 | 44 动物园到处是这声音。。。
【在 t****g 的大作中提到】
: 这是啥意思？
|
x****c
发帖数: 25662 | 45 萌
【在 t****g 的大作中提到】
: 这是啥意思？
|
d*******g
发帖数: 8992 | |
h******g
发帖数: 11250 | 47 只有一半，从0到正无穷
用对称性搞得定么？
【在 t****g 的大作中提到】
: 你这个不就是后面那个t或者t^2乘Gaussian的复利也变换么？
|
h******g
发帖数: 11250 | 48 被bs了。。。。。。
【在 d*******g 的大作中提到】
: 这绝对是在调戏文科生的，一进来就绕道。
|
v***a
发帖数: 23651 | 49 你老说的正确
【在 t****g 的大作中提到】
: 你是说x不是y的函数就行了？
: 我觉得要seperatable才行吧。 u(x,y）= v(x)w(y)才行。
|
t****g
发帖数: 35582 | 50 可以呀，你那个很明显一个是odd一个是even。
【在 h******g 的大作中提到】
: 只有一半，从0到正无穷
: 用对称性搞得定么？
|
|
|
t****g
发帖数: 35582 | 51 擦，一帮子人居然还是我们两个外行数学靠谱。
【在 v***a 的大作中提到】
: 你老说的正确
|
h******g
发帖数: 11250 | 52 贪贪数学最内行
【在 t****g 的大作中提到】
: 擦，一帮子人居然还是我们两个外行数学靠谱。
|
c********y
发帖数: 30813 | 53 peng，就你们俩不靠谱。
你对y做积分，x和y如果是独立的话，x在这个积分里就是当常数处理的啊。
【在 t****g 的大作中提到】
: 擦，一帮子人居然还是我们两个外行数学靠谱。
|
v***a
发帖数: 23651 | 54 贪贪•辛普森比Bart Simpson还调皮
【在 h******g 的大作中提到】
: 贪贪数学最内行
|
b*******t
发帖数: 33714 | |
h******g
发帖数: 11250 | 56 搞出来本贴到此为止每人都有包子
【在 b*******t 的大作中提到】
: 洪老师我有没有包子阿
|
b*******t
发帖数: 33714 | 57 我那个方法不能用么
【在 h******g 的大作中提到】
: 搞出来本贴到此为止每人都有包子
|
c********y
发帖数: 30813 | 58 我明白你们俩不靠谱的原因了，你们和二元积分搞混了。
如果求
\int \int f(x,y) dx dy， 如果f(x,y)可以分离成u(x)v(y)的话，能把这个二元积分
拆成俩一元积分的乘积 \int u(x)dx * \int v(y) dy
【在 t****g 的大作中提到】
: 擦，一帮子人居然还是我们两个外行数学靠谱。
|
h******g
发帖数: 11250 | 59 你那个我看不懂，要先到matlab里试一下
【在 b*******t 的大作中提到】
: 我那个方法不能用么
|
b*******t
发帖数: 33714 | 60 前几行是算symbolic的indefinite integral，int也可以specify上下限可以用inf，
sub那句是带入求值
【在 h******g 的大作中提到】
: 你那个我看不懂，要先到matlab里试一下
|
|
|
v***a
发帖数: 23651 | 61 我来告诉你答案
>> syms x y
>> u=2*x*y;
>> v=4+x;
>> F=int(u*v,y);
>> subs(F,3)
ans =
21*y^2
【在 h******g 的大作中提到】
: 你那个我看不懂，要先到matlab里试一下
|
h******g
发帖数: 11250 | 62 如果u(x,y)=v(x*y)，即使x和y独立也不行啊
【在 c********y 的大作中提到】
: peng，就你们俩不靠谱。
: 你对y做积分，x和y如果是独立的话，x在这个积分里就是当常数处理的啊。
|
b*******t
发帖数: 33714 | 63 甜甜太体贴了
【在 v***a 的大作中提到】
: 我来告诉你答案
: >> syms x y
: >> u=2*x*y;
: >> v=4+x;
: >> F=int(u*v,y);
: >> subs(F,3)
: ans =
: 21*y^2
|
h******g
发帖数: 11250 | 64 ft
这得是啥版本的matlab
我的2009a 不认syms x y
【在 v***a 的大作中提到】
: 我来告诉你答案
: >> syms x y
: >> u=2*x*y;
: >> v=4+x;
: >> F=int(u*v,y);
: >> subs(F,3)
: ans =
: 21*y^2
|
x5
发帖数: 27871 | 65 红老师，你太让我失望了。。。
【在 h******g 的大作中提到】
: 如果u(x,y)=v(x*y)，即使x和y独立也不行啊
|
b*******t
发帖数: 33714 | 66 你装的时候没装symbolic toolbox吧
葛
【在 h******g 的大作中提到】
: ft
: 这得是啥版本的matlab
: 我的2009a 不认syms x y
|
v***a
发帖数: 23651 | 67 R2009a
【在 h******g 的大作中提到】
: ft
: 这得是啥版本的matlab
: 我的2009a 不认syms x y
|
c********y
发帖数: 30813 | 68 举个不行的例子来？
【在 h******g 的大作中提到】
: 如果u(x,y)=v(x*y)，即使x和y独立也不行啊
|
c********y
发帖数: 30813 | 69 完了，五角兽要写不好的推荐信了。。。
【在 x5 的大作中提到】
: 红老师，你太让我失望了。。。
|
v***a
发帖数: 23651 | 70 据我观察 班上诸角兽中做你学生最难混毕业。。。。
【在 x5 的大作中提到】
: 红老师，你太让我失望了。。。
|
|
|
h******g
发帖数: 11250 | 71 sin(x*y)?
【在 c********y 的大作中提到】
: 举个不行的例子来？
|
t****g
发帖数: 35582 | 72 你做富立叶变换，求某个频率点的值能够先把omega0代进去然后对t积分？
不行吧。
【在 c********y 的大作中提到】
: 举个不行的例子来？
|
C****c
发帖数: 9157 | 73 我也在想{这货是什么呢
【在 c********y 的大作中提到】
: 灵童很生气，你应该写成
: u(x_0) \int v(y)dy
|
c********y
发帖数: 30813 | 74 为什么不行？
FT的问题是，你不是只求这一个点，你要拿到所有的\omega_i的值(比如i从1到1000）
，一个一个算，得算1000次积分，很不efficient，所以需要别的方法来优化。
【在 t****g 的大作中提到】
: 你做富立叶变换，求某个频率点的值能够先把omega0代进去然后对t积分？
: 不行吧。
|
t****g
发帖数: 35582 | 75 嗯，好象是可以。
发10个露怯包子吧。
【在 c********y 的大作中提到】
: 为什么不行？
: FT的问题是，你不是只求这一个点，你要拿到所有的\omega_i的值(比如i从1到1000）
: ，一个一个算，得算1000次积分，很不efficient，所以需要别的方法来优化。
|
h******g
发帖数: 11250 | 76 我能排么》？
【在 t****g 的大作中提到】
: 嗯，好象是可以。
: 发10个露怯包子吧。
|
b*******t
发帖数: 33714 | 77 pai
【在 t****g 的大作中提到】
: 嗯，好象是可以。
: 发10个露怯包子吧。
|
c********y
发帖数: 30813 | 78 re
【在 t****g 的大作中提到】
: 嗯，好象是可以。
: 发10个露怯包子吧。
|
b*******t
发帖数: 33714 | 79 哎 泰老我那么看好你
【在 t****g 的大作中提到】
: 嗯，好象是可以。
: 发10个露怯包子吧。
|
t****g
发帖数: 35582 | 80 我觉得你这个其实最efficient的算法是直接做那两个函数的FFT，然后根据你要算得x0
在频域插值。
【在 h******g 的大作中提到】
: 我能排么》？
|
|
|
t****g
发帖数: 35582 | 81 唉，獭老师我让你失望了。
【在 b*******t 的大作中提到】
: 哎 泰老我那么看好你
|
h******g
发帖数: 11250 | 82 这个我想过，如果最后解决不了就这么做
x0
【在 t****g 的大作中提到】
: 我觉得你这个其实最efficient的算法是直接做那两个函数的FFT，然后根据你要算得x0
: 在频域插值。
|
v***a
发帖数: 23651 | 83 我想知道dblquad有什么局限性？
【在 c********y 的大作中提到】
: 为什么不行？
: FT的问题是，你不是只求这一个点，你要拿到所有的\omega_i的值(比如i从1到1000）
: ，一个一个算，得算1000次积分，很不efficient，所以需要别的方法来优化。
|
t****g
发帖数: 35582 | 84 这个是最快的办法，只不过你要算一下误差的upper bond。
【在 h******g 的大作中提到】
: 这个我想过，如果最后解决不了就这么做
:
: x0
|
h******g
发帖数: 11250 | 85 这个问题是在于知不知道x的值
如果知道x的值
不管f(x,y)是什么形式，都可以把x代进去当成y的一元积分
如果不知道x的值
需要f(x,y)=u(x)*v(y)，才能当成y的一元积分
【在 t****g 的大作中提到】
: 嗯，好象是可以。
: 发10个露怯包子吧。
|
c********y
发帖数: 30813 | 86 google了一下，二维的simpson啊
【在 v***a 的大作中提到】
: 我想知道dblquad有什么局限性？
|
h******g
发帖数: 11250 | 87 dblquad(F, x0, x0, y1, y2)=0
【在 v***a 的大作中提到】
: 我想知道dblquad有什么局限性？
|
v***a
发帖数: 23651 | 88 凹
贪贪辛普森
【在 c********y 的大作中提到】
: google了一下，二维的simpson啊
|
c********y
发帖数: 30813 | 89 不知道也可以啊，x在那个积分里就是个常数。
在你evaluate 积分的时候，你不是已经知道x的值了么？
我一开始就问你，是只算一个x_0点，还是算所有的x点。
如果只算一个，就直接当成y的一元积分就完事了。
要算所有的x，也可以对每个x_i都算一次积分，不经济罢了，不是说数学上不对。
【在 h******g 的大作中提到】
: 这个问题是在于知不知道x的值
: 如果知道x的值
: 不管f(x,y)是什么形式，都可以把x代进去当成y的一元积分
: 如果不知道x的值
: 需要f(x,y)=u(x)*v(y)，才能当成y的一元积分
|
h******g
发帖数: 11250 | 90 好像可以。。。。
起码sin(xy)我积出来了。。。。
//flush
【在 c********y 的大作中提到】
: 不知道也可以啊，x在那个积分里就是个常数。
: 在你evaluate 积分的时候，你不是已经知道x的值了么？
: 我一开始就问你，是只算一个x_0点，还是算所有的x点。
: 如果只算一个，就直接当成y的一元积分就完事了。
: 要算所有的x，也可以对每个x_i都算一次积分，不经济罢了，不是说数学上不对。
|
|
|
d*******g
发帖数: 8992 | 91 进来排包，对matlab一窍不通，mathmetica很久之前有一点点经验，画图的。 |
