|
tobikohl |
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 26.02.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 26.02.2008, 17:27
Titel: Ableiten mit MATLAB
|
|
|
|
|
Hi,
ich bin noch nicht sehr vertraut mit Matlab und möchte gerne eine Funktion EK nach einer Variablen ableiten. Mit "diff" bin ich nicht weitergekommen, was wahrscheinlich daran liegt, dass in der Funktion auch WK-Verteilungsfunktionen (normcdf, normpdf) vorkommen. Wenn ich jede Variable in der Funktion mit "syms" beschreibe, erhalte ich folgende Fehlermeldung:
"Input must be single or double"
Könnte mir bitte jemand weiterhelfen?
LG,
Tobi
Die Funktion beginnt wiefolgt:
[EK = V0 - [(1-s)*c*F]/r - ([s*c]/r + d*([(c/r) * ((2*((r - m )/(std^2)- 0.5)*exp(-r*TF) * normcdf(((r - m )/(std^2)- 0.5)*std*(TF^0.5)) - ...
2*((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)* ...
normcdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)*std*(TF^0.5)) ...
-2/(std * (TF^0.5)) * normpdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2) ...
* std * (TF^0.5)) + 2* exp(-r*TF)/(std*(TF^0.5)) ...
* normpdf(((r - m )/(std^2)- 0.5)*std*(TF^0.5)) + ...
(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)-((r - m )/(std^2)- 0.5)))/(r*TF) - ...
(-((2*((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2) + 2/ ...
(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)*(std^2)*TF))) * ...
normcdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)*std*(TF^0.5))...
- (2/(std*(TF^0.5))) * normpdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2) ...
*std*TF^0.5) + (((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)-((r - m )/(std^2)- 0.5)) + ...
(1/(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2) * (std^2)*TF))))- ...
(2*((r - m )/(std^2)- 0.5)*exp(-r*TF) * normcdf(((r - m )/(std^2)- 0.5)*std*(TF^0.5)) - ...
2*((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)* ...
normcdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)*std*(TF^0.5)) ...
-2/(std * (TF^0.5)) * normpdf(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2) ...
* std * (TF^0.5)) + 2* exp(-r*TF)/(std*(TF^0.5)) ...
* normpdf(((r - m )/(std^2)- 0.5)*std*(TF^0.5)) + ...
(((((((r - m )/(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2)-((r - m )/(std^2)- 0.5)))/(r*TF) - ...
(((r - m )/(std^2)- 0.5)+((((((r - m )/ ...
(std^2)- 0.5) * std^2)^2 + 2*r*(std^2))^0.5)/std^2))*((s*c)/r)]))][/code]
|
|
|
|
|
nschlange |
Ehrenmitglied
|
|
Beiträge: 1.320
|
|
|
|
Anmeldedatum: 06.09.07
|
|
|
|
Wohnort: NRW
|
|
|
|
Version: R2007b
|
|
|
|
|
|
Verfasst am: 26.02.2008, 23:05
Titel:
|
|
Hi,
stell doch mach das m-File rein, oder korrigiere den Code
und liefere ein ausführbares Beispiel, inkl. der Definition der
Konstanten und Variablen.
Dann kann man ein bisschen rumprobieren.
_________________
Viele Grüße
nschlange
"Chuck Norris ejakuliert fluessigen Stahl!"
|
|
|
tobikohl |
Themenstarter
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 26.02.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 27.02.2008, 16:01
Titel:
|
|
Hi,
Danke für die schnelle Antwort! Ich habe die Ableitung von der Funktion selbst rausbekommen. Allerdings werde ich damit nicht ganz glücklich. Beim Einsetzen der Ableitung in eine Gleichung, erhalte ich beim Auflösen folgende Fehlermeldung:
'()-indexing must appear last in an index expression'
Wäre klasse, wenn mir jemand diese Fehlermeldung erklären könnte und ggf. einen Lösungsvorschlag parat hätte!!
Vielen Dank und Grüße
Tobi
Beschreibung: |
|
Download |
Dateiname: |
stdEK.m |
Dateigröße: |
152.87 KB |
Heruntergeladen: |
1126 mal |
|
|
|
nschlange |
Ehrenmitglied
|
|
Beiträge: 1.320
|
|
|
|
Anmeldedatum: 06.09.07
|
|
|
|
Wohnort: NRW
|
|
|
|
Version: R2007b
|
|
|
|
|
|
Verfasst am: 27.02.2008, 21:46
Titel:
|
|
Oh man,
das hast Du doch nicht eingetippt, oder?
Wo kommt das her?
Dieser Fehler entsteht z.B., wenn Du so etwas versuchst:
Indizierung einer temporären Variable wird nicht unterstützt.
Eventuell handelt es sich um falsch gesetzte Klammern.
_________________
Viele Grüße
nschlange
"Chuck Norris ejakuliert fluessigen Stahl!"
|
|
|
tobikohl |
Themenstarter
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 26.02.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 28.02.2008, 22:50
Titel:
|
|
Vielen Dank für deine Antwort. Das Ganze ist eine Ableitung einer Funktion. Die Ableitung hat mir Matlab mithilfe von 'diff' so berechnet. Ich denke mal, dass die Fehlermeldung daher kommt, dass man eine Wahrscheinlichkeitsfunktion nicht so einfach ableiten kann (normcdf). Für Ratschläge bin ich natürlich weiterhin sehr dankbar...
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|