[問題] 請問有關quadgk積分的問題已回收

看板MATLAB作者 (就是愛打球)時間14年前 (2010/05/14 03:50), 編輯推噓0(0017)
留言17則, 2人參與, 最新討論串1/1
各位前輩好: 因為小弟目前遇有關積分的問題如下 在第二個積分中,積分源是一個向量,不知是否跟迴圈有關 明明寫的是迴圈,卻只跑出第一筆數據的積分 便跑出 ??? Error using ==> mtimes Inner matrix dimensions must agree.不知是何故? 以下是小弟的程式碼,不知可否請前輩賜教,問題是在第二個積分的quadgk有問題 omega=[-1000:1000]; for w=1:2001 alphaintegral=(@(u)(-i)*omega(w)*exp(revert*(t-u))*theta*h-(h/k)*log(1+(vol^2* G1(w)=quadgk(alphaintegral,0,t); end for j=1:5458 %country6 charact = (@(omega) exp(-i*omega.*X6(j))*exp(G1)); V(j)=(0.5/pi)*quadgk(charact,-1000,1000); end 謝謝 懇請會的前輩能施予援手 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 163.18.96.5

05/14 05:22, , 1F
*? .*? exp(G1)? exp(G1(j))?
05/14 05:22, 1F

05/14 07:51, , 2F
補充一下 quadgk 吃的只能是 scalar-valued function
05/14 07:51, 2F

05/14 07:52, , 3F
若是向量可能得把各分量分開做。
05/14 07:52, 3F

05/14 11:37, , 4F
G1是第一個積分的結果 會是一個1乘2001的矩陣
05/14 11:37, 4F

05/14 11:43, , 5F
Omega 也是一個1乘2001的矩陣
05/14 11:43, 5F

05/14 11:52, , 6F
Omega 是 1*2001 的矩陣,exp(G1) 也是 1*2001 的矩陣
05/14 11:52, 6F

05/14 11:53, , 7F
你覺得他們兩個要怎麼相乘...
05/14 11:53, 7F

05/14 22:00, , 8F
可是我用點乘,不是可以向量相乘,一個元素對一個元素乘
05/14 22:00, 8F

05/15 02:19, , 9F
...)*exp(G1) 你是用乘啊,而且就算改成點乘應該也是
05/15 02:19, 9F

05/15 02:20, , 10F
不行,前面說過他只能積 "scalar-valued" function...
05/15 02:20, 10F

05/15 12:23, , 11F
阿 不小心用掉那個點...是喔 那可以請問一下
05/15 12:23, 11F

05/15 12:24, , 12F
所謂的"scalar-valued" function是什麼意思嗎?
05/15 12:24, 12F

05/15 15:18, , 13F
函數值是純量的意思
05/15 15:18, 13F

05/17 10:45, , 14F
是指說只能一次積一個值 不能積向量嗎?
05/17 10:45, 14F

05/17 14:42, , 15F
之前就說各分量要分開做啦
05/17 14:42, 15F

05/18 12:26, , 16F
不太了解各分量做的意思 是指將向量中的數一個一個取出
05/18 12:26, 16F

05/18 12:27, , 17F
做嗎? 那是用迴圈的方式嗎?
05/18 12:27, 17F
文章代碼(AID): #1Bx5Vkgo (MATLAB)