[問題] 矩陣的維度不一致問題已回收
想請問一下!由於我的式子會用到兩個simpson積分
不過第一個積分中的積分元u包含了下一個積分的積分元omega
但兩者的積分範圍不一樣,導致程式出現下列問題
??? Error using ==> mtimes
Inner matrix dimensions must agree.
以下是我的原始程式碼:
function ft = paperfun(xxx)
load T1
load S1
load R1
global T1 S1 R1
%xxx=[0.01 0.01 0.01]
kappa=xxx(1);
vol=xxx(2);
theta=xxx(3);
sigma=3.1819;
sum=0;
for j=2:5551
omega=[-500:0.5:500];
i=sqrt(-1);
u=[1:1:5551];
G1=simpson5_20091209(u,exp((1+vol^2*(omega.*sigma*exp(R1(1)*(j-u))).^2*kappa*0.5-i*theta*kappa*(omega.*sigma*exp(R1(1)*(j-u)))).^(-j/kappa)));
Q1=(0.5/pi)*simpson5_20091209(omega,exp(-i*omega.*(T1(j)-S1(j))).*exp(i*omega.*exp(R1(1)*(j-u))*(T1(1)-S1(1))).*G1);
sum=sum+log(Q1);
end
請問可有前輩有解決的方法呢?
另外想請問印象中有一個指令還是按鍵可以看是從哪一行開始就錯了
是哪個呢?還是我記錯了?
謝謝
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 163.18.31.28