Matlab- Problema con ciclo for
Salve a tutti,
Premetto che uso da pochissimo Matlab. Ho un problema con un ciclo for in quanto dovrei mostrare come la funzione seguente (lo so è lunghissima):$ KAP=(F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k^2*exp(-k*(F + c - x))*(F*exp(-x) - 1) - F*exp(-k*(F + c - x))*exp(-x) + 2*F*k*exp(-k*(F + c - x))*exp(-x) + F*k^2*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1)^2 + 2*F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1))/(k*exp(-k*(F + c - x))*(F*exp(-x) - 1) + F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - F*exp(-k*(F + c - x))*exp(-x) + F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1))$
varia al variare di k tenendo costanti le altre variabili della funzione per poi ottenere un vettore di risultati da poter rapidamente inserire in un grafico. Vi in invio il codice che sto usando (ovviamente questo codice non mi da i risultati che vorrei):
Premetto che uso da pochissimo Matlab. Ho un problema con un ciclo for in quanto dovrei mostrare come la funzione seguente (lo so è lunghissima):$ KAP=(F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k^2*exp(-k*(F + c - x))*(F*exp(-x) - 1) - F*exp(-k*(F + c - x))*exp(-x) + 2*F*k*exp(-k*(F + c - x))*exp(-x) + F*k^2*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1)^2 + 2*F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1))/(k*exp(-k*(F + c - x))*(F*exp(-x) - 1) + F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - F*exp(-k*(F + c - x))*exp(-x) + F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1))$
varia al variare di k tenendo costanti le altre variabili della funzione per poi ottenere un vettore di risultati da poter rapidamente inserire in un grafico. Vi in invio il codice che sto usando (ovviamente questo codice non mi da i risultati che vorrei):
syms F x c k KAP=(F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - k^2*exp(-k*(F + c - x))*(F*exp(-x) - 1) - F*exp(-k*(F + c - x))*exp(-x) + 2*F*k*exp(-k*(F + c - x))*exp(-x) + F*k^2*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1)^2 + 2*F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1))/(k*exp(-k*(F + c - x))*(F*exp(-x) - 1) + F*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F)) - F*exp(-k*(F + c - x))*exp(-x) + F*k*exp(-x)*exp(-k*(F + c - x - (exp(F) - exp(x))/F))*(exp(x)/F - 1)) KAP=subs(KAP,{F,x,c,k},{5,4,4,0.6}) KAP=vpa(KAP) k=0.1:0.1:0.6 for i=1:length(k) KAP(i)=subs(KAP,{F,x,c,k(i)},{5,4,4}) KAP(i)=vpa(KAP(i)) end
Risposte
grazie x la risposta ma sono riuscito a risolvere!
"Pier24":
qualcuno può aiutarmi?
che errore ti dà? Dove e come hai definito le function subs e vpa?
qualcuno può aiutarmi?