Matlab dll и Delphi

procedure _additorInitialize; stdcall; external 'additor.dll'; 
 function _mlxAdditor(x,y:real):real; stdcall; external 'additor.dll'; 
 
 procedure TForm1.Button1Click(Sender: TObject); 
 var a,b,c:real; 
 begin 
 _additorInitialize; 
 a:=1; 
 b:=2; 
 c:=_mlxAdditor(a,b); 
 Label1.Caption:=floattostr(c); 
 end;

unit Unit2;
 
interface
 
type
    pDouble = ^Double;
    mxComplexty = (mxREAL, mxCOMPLEX);
    pmxArray = ^mxArray;
    mxArray = record end;
 
procedure mlfMatrixMULTIPLICATION(i:Integer; var pout:pmxArray; A,B:pmxArray);cdecl; external 'matr_dll.dll';

mlfMatrixMULTIPLICATION(1,C,A,B);

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Unit2;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    Button2: TButton;
    Button3: TButton;
    Memo1: TMemo;
    Memo2: TMemo;
    Label1: TLabel;
    Label2: TLabel;
    procedure Button1Click(Sender: TObject);
    procedure Button3Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
    procedure FormCreate(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
begin
   if Unit2.Solver11Initialize() then
   Application.MessageBox('Good!','Initialized',MB_OK+MB_ICONINFORMATION)
   else
   Application.MessageBox('Fail!',' Not initialized',MB_OK+MB_ICONERROR);
end;
 
procedure TForm1.Button3Click(Sender: TObject);
begin
   Close();
end;
 
procedure TForm1.Button2Click(Sender: TObject);
var A,B:array of Double;
    AA,BB,CC:pmxArray;
    i:Integer;
    St:String;
begin
   //Ïèøåì óðàâíåíèÿ â ôàéë Solver.txt
   Memo1.Lines.SaveToFile('Solver.txt');
 
   //Âåêòîð íà÷àëüíûõ ïðèáëèæåíèé è ðåçóëüòàòîâ:
   SetLength(A,Memo2.Lines.Count);
   SetLength(B,Memo2.Lines.Count);
   for i:=0 to Memo2.Lines.Count-1 do
   begin
      A[i]:=StrToFloat(Memo2.Lines.Strings[i]);
   end;
 
   AA:=mxCreateDoubleMatrix(1,Length(A),mxReal);
   Move(A[0],mxGetPr(AA)^,Length(A)*SizeOf(Double));
 
   BB:=nil; CC:=nil;
 
   mlfSolver11(2,BB,CC,AA);
 
   Move(mxGetPr(BB)^,B[0],Length(A)*SizeOf(Double));
   St:='';
   for i:=0 to Length(B)-1 do
   begin
      St:=St+FormatFloat('0.0000',B[i]);
      if i<Length(B)-1 then
      St:=St+#$A;
   end;
 
   ShowMessage(St);
 
end;
 
procedure TForm1.FormCreate(Sender: TObject);
begin
   if FileExists('Solver.txt') then
   Memo1.Lines.LoadFromFile('Solver.txt');
end;
 
end.

unit Unit2;
 
interface
 
type
    pDouble = ^Double;
    mxComplexty = (mxREAL, mxCOMPLEX);
    pmxArray = ^mxArray;
    mxArray = record end;
    pString = ^String;
 
//Libmx.dll
function mxCreateDoubleMatrix(m,n:Integer;mxData:mxComplexty):pmxArray; cdecl;
function mxGetPr(arr_ptr:pmxArray):pDouble; cdecl;
function mxGetPi(arr_ptr:pmxArray):pDouble;cdecl;
function mxTranspose(x:pmxArray):pmxArray; cdecl;
procedure mxDestroyArray(arr_ptr:pmxArray); cdecl;
//===
function mxCreateString(S:String):pString; cdecl;
 
function Solver11Initialize():Boolean; cdecl; external 'Solver11.dll';
 
procedure mlfSolver11(i:Integer; var poutX:pmxArray; var poutFVal:pmxArray;  N:pmxArray); cdecl; external 'Solver11.dll';
 
implementation
 
function mxCreateDoubleMatrix; external 'libmx.dll' name 'mxCreateDoubleMatrix';
function mxGetPr; external 'libmx.dll' name 'mxGetPr';
function mxGetPi; external 'libmx.dll' name 'mxGetPi';
function mxTranspose; external 'libmx.dll' name 'mxTranspose';
procedure mxDestroyArray; external 'libmx.dll' name 'mxDestroyArray';
 
//===
function mxCreateString; external 'libmx.dll' name 'mxCreateString';
 
end.

function a = Plusvak(a1, a2)
a = a1 + a2;

unit main;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    Edit1: TEdit;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
  type
    pDouble = ^Double;
    mxComplexty = (mxREAL, mxCOMPLEX);
    pmxArray = ^mxArray;
    mxArray = record end;
 
procedure mlfPlusvak(i:Integer; var pout:pmxArray; A,B:pmxArray);cdecl; external 'my2.dll';
procedure my2Initialize; cdecl; external 'my2.dll';
// procedure _myTerminate; cdecl; external 'testc.dll';
 
 
{*************** .......................... libmx.dll ****************}
{9} function mxCreateDoubleMatrix(m,n:Integer;mxData:mxComplexty):pmxArray; cdecl; external 'libmx.dll';
{10}procedure mxSetName(arr_ptr:pmxArray; const name:String);cdecl; external 'libmx.dll';
{11}function mxGetPr(arr_ptr:pmxArray):pDouble;cdecl; external 'libmx.dll';
{12}function mxGetPi(arr_ptr:pmxArray):pDouble;cdecl; external 'libmx.dll';
{13}procedure mxDestroyArray(arr_ptr:pmxArray); cdecl; external 'libmx.dll';
{14}procedure mxFree(var ram:THandle); cdecl; external 'libmx.dll';
{15}function mxTranspose(x:pmxArray):pmxArray; cdecl; external 'libmx.dll';
{*************** .......................... mclmcr.dll ****************}
{16}function mclInitializeApplication(A:THandle;B:Integer):Boolean;cdecl; external 'mclmcr.dll';
{17}function mclTerminateApplication:Boolean; cdecl; external 'mclmcr.dll';
 
 
var
  Form1: TForm1;
  Flag_MCR_Initialize : Boolean = False;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
const col = 1; row = 1;
var A,B,C: pmxArray;
    AMass,BMass,CMass: array[1..row,1..col] of Double;
    i,j: integer;
begin
  if not Flag_MCR_Initialize then
  begin
   mclInitializeApplication(0,0);
   my2Initialize;
   Flag_MCR_Initialize:=true;
  end;
  randomize;
  for i:=1 to row do
  for j:=1 to col do
  begin
   AMass[i,j]:=random(10)+1;
   BMass[i,j]:=random(10)+1;
  end;
  A := mxCreateDoubleMatrix(row, col, mxREAL);
  Move(AMass, mxGetPr(A)^, row*col*sizeof(double));
  A:=mxTranspose(A);
  B := mxCreateDoubleMatrix(row, col, mxREAL);
  Move(BMass, mxGetPr(B)^, row*col*sizeof(double));
  B:=mxTranspose(B);
  C:=nil;
  i:=StrToInt(Edit1.Text);
  mlfPlusvak(i,C,A,B);
  if C<>nil then
  begin
   C:=mxTranspose(C);
   Move(mxGetPr(C)^, CMass, row*col*sizeof(double));
   ShowMessage('Успешно!');
  end;
end;
 
end.

x=1:100;
y=sin(x);
x=x./y;
y=sin(x);
tic
[a,b]=dst1(y,1);
toc

function [out,nimf] = dst1(y,nn)
if nn==1 % HHT
    options.maxmodes = inf;
    out=emd(y,options);
    nimf=size(out,1);
end
if nn==2 % FFT
    L=max(size(y));
    NFFT = 2^nextpow2(L); % Next power of 2 from length of y
    out = fft(y,NFFT)/L;
    nimf=L;
end
end

  {*************** .......................... mclmcr.dll ****************}
    function mclInitializeApplication(A:THandle;B:Integer):Boolean; cdecl; external 'mclmcr.dll';
    function mclTerminateApplication:Boolean; cdecl; external 'mclmcr.dll';
{******** ........ dll-.................... Untitled1.dll .................... Matlab *******}
    function _Untitled1Initialize:Boolean; cdecl; external 'Untitled1.dll';
    function _Untitled1Terminate :Boolean; cdecl; external 'Untitled1.dll';
    procedure _mlfEdges(i:Integer; var d1:integer; d2:integer); cdecl; external 'Untitled1.dll';

procedure TMainForm.ComputeButtonClick(Sender: TObject);
var
  d1,d2:integer;
begin
  d1:=1;
  d2:=1;
  _mlfEdges(1,d1,d2);
end;