Redian新闻
>
问獭灵童或其他matlab大牛
avatar
问獭灵童或其他matlab大牛# PhotoGear - 摄影器材
P*t
1
I am seeing wall street is piling their year end bonus...
avatar
h*g
2
matlab怎么做积分变换?
比如说
f(x)={u(x)v(y)dy,对dy从0到正无穷积分
怎么求f(x0)?
avatar
g*i
3
you are right。

【在 P*****t 的大作中提到】
: I am seeing wall street is piling their year end bonus...
avatar
g*r
4
獭都成灵童啦?哪位獭活佛转世的?
avatar
h*g
5
牛爱

【在 g********r 的大作中提到】
: 獭都成灵童啦?哪位獭活佛转世的?
avatar
t*g
6
我每次看见这个,加上獭老师住山里就想起了天山童獭。

【在 g********r 的大作中提到】
: 獭都成灵童啦?哪位獭活佛转世的?
avatar
h*g
7
我看行

【在 t****g 的大作中提到】
: 我每次看见这个,加上獭老师住山里就想起了天山童獭。
avatar
d*0
8
长白山还是昆仑山?

【在 t****g 的大作中提到】
: 我每次看见这个,加上獭老师住山里就想起了天山童獭。
avatar
g*r
9
天山童獭那是老妖怪。

【在 t****g 的大作中提到】
: 我每次看见这个,加上獭老师住山里就想起了天山童獭。
avatar
g*r
10
我一直以为天山是另一座山。。。

【在 d*****0 的大作中提到】
: 长白山还是昆仑山?
avatar
c*y
11
落基山啊

【在 d*****0 的大作中提到】
: 长白山还是昆仑山?
avatar
h*g
12
天阴昆仑秦南岭
大太巫雪长武台
还有祁连横断山

【在 g********r 的大作中提到】
: 我一直以为天山是另一座山。。。
avatar
b*e
13

楼上都在打酱油,我开宣布答案是u(x0){v(y)dy,哈哈哈。

【在 h******g 的大作中提到】
: matlab怎么做积分变换?
: 比如说
: f(x)={u(x)v(y)dy,对dy从0到正无穷积分
: 怎么求f(x0)?

avatar
b*t
14
numerical还是symbolic?

【在 h******g 的大作中提到】
: matlab怎么做积分变换?
: 比如说
: f(x)={u(x)v(y)dy,对dy从0到正无穷积分
: 怎么求f(x0)?

avatar
b*t
15
。。。
我还想好心回答问题呢
楼上各位发包子吧
avatar
g*r
16
chi

【在 b*******t 的大作中提到】
: 。。。
: 我还想好心回答问题呢
: 楼上各位发包子吧

avatar
c*y
17
灵童很生气,你应该写成
u(x_0) \int v(y)dy

【在 b*****e 的大作中提到】
:
: 楼上都在打酱油,我开宣布答案是u(x0){v(y)dy,哈哈哈。

avatar
h*g
18
numerical

【在 b*******t 的大作中提到】
: numerical还是symbolic?
avatar
h*g
19
看来我问题没问好
应该是
f(x)={F(x*y)dy

【在 c********y 的大作中提到】
: 灵童很生气,你应该写成
: u(x_0) \int v(y)dy

avatar
b*e
20

那还是把x0代入,该咋积咋积。

【在 h******g 的大作中提到】
: 看来我问题没问好
: 应该是
: f(x)={F(x*y)dy

avatar
h*g
21
问题是x0要根据前面的结果得到,不是个常量
本质问题就是matlab怎么对一个二元函数中的某一元积分

【在 b*****e 的大作中提到】
:
: 那还是把x0代入,该咋积咋积。

avatar
b*e
22

您老还是求mathmatica吧。

【在 h******g 的大作中提到】
: 问题是x0要根据前面的结果得到,不是个常量
: 本质问题就是matlab怎么对一个二元函数中的某一元积分

avatar
h*g
23
不行啊,都有1000行code了,还zt到mathmatica里去?
或者您老来个matlab调用mathmatica?
顺便问一句,您老那俱乐部没通过审批,原因是没有详细描述
该怎么描述您老发句话?

【在 b*****e 的大作中提到】
:
: 您老还是求mathmatica吧。

avatar
c*y
24
numerical的话,不就是和一元积分一样嘛。
除非你想一下子拿到一堆x0下的值,一个个算太慢,可以考虑一下优化的方法

【在 h******g 的大作中提到】
: 问题是x0要根据前面的结果得到,不是个常量
: 本质问题就是matlab怎么对一个二元函数中的某一元积分

avatar
x5
25
把你的函数写出来把,说不定能直接积出来呢

【在 h******g 的大作中提到】
: 不行啊,都有1000行code了,还zt到mathmatica里去?
: 或者您老来个matlab调用mathmatica?
: 顺便问一句,您老那俱乐部没通过审批,原因是没有详细描述
: 该怎么描述您老发句话?

avatar
h*g
26
怎么积?
你的意思是
y=1:1:10000000
求F(x0*y)
然后再积分?
用dblquad不行

【在 c********y 的大作中提到】
: numerical的话,不就是和一元积分一样嘛。
: 除非你想一下子拿到一堆x0下的值,一个个算太慢,可以考虑一下优化的方法

avatar
t*g
27
我觉得洪老师是要问
if
f(x) = \int u(x,y)dy
f(x0) ?= \int u(x0, y)dy

【在 c********y 的大作中提到】
: numerical的话,不就是和一元积分一样嘛。
: 除非你想一下子拿到一堆x0下的值,一个个算太慢,可以考虑一下优化的方法

avatar
x5
28
这得看x和y是不是independent的了

【在 t****g 的大作中提到】
: 我觉得洪老师是要问
: if
: f(x) = \int u(x,y)dy
: f(x0) ?= \int u(x0, y)dy

avatar
t*g
29
估计积不出来,否则洪老师也不会郁闷了。

【在 x5 的大作中提到】
: 把你的函数写出来把,说不定能直接积出来呢
avatar
h*g
30
f(w)={exp(-iwt)*exp(-t^2/2)*t dt
f(w)={exp(-iwt)*exp(-t^2/2)*t^2 dt

【在 x5 的大作中提到】
: 把你的函数写出来把,说不定能直接积出来呢
avatar
h*g
31
跟傅立叶变换差不多,但不是

【在 t****g 的大作中提到】
: 估计积不出来,否则洪老师也不会郁闷了。
avatar
t*g
32
你是说x不是y的函数就行了?
我觉得要seperatable才行吧。 u(x,y)= v(x)w(y)才行。

【在 x5 的大作中提到】
: 这得看x和y是不是independent的了
avatar
t*g
33
原来是在搞复利也变缓?

【在 h******g 的大作中提到】
: f(w)={exp(-iwt)*exp(-t^2/2)*t dt
: f(w)={exp(-iwt)*exp(-t^2/2)*t^2 dt

avatar
v*a
34
hahahah
我正想说

【在 b*****e 的大作中提到】
:
: 您老还是求mathmatica吧。

avatar
h*g
35
差不多,但不是啊
fourier transform matlab里有函数的

【在 t****g 的大作中提到】
: 原来是在搞复利也变缓?
avatar
c*y
36
自己写一个呗,simpson就挺够用了啊

【在 h******g 的大作中提到】
: 怎么积?
: 你的意思是
: y=1:1:10000000
: 求F(x0*y)
: 然后再积分?
: 用dblquad不行

avatar
t*g
37
你这个不就是后面那个t或者t^2乘Gaussian的复利也变换么?

【在 h******g 的大作中提到】
: 差不多,但不是啊
: fourier transform matlab里有函数的

avatar
b*t
38
syms x y
u=2*x*y;
v=4+x;
F=int(u*v,y);
subs(F,3)
这种效果?

【在 h******g 的大作中提到】
: numerical
avatar
x5
39
这个可积吧
第一个,用 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

avatar
c*y
40
凹凹

【在 t****g 的大作中提到】
: 你是说x不是y的函数就行了?
: 我觉得要seperatable才行吧。 u(x,y)= v(x)w(y)才行。

avatar
v*a
41
还是獭獭简单明了

【在 b*******t 的大作中提到】
: syms x y
: u=2*x*y;
: v=4+x;
: F=int(u*v,y);
: subs(F,3)
: 这种效果?

avatar
t*g
42
数学大牛教一下我,都快忘光了
f(x) = \int u(x,y) dy
求f(x0)可以把x0代进去再做一元积分?不行吧。

【在 c********y 的大作中提到】
: 自己写一个呗,simpson就挺够用了啊
avatar
t*g
43
这是啥意思?

【在 c********y 的大作中提到】
: 凹凹
avatar
b*t
44
因为我不会手动积。。

【在 v***a 的大作中提到】
: 还是獭獭简单明了
avatar
b*t
45
贪贪卖萌

【在 t****g 的大作中提到】
: 这是啥意思?
avatar
x5
46
动物园到处是这声音。。。

【在 t****g 的大作中提到】
: 这是啥意思?
avatar
x*c
47


【在 t****g 的大作中提到】
: 这是啥意思?
avatar
d*g
48
这绝对是在调戏文科生的,一进来就绕道。
avatar
h*g
49
只有一半,从0到正无穷
用对称性搞得定么?

【在 t****g 的大作中提到】
: 你这个不就是后面那个t或者t^2乘Gaussian的复利也变换么?
avatar
h*g
50
被bs了。。。。。。

【在 d*******g 的大作中提到】
: 这绝对是在调戏文科生的,一进来就绕道。
avatar
v*a
51
你老说的正确

【在 t****g 的大作中提到】
: 你是说x不是y的函数就行了?
: 我觉得要seperatable才行吧。 u(x,y)= v(x)w(y)才行。

avatar
t*g
52
可以呀,你那个很明显一个是odd一个是even。

【在 h******g 的大作中提到】
: 只有一半,从0到正无穷
: 用对称性搞得定么?

avatar
t*g
53
擦,一帮子人居然还是我们两个外行数学靠谱。

【在 v***a 的大作中提到】
: 你老说的正确
avatar
h*g
54
贪贪数学最内行

【在 t****g 的大作中提到】
: 擦,一帮子人居然还是我们两个外行数学靠谱。
avatar
c*y
55
peng,就你们俩不靠谱。
你对y做积分,x和y如果是独立的话,x在这个积分里就是当常数处理的啊。

【在 t****g 的大作中提到】
: 擦,一帮子人居然还是我们两个外行数学靠谱。
avatar
v*a
56
贪贪•辛普森比Bart Simpson还调皮

【在 h******g 的大作中提到】
: 贪贪数学最内行
avatar
b*t
57
洪老师我有没有包子阿
avatar
h*g
58
搞出来本贴到此为止每人都有包子

【在 b*******t 的大作中提到】
: 洪老师我有没有包子阿
avatar
b*t
59
我那个方法不能用么

【在 h******g 的大作中提到】
: 搞出来本贴到此为止每人都有包子
avatar
c*y
60
我明白你们俩不靠谱的原因了,你们和二元积分搞混了。
如果求
\int \int f(x,y) dx dy, 如果f(x,y)可以分离成u(x)v(y)的话,能把这个二元积分
拆成俩一元积分的乘积 \int u(x)dx * \int v(y) dy

【在 t****g 的大作中提到】
: 擦,一帮子人居然还是我们两个外行数学靠谱。
avatar
h*g
61
你那个我看不懂,要先到matlab里试一下

【在 b*******t 的大作中提到】
: 我那个方法不能用么
avatar
b*t
62
前几行是算symbolic的indefinite integral,int也可以specify上下限可以用inf,
sub那句是带入求值

【在 h******g 的大作中提到】
: 你那个我看不懂,要先到matlab里试一下
avatar
v*a
63
我来告诉你答案
>> 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里试一下
avatar
h*g
64
如果u(x,y)=v(x*y),即使x和y独立也不行啊

【在 c********y 的大作中提到】
: peng,就你们俩不靠谱。
: 你对y做积分,x和y如果是独立的话,x在这个积分里就是当常数处理的啊。

avatar
b*t
65
甜甜太体贴了

【在 v***a 的大作中提到】
: 我来告诉你答案
: >> syms x y
: >> u=2*x*y;
: >> v=4+x;
: >> F=int(u*v,y);
: >> subs(F,3)
: ans =
: 21*y^2

avatar
h*g
66
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

avatar
x5
67
红老师,你太让我失望了。。。

【在 h******g 的大作中提到】
: 如果u(x,y)=v(x*y),即使x和y独立也不行啊
avatar
b*t
68
你装的时候没装symbolic toolbox吧


【在 h******g 的大作中提到】
: ft
: 这得是啥版本的matlab
: 我的2009a 不认syms x y

avatar
v*a
69
R2009a

【在 h******g 的大作中提到】
: ft
: 这得是啥版本的matlab
: 我的2009a 不认syms x y

avatar
c*y
70
举个不行的例子来?

【在 h******g 的大作中提到】
: 如果u(x,y)=v(x*y),即使x和y独立也不行啊
avatar
c*y
71
完了,五角兽要写不好的推荐信了。。。

【在 x5 的大作中提到】
: 红老师,你太让我失望了。。。
avatar
v*a
72
据我观察 班上诸角兽中做你学生最难混毕业。。。。

【在 x5 的大作中提到】
: 红老师,你太让我失望了。。。
avatar
h*g
73
sin(x*y)?

【在 c********y 的大作中提到】
: 举个不行的例子来?
avatar
t*g
74
你做富立叶变换,求某个频率点的值能够先把omega0代进去然后对t积分?
不行吧。

【在 c********y 的大作中提到】
: 举个不行的例子来?
avatar
C*c
75
我也在想{这货是什么呢

【在 c********y 的大作中提到】
: 灵童很生气,你应该写成
: u(x_0) \int v(y)dy

avatar
c*y
76
为什么不行?
FT的问题是,你不是只求这一个点,你要拿到所有的\omega_i的值(比如i从1到1000)
,一个一个算,得算1000次积分,很不efficient,所以需要别的方法来优化。

【在 t****g 的大作中提到】
: 你做富立叶变换,求某个频率点的值能够先把omega0代进去然后对t积分?
: 不行吧。

avatar
t*g
77
嗯,好象是可以。
发10个露怯包子吧。

【在 c********y 的大作中提到】
: 为什么不行?
: FT的问题是,你不是只求这一个点,你要拿到所有的\omega_i的值(比如i从1到1000)
: ,一个一个算,得算1000次积分,很不efficient,所以需要别的方法来优化。

avatar
h*g
78
我能排么》?

【在 t****g 的大作中提到】
: 嗯,好象是可以。
: 发10个露怯包子吧。

avatar
b*t
79
pai

【在 t****g 的大作中提到】
: 嗯,好象是可以。
: 发10个露怯包子吧。

avatar
c*y
80
re

【在 t****g 的大作中提到】
: 嗯,好象是可以。
: 发10个露怯包子吧。

avatar
b*t
81
哎 泰老我那么看好你

【在 t****g 的大作中提到】
: 嗯,好象是可以。
: 发10个露怯包子吧。

avatar
t*g
82
我觉得你这个其实最efficient的算法是直接做那两个函数的FFT,然后根据你要算得x0
在频域插值。

【在 h******g 的大作中提到】
: 我能排么》?
avatar
t*g
83
唉,獭老师我让你失望了。

【在 b*******t 的大作中提到】
: 哎 泰老我那么看好你
avatar
h*g
84
这个我想过,如果最后解决不了就这么做

x0

【在 t****g 的大作中提到】
: 我觉得你这个其实最efficient的算法是直接做那两个函数的FFT,然后根据你要算得x0
: 在频域插值。

avatar
v*a
85
我想知道dblquad有什么局限性?

【在 c********y 的大作中提到】
: 为什么不行?
: FT的问题是,你不是只求这一个点,你要拿到所有的\omega_i的值(比如i从1到1000)
: ,一个一个算,得算1000次积分,很不efficient,所以需要别的方法来优化。

avatar
t*g
86
这个是最快的办法,只不过你要算一下误差的upper bond。

【在 h******g 的大作中提到】
: 这个我想过,如果最后解决不了就这么做
:
: x0

avatar
h*g
87
这个问题是在于知不知道x的值
如果知道x的值
不管f(x,y)是什么形式,都可以把x代进去当成y的一元积分
如果不知道x的值
需要f(x,y)=u(x)*v(y),才能当成y的一元积分

【在 t****g 的大作中提到】
: 嗯,好象是可以。
: 发10个露怯包子吧。

avatar
c*y
88
google了一下,二维的simpson啊

【在 v***a 的大作中提到】
: 我想知道dblquad有什么局限性?
avatar
h*g
89
dblquad(F, x0, x0, y1, y2)=0

【在 v***a 的大作中提到】
: 我想知道dblquad有什么局限性?
avatar
v*a
90

贪贪辛普森

【在 c********y 的大作中提到】
: google了一下,二维的simpson啊
avatar
c*y
91
不知道也可以啊,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的一元积分

avatar
h*g
92
好像可以。。。。
起码sin(xy)我积出来了。。。。
//flush

【在 c********y 的大作中提到】
: 不知道也可以啊,x在那个积分里就是个常数。
: 在你evaluate 积分的时候,你不是已经知道x的值了么?
: 我一开始就问你,是只算一个x_0点,还是算所有的x点。
: 如果只算一个,就直接当成y的一元积分就完事了。
: 要算所有的x,也可以对每个x_i都算一次积分,不经济罢了,不是说数学上不对。

avatar
d*g
93
进来排包,对matlab一窍不通,mathmetica很久之前有一点点经验,画图的。
相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。