[求助]MATLAB求助:psdchk函数的相关帮助

pheigenbau · 发表于 2006-6-18 10:28

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

在看一些函数的代码时，经常看见一些＂内部＂命令，比如函数psd(),里面就有一个函数是psdchk()，可是却找不到这个函数的相关帮助，不知这是怎么回事．类似的例子还不少，一般情况是母函数func(),内部有一子函数为funcchk(),当作是参数检测，求各位大虾指点迷津

[此贴子已经被ericlin于2006-6-18 14:37:46编辑过]

happy · 发表于 2006-6-18 22:35

function [msg,nfft,Fs,window,noverlap,p,dflag] = psdchk(P,x,y) %PSDCHK Helper function for PSD, CSD, COHERE, and TFE. % [msg,nfft,Fs,window,noverlap,p,dflag]=PSDCHK(P,x,y) takes the cell % array P and uses each element as an input argument. Assumes P has % between 0 and 7 elements which are the arguments to psd, csd, cohere % or tfe after the x (psd) or x and y (csd, cohere, tfe) arguments. % y is optional; if given, it is checked to match the size of x. % x must be a numeric vector. % Outputs: % msg - error message, [] if no error % nfft - fft length % Fs - sampling frequency % window - window vector % noverlap - overlap of sections, in samples % p - confidence interval, [] if none desired % dflag - detrending flag, 'linear' 'mean' or 'none'
% Author(s): T. Krauss, 10-28-93
msg = [];
if length(P) == 0 % psd(x) nfft = min(length(x),256); window = hanning(nfft); noverlap = 0; Fs = 2; p = []; dflag = 'none'; elseif length(P) == 1 % psd(x,nfft) % psd(x,dflag) if isempty(P{1}), dflag = 'none'; nfft = min(length(x),256); elseif isstr(P{1}), dflag = P{1}; nfft = min(length(x),256); else dflag = 'none'; nfft = P{1}; end Fs = 2; window = hanning(nfft); noverlap = 0; p = []; elseif length(P) == 2 % psd(x,nfft,Fs) % psd(x,nfft,dflag) if isempty(P{1}), nfft = min(length(x),256); else nfft=P{1}; end if isempty(P{2}), dflag = 'none'; Fs = 2; elseif isstr(P{2}), dflag = P{2}; Fs = 2; else dflag = 'none'; Fs = P{2}; end window = hanning(nfft); noverlap = 0; p = []; elseif length(P) == 3 % psd(x,nfft,Fs,window) % psd(x,nfft,Fs,dflag) if isempty(P{1}), nfft = min(length(x),256); else nfft=P{1}; end if isempty(P{2}), Fs = 2; else Fs = P{2}; end if isstr(P{3}) dflag = P{3}; window = hanning(nfft); else dflag = 'none'; window = P{3}; if length(window) == 1, window = hanning(window); end if isempty(window), window = hanning(nfft); end end noverlap = 0; p = []; elseif length(P) == 4 % psd(x,nfft,Fs,window,noverlap) % psd(x,nfft,Fs,window,dflag) if isempty(P{1}), nfft = min(length(x),256); else nfft=P{1}; end if isempty(P{2}), Fs = 2; else Fs = P{2}; end window = P{3}; if length(window) == 1, window = hanning(window); end if isempty(window), window = hanning(nfft); end if isstr(P{4}) dflag = P{4}; noverlap = 0; else dflag = 'none'; if isempty(P{4}), noverlap = 0; else noverlap = P{4}; end end p = []; elseif length(P) == 5 % psd(x,nfft,Fs,window,noverlap,p) % psd(x,nfft,Fs,window,noverlap,dflag) if isempty(P{1}), nfft = min(length(x),256); else nfft=P{1}; end if isempty(P{2}), Fs = 2; else Fs = P{2}; end window = P{3}; if length(window) == 1, window = hanning(window); end if isempty(window), window = hanning(nfft); end if isempty(P{4}), noverlap = 0; else noverlap = P{4}; end if isstr(P{5}) dflag = P{5}; p = []; else dflag = 'none'; if isempty(P{5}), p = .95; else p = P{5}; end end elseif length(P) == 6 % psd(x,nfft,Fs,window,noverlap,p,dflag) if isempty(P{1}), nfft = min(length(x),256); else nfft=P{1}; end if isempty(P{2}), Fs = 2; else Fs = P{2}; end window = P{3}; if length(window) == 1, window = hanning(window); end if isempty(window), window = hanning(nfft); end if isempty(P{4}), noverlap = 0; else noverlap = P{4}; end if isempty(P{5}), p = .95; else p = P{5}; end if isstr(P{6}) dflag = P{6}; else msg = 'DFLAG parameter must be a string.'; return end end
% NOW do error checking if (nfft<length(window)), msg = 'Requires window''s length to be no greater than the FFT length.'; end if (noverlap >= length(window)), msg = 'Requires NOVERLAP to be strictly less than the window length.'; end if (nfft ~= abs(round(nfft)))|(noverlap ~= abs(round(noverlap))), msg = 'Requires positive integer values for NFFT and NOVERLAP.'; end if ~isempty(p), if (prod(size(p))>1)|(p(1,1)>1)|(p(1,1)<0), msg = 'Requires confidence parameter to be a scalar between 0 and 1.'; end end if min(size(x))~=1 | ~isnumeric(x) | length(size(x))>2 msg = 'Requires vector (either row or column) input.'; end if (nargin>2) & ( (min(size(y))~=1) | ~isnumeric(y) | length(size(y))>2 ) msg = 'Requires vector (either row or column) input.'; end if (nargin>2) & (length(x)~=length(y)) msg = 'Requires X and Y be the same length.'; end
dflag = lower(dflag); if strncmp(dflag,'none',1) dflag = 'none'; elseif strncmp(dflag,'linear',1) dflag = 'linear'; elseif strncmp(dflag,'mean',1) dflag = 'mean'; else msg = 'DFLAG must be ''linear'', ''mean'', or ''none''.'; end 是不是这个？

happy · 发表于 2006-6-18 22:39

内部函数是看不了的，可能是技术保密吧

pheigenbau · 发表于 2006-6-18 23:21

我想应该就是了,只是老师怎么得到这些信息的?

happy · 发表于 2006-6-23 20:40

<DIV class=quote>以下是引用pheigenbau在2006-6-18 23:21:06的发言： 我想应该就是了,只是老师怎么得到这些信息的?</DIV>
这些东西有些可以在网上搜到

账号		自动登录	找回密码
密码			我要加入

[综合讨论] [求助]MATLAB求助:psdchk函数的相关帮助

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（pheigenbau）[求助]MATLAB求助:psdchk函数的相...

回复：（pheigenbau）[求助]MATLAB求助:psdchk函数的相...

回复：（pheigenbau）我想应该就是了,只是老师怎么得到...

浏览过的版块