1 files changed, 84 insertions, 0 deletions
diff --git a/frd.sci b/frd.sci
new file mode 100644
index 0000000..c3838b4
--- /dev/null
+++ b/frd.sci
@@ -0,0 +1,84 @@
+function varargout = frd(varargin)
+    [lhs,rhs] = argn(0)
+    if rhs < 2 || rhs > 4 then
+        errmsg = msprintf(gettext("%s: Wrong numbers of input arguments."), "frd");
+        error(errmsg)
+    end
+    frequency = varargin(2)
+    freqUnit = 'rad/TimeUnit'
+    if ~iscolumn(frequency) then
+        errmsg = msprintf(gettext("%s: frequency must be a finite row vector."), "frd");
+        error(errmsg)
+    end
+    respData = varargin(1)
+//    pause
+    if size(frequency,'r') <> size(respData,'r') then
+        errmsg = msprintf(gettext("%s: input output matrix dimension must be equal."), "frd");
+        error(errmsg)
+    end
+    if rhs == 2 then
+        Ts = 0
+    elseif rhs >2 then
+        Ts = varargin(3)
+    end
+    if Ts < 0 || size(Ts,'*') <> 1 || typeof(Ts) <> 'constant' then
+        errmsg = msprintf(gettext("%s: Sampling time must be a scalar non negative real number."), "frd");
+        error(errmsg)
+    end
+    // saving the spectrum value
+    if rhs == 4 then
+        spect = varargin(4)
+    else
+        spect = []
+    end
+    /// matching its dimensions
+    if ~size(spect) then
+    elseif size(frequency,'r') <> size(spect,'r') then
+        errmsg = msprintf(gettext("%s: Numbers of power spectra must be equal to the numbers of frequency."), "frd");
+        error(errmsg)
+    end
+    TUnit = 'seconds'
+    t = tlist(['frd','Frequency','FrequencyUnit','ResponseData','Ts','TimeUnit','Spect'],frequency,freqUnit,respData,Ts,TUnit,spect)
+    varargout(1) = t
+endfunction
+
+//overloading
+function %frd_p(varargin)
+    myTlist = varargin(1)
+    f = fieldnames(myTlist)
+    freqData = myTlist.Frequency
+    tempRespData= myTlist.ResponseData
+    for jj = 1:size(tempRespData,'c')
+        respData = tempRespData(:,jj)
+        mprintf("\t -------------------------")
+        mprintf("\n")
+        mprintf("\t Frequency \t  Response")
+        mprintf("\n")
+        mprintf("\t -------------------------")
+        mprintf("\n")
+        for ii = 1:length(myTlist.Frequency)
+            temp = ''
+            if real(respData(ii))>=0 then
+                temp = temp + ' '
+            end
+            temp = temp + string(real(respData(ii)))
+    //        temp = string(real(respData(ii)))
+            if imag(respData(ii)) > 0 then
+                temp = temp +"+"
+            end
+            if ~imag(respData(ii)) then
+            else
+                temp = temp + string(imag(respData(ii))) +"i"
+            end
+    //        temp = temp + string(imag(respData(ii))) + " i"
+            mprintf("\n\t %f \t %s",freqData(ii),temp)//real(respData(ii)),imag(respData(ii)))
+        end
+        mprintf("\n\n")
+    end
+    if ~myTlist.Ts then
+        mprintf("\n  Continuous Domain frequency response.")
+    else
+        mprintf("\n  Sampling Time = "+string(myTlist.Ts)+" "+myTlist.TimeUnit)
+        mprintf("\n  Discrete Domain frequency response.")
+    end
+endfunction