προβλημα σε script του matlab

Γενικά θέματα για τις γλώσσες προγραμματισμού που δεν καλύπτονται από τις άλλες περιοχές της κατηγορίας.

Συντονιστές: WebDev Moderators, Super-Moderators

Απάντηση
theo2008
Δημοσιεύσεις: 2
Εγγραφή: 28 Νοέμ 2008 18:43

προβλημα σε script του matlab

Δημοσίευση από theo2008 » 28 Νοέμ 2008 18:47

καλησπερα παιδια.εχο ενα προβλημα σε σκριπτ στο ματλαμπ και θα ηθελα την βοηθεια σας(οχι απαραιτητα την λυση)
Εχω τον παρακατω κωδικα στον οποιο θελω για v=-12:0.05:2 να υπολογισω το ολοκληρωμα teliki για καθε τιμη του v και μετα να υπολογισω το total παλι για καθε τιμη του v και τελος να κανω σε μια plot το v και το total.Το ολοκληρωμα μου θελω να ειναι απο 0 εως +απειρο.Παρακατω παραθετω τον κωδικα καθως και την σκεψη μου για το ολοκληρωμα.

A=12.56e-4;
L=0.44;
vdot=v/L;
G=2;
T=300;
q=1.6e-19;
k=8.62e-5;
%Edot μεταξυ 0.01-0.7
Edot=0.01;
%Ndot μεταξύ 10^8-10^15
Ndot=10^8;
%diaf μεταξύ 0.01-0.5
diaf=0.05;

paronomasths=1+exp((E-q*vdot)/(k*T));
sinartisi1=1./paronomasths;
sta8eroioroi=(G*Ndot)/(sqrt(pi/2)*diaf);
ek8eths=((E+Edot+q*vdot)/diaf).^2;
sinartisi2=exp(-2*ek8eths);
sinartisi3=sinartisi1.*sta8eroioroi;
teliki=sinartisi3.*sinartisi2;

Η σκεψη μου για το ολοκληρωμα ειναι η εξης:αρχικα πριν το paronomasths να βαλω μια syms Ε και μετα την teliki να βαλω το
%ipologismos olokliromatos
subtotal=double(int(teliki,E,0,inf))
total=q*A*L*subtotal

Δηλαδη ο κωδικας μου να γινει ως εξης:

syms E
paronomasths=1+exp((E-q*vdot)/(k*T));
sinartisi1=1./paronomasths;
sta8eroioroi=(G*Ndot)/(sqrt(pi/2)*diaf);
ek8eths=((E+Edot+q*vdot)/diaf).^2;
sinartisi2=exp(-2*ek8eths);
sinartisi3=sinartisi1.*sta8eroioroi;
teliki=sinartisi3.*sinartisi2;

%ipologismos olokliromatos
subtotal=double(int(teliki,E,0,inf))
total=q*A*L*subtotal

Αυτο παρατηρησα πως για μια τιμη του v πχ-2 δινει καποιο αποτελεσμα το οποιο το βρισκω λογικο νουμερο.Οταν αλλαζω ομως τιμη στο v πχ-3 τοτε η τιμη του ολοκληρωματος παραμενει η ιδια και συνεπως και του total.Κανω κατι λάθος?Υπαρχει μηπως κανας καλυτερος τροπος υπολογισμου ολοκληρωματος?
Ευχαριστω.

Απάντηση

Επιστροφή στο “γλώσσες προγραμματισμού - γενικά”

Μέλη σε σύνδεση

Μέλη σε αυτήν τη Δ. Συζήτηση: Δεν υπάρχουν εγγεγραμμένα μέλη και 0 επισκέπτες