mit randn erzeugst Du nur Rauschen, dessen Werte normalverteilt sind.
Sowwohl das rand als auch das randn-Rauschen sind zunächst weiß (mehr oder weniger) und haben daher eine konstante Rauschleistungsdichte.
Um bandbegrenztes Rauschen zu erzeugen, probier mal folgendes aus:
Code:
x=randn(1,1000);
[B,A]=butter(8,0.2); % Filter Koeffizienten für Tiefpassfilter 8. Ordnung % Bandbreite = 0.2 * halber Abtastfrequenz
y=filter(B,A,x); % Das weiße Rauschen wird nun gefiltert.
Anstelle von tiefpaßgefiltertem Rauschen könntest Du auch bandpass-gefiltertes Rauschen erzeugen, indem Du
[B,A]=butter(8,[0.2,0.25])
verwendest.
Andere Filter sind natürlich auch möglich.
Nicht vergessen: Die spektrale Leistungsverteilung hat erst mal nichts mit der Verteilung der Amplitudenwerte zu tun.
Grüße
_________________
Andreas Geißler
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 10.02.2010, 19:33
Titel:
Danke erstmal für den Ansatz.
Das Problem ist ich muss da ne Frequenzabhängigkeit mit rein bringen.
Ich muss eine Bandbreiten abhängige Rauschleistung berechnen die
der Leistung P=k*T*B enspricht.
k=1.38*e-23(Joule/K )
T=290 K
B= Variable Bandbreite
Hast du vielleicht eine Idee,lies mich gerade duch Literatur durch hab aber noch keine Lösung gefunden.
sorry das ich jetzt antworte , hatte die letzten Tage kein Internet Zugang.
Hab mir die Datei Rauschleistung angesehen,hab aber immer noch Probleme die Leistung zu berechnen.
Wieso enspricht k*T der Varianz? annahme?
N=10000
X=sqrt(k*T)*rand(N)
fs=100000;
fgu=10000 ;fgo=15000; % Grenzfrequenzen
B=fgo-fgu % Bandbreite
[b,a]=fir1(300,[fgu fgo]/(fs/2)); %Koeffizienten Bestimmung
y=filter(b,a,X);
Wie kann ich aus y eine Leistung berechen die P=k*T*B enspricht ? B=fgo-fgu .
P=1/N*(sum(y(k).^2)) enspricht auch nicht P=k*T*B
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
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.
RSS