formaprecios.m

function Fn = formaprecios(F, Pc, Pv)
%formaprecios
%   Fn = formaprecios(F, Pc, Pv)
%
%   A�ade a unas recetas en tiempo discreto F procesos de intercambio con
%   los precios de compra Pc y de venta Pv correspondientes. 
%           |  F     0 |
%      Fn = |  I   -Pc |
%           | -I    Pv |

if nargin < 3, Pv = Pc; end
[filas, columnas, ancho] = size(F);
Fn = zeros(filas+columnas*2, columnas+1, ancho);
Fn(1:filas, 1:columnas, :) = F;
Fn(filas+1:filas+columnas, 1:columnas, 1) = eye(columnas);
Fn(filas+1:filas+columnas, columnas+1, 1) = -Pc;
Fn(filas+columnas+1:filas+columnas+columnas, 1:columnas, 1) = -eye(columnas);
Fn(filas+columnas+1:filas+columnas+columnas, columnas+1, 1) = Pv;