WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Werte zu klein für quiver?

 

Snupyman
Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.05.2015, 12:42     Titel: Werte zu klein für quiver?
  Antworten mit Zitat      
Hallo Leute,

ich habe ein Problem mit quiver. Ich möchte eine Funktion, die in Polarkoordinaten berechnet wurde nach dR ableiten.
Code:

Radius1=10;
Radius =linspace(0,Radius1,241);
Theta = linspace(0,2*pi,241);
[Phi,R]=meshgrid(Theta,Radius);
U0 = zeros(length(Theta),length(Radius));
U0 =  ...
[Uphi,UR]=gradient(U0);
[x,y]=pol2cart(Phi,R);

figure(1)
quiver(x,y,Uphi,UR)
 


Obwohl er mir keinen Fehler anzeigt, plottet er nur die Achsen, keine Pfeile.
Die Were von U0 sind im 10^-5 Bereich und die von Uphi und UR im Bereich von 10^-6. Sind die Werte zu klein um mit quiver geplottet zu werden?
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.495
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 27.05.2015, 12:53     Titel:
  Antworten mit Zitat      
Hallo,

auf welcher Skala sind denn Uphi und UR?
Du kannst z.B. das optionale fünfte Argument nutzen, um die Pfeile zu verlängern.
Für weitere Unterstützung hilfreich wäre ein .mat-File mit x,y,Uphi,UR, oder ein anderweitig vollständig reproduzierbarer Code.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Snupyman
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.05.2015, 13:25     Titel:
  Antworten mit Zitat      
Hier mein ganzer Quellcode. Es geht dabei um die Potentialverteilung wenn ein Elektron bei (Radius1 -20Grad) und ein Proton bei (10 0Grad) sitzt. Und wie gesagt ich möchte die Potentialverteilung nach dR ableiten bzw. die Normalenableitung bilden aber das entspricht in Polarkoordinaten ja der Ableitung nach dR. Das verlängern der Pfeile hat nichts gebracht.

Code:
clear all; close all; clc;
Radius1=10;
Radius =linspace(0,Radius1,241);
Theta = linspace(0,2*pi,241);
[Phi,R]=meshgrid(Theta,Radius);
U0 = zeros(length(Theta),length(Radius));

for m=1:size(Phi,2)
for n=1:size(R,1)
U0(n,m)=5e-3/(2*pi*0.432)*log(((R(n,1)*cos(Phi(1,m))-9.2388)^2+(R(n,1)...
    *sin(Phi(1,m))+3.8268)^2)/((R(n,1)*cos(Phi(1,m))-10)^2+(R(n,1)*sin(Phi(1,m))+2.4493e-15)^2));
end
end

[Uphi,UR]=gradient(U0);
[x,y]=pol2cart(Phi,R);

figure(1)
quiver(x,y,Uphi,UR)
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 27.05.2015, 13:46     Titel:
  Antworten mit Zitat      
könnte mir vorstellen das es daran ligt das sich deine gradienten betragsmäßig um den faktor 3*10^6 unterscheiden.
Code:
test=sqrt(Uphi.^2+UR.^2)
max(test(:))/min(test(:))
 

Zitat:
Das verlängern der Pfeile hat nichts gebracht.

das ist nicht sehr aufschlussreich wenn du nicht dazu sagst warum das nix gebracht hat. in deinem code taucht der versuch nicht auf.
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Snupyman
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.05.2015, 07:47     Titel:
  Antworten mit Zitat      
Mit dem verlängern habe ich es so versucht
Code:
quiver(x,y,Uphi,UR,10000000)

und mir wurden trotzdem keine Pfeile angezeigt.
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 28.05.2015, 08:17     Titel:
  Antworten mit Zitat      
Snupyman hat Folgendes geschrieben:
Mit dem verlängern habe ich es so versucht
Code:
quiver(x,y,Uphi,UR,10000000)

und mir wurden trotzdem keine Pfeile angezeigt.

das kann ich so nicht nachvollzihen. bei mir werden da sehr viele riesige hässliche pfeile angezeigt.
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Snupyman
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.05.2015, 09:05     Titel:
  Antworten mit Zitat      
Ok, ich weiß nicht warum aber jetzt sehe ich sie auch.
Da stellt sich aber gleich das nächste Problem.
Wenn ich die Verlängerung auf einen Faktor 100 einstelle sehe ich sehr viele winzige und ein paar große Pfeile. Die großen interessieren mich aber nicht so sehr, sondern viel mehr die kleinen und die hätte ich gerne in einer größe das man die auch gut erkennen kann. Dazu müssten die großen aber verschinden, da die Matlab sozusagen die skalierung vorgeben. Kann ich die irgendwie ausblenden oder die skalierung selber bestimmen?
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 28.05.2015, 09:09     Titel:
  Antworten mit Zitat      
wenn du nur die richtung willst klannst du auch alle auf 1 normieren. wenn du die großen weghaben willst must du hallt mittels logischer indizierung und einem schwellwert die großen pfeile wegmachen.
grüße
_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Snupyman
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.05.2015, 09:35     Titel:
  Antworten mit Zitat      
Und wie normier ich die alle auf 1?
Private Nachricht senden Benutzer-Profile anzeigen
 
Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 28.05.2015, 09:54     Titel:
  Antworten mit Zitat      
Snupyman hat Folgendes geschrieben:
Und wie normier ich die alle auf 1?

mit mathe Smile

Code:

quiver(x,y,Uphi./sqrt(Uphi.^2+UR.^2) ,UR./sqrt(Uphi.^2+UR.^2) )

_________________

richtig Fragen
Private Nachricht senden Benutzer-Profile anzeigen
 
Snupyman
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 18
Anmeldedatum: 12.05.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.05.2015, 10:22     Titel:
  Antworten mit Zitat      
Ah ja. Vielen Dank, klappt jetzt alles
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 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.