What should I do to correct user defined function?

Question

Explorer 2016-3-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/273495-what-should-i-do-to-correct-user-defined-function

编辑： Explorer 2016-3-15

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

I have written code for function file i.e nsr_f.m

Getting following error when calling it:

Error using arburg (line 72)
The number of rows of X must at least 5.
Error in nsr_f (line 33)
[d1,p1,rc] = arburg(seg,4);
Error in feature_matching (line 3)
nsr_f(fs,Ekg,Gain)

Code to call function:

Gain = 200; fs=128;
d = load('PhysioNet_Database\NSR\16272m.mat'); Ekg = d.val(1,:);
nsr_f(fs,Ekg,Gain)

2 个评论
显示无隐藏无

Explorer 2016-3-14

编辑：Explorer 2016-3-14

16272m.mat

MAT file also attached

Image Analyst 2016-3-14

I'm pretty sure this is the route to solving it: Click me

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2016-3-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/273495-what-should-i-do-to-correct-user-defined-function#answer_213602

You do not pre-initialize seg, so it is growing as you execute your loop. You start with i = 1 and assign something to seg(i,:) . That would assign to seg(1,:) the first time, so your seg is going to start out as a single row. You then pass that single row to arburg(), but arburg() requires at least 5 rows to work on, so it fails.

If you somehow succeeded on that first call, then you would go to the next iteration of the loop and would assign to seg(i,:) which would be seg(2,:) so you would be growing seg to two rows, and passing that to arburg(), which would then be processing not only this new row but the previous row. You would thus be incrementally refining your coefficients as you want through the entire EEG, first output based upon the first segment, second output based upon the first two segments, third output based upon the first three segments, and so on. Are you sure that is what you want?

I suspect you want to populate the seg matrix completely before calling arburg()

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Explorer 2016-3-15

编辑：Explorer 2016-3-15

在 MATLAB Online 中打开

I found the solution of my problem because of what your wrote above. Thank you!

I replaced following things:

 % function [seg]=nsr_f(fs,Ekg,Gain) with 
function [rc_seg]=nsr_f(fs,Ekg,Gain) 
 % and [d1,p1,rc] = arburg(seg,4) with
 [d1,p1,rc] = arburg(seg(i,:),4);

请先登录，再进行评论。