октава работает нормально, но сюжет не отображается? - PullRequest
Новая
       0

октава работает нормально, но сюжет не отображается?

5 голосов
/ 14 мая 2011

Привет Я использую Octave 2.3.4 с командой plot. Я новичок в Октаве. Этот сюжет не отображается по какой-то причине. Вот мой пример файла M: 

        1;

clear all;

%%%%%%%%% parameters setting  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=0.01; %risk free rate
S0=50; %underlying price 1

%create an implied volatiltiy surface using below parameters:
basevol=0.25; %implied volatility at time t=0 and in center of strike axis
skewT=-0.001; %icrease in vol for one unit increase in maturity
v1=0.1; %defines how much a smile is raised at left end from base vol
v3=0.2; %defines how much a smile is raised at right end from base vol

nK=100; %no. of strike steps
nT=10; %no. of time steps
Tmax=1; %maximum value in time axis
Kmin=1; %minimum value in strike price axis 
Kmax=150; %maximum value of strike price axis
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


dt=Tmax/(nT-1);
Tvec=0:dt:1;
dk=(Kmax-Kmin)/(nK-1);
Kvec=Kmin:dk:Kmax;
Tvec=Tvec';
Kvec=Kvec';
nK=size(Kvec,1);
nT=size(Tvec,1);
dvolT=ones(nK,nT)*(skewT*dt);
dvolT=cumsum(dvolT,2);

SmileVec=GetSmile(Kvec,v1,0,v3);
dvolK=ones(nK,nT);
dvolK=repmat(SmileVec,1,nT);

ImpliedVolSurface=ones(nK,nT)*basevol+dvolT+dvolK;

%use formula mentioned by John Elder in "Hedging for Financial Derivatives"
%this formula gives local volatility using implied volatility 
function ret=GetLocalVolSurface(ImpliedVolSurface, S, r, Kvec, Tvec)
  [m,n]=size(ImpliedVolSurface);
  LocalVolSurface=zeros(m,n);
  dk=Kvec(2)-Kvec(1);
  dt=Tvec(2)-Tvec(1);
  x=ImpliedVolSurface;
  for i=3:m-2, %loop over strikes
    for j=1:n-1, %loop over time steps
      dv_dk=(x(i+1,j)-x(i-1,j))/(2*dk);
      dv2_dk2=(x(i-1,j)-2*x(i,j)+x(i+1,j))/(dk*dk);
      dv_dt=(x(i,j+1)-x(i,j))/dt;
      K=Kvec(i);
      T=Tvec(j+1);
      rT=T^0.5;
      sig=x(i,j);
      h1=(log(S/K)+r*T+0.5*sig*sig*T)/(sig*rT);
      numer=sig*sig + 2*T*sig*dv_dt + 2*r*K*T*sig*dv_dk;
      denom=(1+K*h1*rT*dv_dk)^2 + K*K*T*sig*sig*(dv2_dk2-h1*dv_dk*dv_dk*rT);
      LocalVolSurface(i,j)=(numer/denom)^0.5;
    end
  end
  ret=LocalVolSurface;
endfunction


LocalVol_Surface=GetLocalVolSurface(ImpliedVolSurface,S0,r,Kvec,Tvec);

AsyImplVols=zeros(nK,1);
T=Tvec(nT-1);
F=S0*exp(r*T);

for i=3:nK-2,
    % use formula sigBS(F,K)=sigLoc( (F+K)/2 )
    K=Kvec(i);
    lookupK=(F+K)/2;
    kdiff=abs(Kvec-lookupK); %try to find nearest point in grid
    kidx=min(find(kdiff==min(kdiff)));
    if ( (kidx > 3) && (kidx < nK-2) ),
        AsyImplVols(i)=LocalVol_Surface(kidx);
    else
        AsyImplVols(i) = NaN;
    end
end
figure(1);
plot(Kvec(3:nK-2),[ImpliedVolSurface(3:nK-2,nT-1) LocalVol_Surface(3:nK-2,nT-1) AsyImplVols(3:nK-2)]);

Когда я запускаю в Octave без ошибок, график никогда не отображается. Он включает в себя gnuplot 1.0.1, который я понимаю, граф? Есть что-то, чего я не делаю или отсутствует? Я также запускаю это на Windows 2003 Server. Спасибо

Ответы [ 3 ]

3 голосов
/ 29 октября 2017

Я получил ответ здесь . Octave по умолчанию использует fltk для построения графиков и т. Д., Который не работает, здесь работает gnuplot. Просто добавьте следующие строки в файл .octaverc в вашем домашнем каталоге. 

graphics_toolkit("gnuplot")

Так что при каждом запуске октавы будет установлен пакет по умолчанию для gnuplot

1 голос
/ 09 ноября 2013

Я знаю, что это старый вопрос, но поскольку вчера я столкнулся с совершенно той же ошибкой, может быть, это может помочь и кому-то еще:

Согласно Октавским вики-страницам , естьпохоже, проблема с построением графиков и библиотекой "oct2mat".Для меня проблема была решена после того, как я запустил это в командной строке октавы: 

pkg rebuild -noauto oct2mat

и перезапустил октаву.Когда вам нужно использовать «oct2mat», введите: 

pkg load oct2mat

Надеюсь, это поможет!

0 голосов
/ 26 августа 2011

У меня такая же проблема при установке на 64-битный ноутбук с Windows 7 (HP 630) с графикой Intel.Каждый раз, когда вы строите график, он ничего не может сделать, но если вы строите график снова, он появляется.Это какая-то ошибка обновления.Это раздражает, но если вы строите дважды, второй раз это работает.

Мне интересно, если это какая-то ошибка двойной буферизации, потому что она работает правильно на моем ноутбуке под управлением Windows 7 с выделенной видеокартой.

В любом случае, попробуйте построить два раза подряд, я готов поспорить, что это работает, и, пожалуйста, дайте мне знать, что такое машина и видеокарта, потому что я сообщил об этом в октавную разработку.

...