for-Schleife < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 09:34 Do 11.12.2014 | | Autor: | calabi |
| Aufgabe | Code mit fester Masse m:
clear;
j=0;
e = 1.602e-19;
w = 1e6;
r0 = 1e-2;
m = 170 * 1.66e-27;
c1 = (m * [mm] r0^2 [/mm] * [mm] w^2) [/mm] / (4 * e);
c2 = (m * [mm] r0^2 [/mm] * [mm] w^2) [/mm] / (2 * e);
set(gca,'FontSize',16)
[x, D] = fourdif(32,2);
for v = c2*[0:0.01:1]
u = c1*eig(2*(v/c2)*diag(cos(2*x))-D);
j=j+1;
V(j)=v;
U1(:,j)=u(1);
U2(:,j)=u(2);
end
hold on
plot(V,U2,'c','LineWidth',2);
plot(V,-U1,'r','LineWidth',2);
axis([0 80 0 12])
ylabel('U[V]')
xlabel('V[V]')
Code mit mehreren Massen:
clear;
j=0;
set(gca,'FontSize',16)
[x, D] = fourdif(32,2);
for q = 0:0.01:1
a = eig(2*q*diag(cos(2*x))-D);
j=j+1;
Q(j)=q;
A1(:,j)=a(1);
A2(:,j)=a(2);
end
e = 1.602e-19;
w = 1e6;
r0 = 1e-2;
for mass = 100:10:180
m = mass * 1.66e-27;
j=j+1;
U1(j,:)= A1 * (m * [mm] r0^2 [/mm] * [mm] w^2) [/mm] / (4 * e);
U2(j,:)= A2 * (m * [mm] r0^2 [/mm] * [mm] w^2) [/mm] / (4 * e);
V(j,:)= Q * (m * [mm] r0^2 [/mm] * [mm] w^2) [/mm] / (2 * e);
end
hold on
plot(V,U2,'c','LineWidth',2);
plot(V,-U1,'r','LineWidth',2);
axis([0 80 0 12])
ylabel('U[V]')
xlabel('V[V]') |
Hallo zusammen,
ich wollte den alten Code so ändern, dass in einem Schaubild die Kurven für mehreren Massen dargestellt werden. Leider ohne Erfolg.
1) Könnt Ihr mir bitte sagen, was ich dafür anders machen muss, damit die Kurven richtig dargestellt werden?
2) Ich möchte dass z.B. für m = 170 das Schaubild nur bis zum Schnittpunkt der roten und cyan gefärbten Kurven gezeichnet wird und nicht weiter. Diese Limitierung soll bei allen Massen m gelten. Wie kann man das einprogrammieren?
Danke!
Gruß calabi
![[a]](/images/paperclip.gif "Dateianhänge") Datei-Anhang
Dateianhänge:
Anhang Nr. 1 (Typ: m) [nicht öffentlich]
|
|
| |
|
Hallo calabi,
Du musst immer sehr sorgfältig abchecken, ob Du mit Zeilen- oder Spaltenvektoren arbeitest.
Ändere drei Zeilen, und es geht:
| 1: | U1(:,j)= A1 * (m * r0^2 * w^2) / (4 * e);
| | 2: | U2(:,j)= A2 * (m * r0^2 * w^2) / (4 * e);
| | 3: | V(:,j)= Q * (m * r0^2 * w^2) / (2 * e);
|
Gruss,
Hanspeter
|
|
|
| |
|
| Status: |
(Frage) beantwortet | | Datum: | 10:41 Do 11.12.2014 | | Autor: | calabi |
Hallo Hanspeter,
vielen Dank!
Kannst du mir auch zu meiner zweiten Frage helfen?
Ich möchte dass die Kurven bis zu den Schnittpunkten gezeichnet werden.
Gruß calabi
|
|
|
| |
|
Hab ich übersehen, sorry!
Füge folgendes vor "hold on" ein:
| 1: | [nx,ny]=size(U1);
| | 2: |
| | 3: | for k=1:nx,
| | 4: | for m=1:ny,
| | 5: | if(U2(k,m)>-U1(k,m))
| | 6: | U2(k,m)=NaN;
| | 7: | elseif(U2(k,m)<-U1(k,m))
| | 8: | U1(k,m)=NaN;
| | 9: | end
| | 10: | end
| | 11: | end |
Vielleicht geht es cooler, aber so geht es sicher ;)
|
|
|
| |
|
| Status: |
(Frage) beantwortet | | Datum: | 14:52 Do 11.12.2014 | | Autor: | calabi |
Hallo Hanspeter,
danke! Es geht fast! 
Die Kurven werden nicht bis zum Schnittpunkt gezeichnet, sondern bis kurz davor (siehe Anhang).
Was muss man ändern, damit bis zum Schnittpunkt gezeichnet wird?
Gruss calabi
[Dateianhang nicht öffentlich]
Dateianhänge:
Anhang Nr. 1 (Typ: png) [nicht öffentlich]
|
|
|
| |
|
Die V-Achse feiner auflösen. Oder kannst Du damit leben, dass alles dieselbe Farbe hat?
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 17:08 Do 11.12.2014 | | Autor: | calabi |
Hallo Hanspeter,
ich habe die Abszisse feiner aufgelöst, allerdings fehlen die "Spitzen" weiterhin. Liegt wohl nicht an der Auflösung.
Besser wäre schon unterschiedliche Farben, wenn es aber anders nicht geht, dann eben alle dieselbe Farbe. Auch nichts so schlimm.
Hoffe, dass es eine Lösung gibt.
Gruß calabi
|
|
|
| |
|
Das Problem bei Spitzen ist: wenn die Spitze nicht auf einem der diskreten Punkte liegt, die Du hast, dann sieht man sie eben nicht.
Du könntest die Punkte, wo die Spitzen sind, berechnen und den Vektoren hinzufügen, aber das gibt dann pötzlich sehr viel zu tun.
|
|
|
|
