Wo liegt mein Fehler?

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

beijo

Forum-Anfänger


	Beiträge: 14

	Anmeldedatum: 14.05.09

	Wohnort: ---

	Version: ---

Verfasst am: 24.06.2009, 16:01 Titel: Wo liegt mein Fehler?

Hallo, ich habe einen pseudoinversen WienerFilter geschrieben. Problem ist, dass ich ihn nicht zum laufen bekomme. Hier mal der Code:

Code:

function invers=pseudowiener(psf,B,alpha,beta,gamma);
tamh=max(size(psf));
tamg=max(size(B));
padd=2.*max(tamh,tamg)-1;
%padd=tamg;
H=psf2otf(psf);
g=B;
g=[g(:,(tamg-1)./2:tamg) zeros([tamg padd-tamg]),g(:,1:(tamg-1)./2-1)];
g=[g((tamg-1)./2:tamg,:);zeros([padd-tamg padd]);g(1:(tamg-1)./2-1,:)];
mask1=(abs(H)<=gamma);
mask2=(abs(H)>gamma);
Be=mask2./(H.*mask2+mask1)+abs(H).*mask1./(H.*gamma.*mask1+mask2);
clear mask1;
clear mask2;
%---
medmask=fspecial(3,3);
filtmed=mfilt.firdecim(medmask);
proxRuido= B-filtmed %annaeherung an das Rauschen
NSR = proxRuido/B;
%---
jan=jansep(NSR,'hamming');
Rn=jan.*proxRuido;
Rn=[Rn(:,(tamRn-1)./2:tamRn) zeros([tamRn padd-tamRn]),Rn(:,1:(tamRn-1)./2-1)];
Rn=[Rn((tamRn-1)./2:tamRn,:);zeros([padd-tamRn padd]);Rn(1:(tamRn-1)./2-1)];
Sn=fftn(Rn,padd,padd)./(tamRn.^2); clear Rn;
Sn=abs(Sn);
jan=jansep(tamRi,'hamming');
Ri=jan.*Ri;
Ri=[Ri(:,(tamRi-1)./2:tamRi) zeros([tamRi padd-tamRi]),Ri(:,1:(tamRi-1)./2-1)];
Ri=[Ri((tamRi-1)./2:tamRi,:);zeros([padd-tamRi padd]);Ri(1:(tamRi-1)./2-1,:)];
Si=fftn(Ri,padd,padd)./(tamRi.^2); clear Ri;
Si=abs(Si);
g=g-(proxRuido+NSR);
G=fftn(g,padd,paddd)./(filtmed^2); clear g
Red=1./(1+alpha.*(Sn./Si).*Be);
Red=Red.^beta;
W=Red.*Be;
Fs=W.*G;
fs=abs(iffts(Fs));
invers=fs+filtmed;

Funktion ohne Link?

Code jansep

Code:

function JAN=jansep(N,X);
if strcmp(X,'triang');
Xx=triang(N);
elseif strcmp (X,'hamming')
Xx=hamming(N);
elseif strcmp (X,'hanning')
Xx=hanning(N);
elseif strcmp(X,'blackman')
Xx=blackman(N)
end;

y=Xx;
Xx=ones(size(y))*Xx';
y=Xx';
jan=Xx.*y;
clear Xx; clear y;

Funktion ohne Link?

Sieht vll jemand meinen Fehler. Ich sitze schon ne Weile vor dem Programm und finde ihn einfach nicht.

beim aufrufen des Programms wird Fehlermeldung gezeigt:

Warning: Integer operands are required for colon operator when used
as index
> In pseudowiener at 8
Warning: Integer operands are required for colon operator when used
as index
> In pseudowiener at 8
Warning: Integer operands are required for colon operator when used
as index
> In pseudowiener at 9
Warning: Integer operands are required for colon operator when used
as index
> In pseudowiener at 9
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.

Error in ==> pseudowiener at 9
g=[g((tamg-1)./2:tamg,: );zeros([padd-tamg
padd]);g(1:(tamg-1)./2-1,: )];

Ich brauche dringend Hilfe!

GruB,

Philip

steve

Ehrenmitglied


	Beiträge: 2.022

	Anmeldedatum: 03.09.07

	Wohnort: Wien

	Version: R2023b

Verfasst am: 24.06.2009, 16:29 Titel:

Moin,

du versuchst in Zeile 8 und 9 auf nicht ganzzahlige Indizies zuzugreifen. (1:x/2:5), da x evtl auch ungerade sein kann, ist das Intervall nicht ganzzahlig und du bekommst diesen Fehler.

Gruß
Alex
_________________

>> I told me to.

____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?

ojieB

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 24.06.2009, 20:30 Titel:

steve hat Folgendes geschrieben:

Hallo Alex, danke für die Hilfe. Hab meine Indizies durchgeschaut. an denen liegts nicht...

Hat vll noch jemand eine idee

Andy386

Forum-Guru


	Beiträge: 485

	Anmeldedatum: 24.06.09

	Wohnort: ---

	Version: 7.1/8

Verfasst am: 25.06.2009, 14:20 Titel:

also, wär das C, würd ich mich über fehlende Absätze beschweren Smile

g=[g((tamg-1)./2:tamg,: );zeros([padd-tamg
padd]);g(1:(tamg-1)./2-1,: )];

ist

g=[
g( (tamg-1)./2:tamg , : );
zeros([padd-tamg padd]);
g( 1:(tamg-1)./2-1 , : )
];

solange padd nicht gleich size(g,2) ist kommt es zu einem Fehler - liegt es daran ? hab k.A. was das psf ist

ojieB

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 26.06.2009, 14:41 Titel:

Hallo, danke fuer die Anzaetze!

Ich bin echt am verzweifeln, hab den Fehler noch nicht gefunden.

Hat jemand vll noch einen anderen Ansatz?

@ Andy, ist deine Frage was PSF ist oder welche Werte darin liegen?

Andy386

Forum-Guru


	Beiträge: 485

	Anmeldedatum: 24.06.09

	Wohnort: ---

	Version: 7.1/8

Verfasst am: 29.06.2009, 19:31 Titel:

Also, ich hab keine Ahnung was PSF ist. auch nicht, was du da für Variablen nimmst und was die für Dimensionen haben.
Ist auch erstmal nicht soo wild, schliesslich wirft dein Code ja Fehler aus, und ich hab angemerkt, wo der Fehler liegen könnte. Sollte er dort liegen, war mein Kommentar von wegen Struktur gerechtfertigt. Wenn nicht, dann nicht.

zitat von mir
"solange padd nicht gleich size(g,2) ist kommt es zu einem Fehler"

Zitat Code:
"tamh=max(size(psf));
tamg=max(size(B));
padd=2.*max(tamh,tamg)-1;"

macht:

Code:

padd=2.*max(max(size(psf)), max(size(B)) )-1;

Funktion ohne Link?

also wird da wohl der fehler sein, oder ? Rolling Eyes


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 |

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.