sjbh1.c

其中包括了matlab中最基本的算法院程序

     */
    mlfAssign(
      &h,
      mlfNLine(
        1,
        _mxarray12_,
        _mxarray14_,
        _mxarray16_,
        _mxarray18_,
        _mxarray20_,
        _mxarray22_,
        _mxarray23_,
        _mxarray25_,
        NULL));
    /*
     * w=1:0.05:2;
     */
    mlfAssign(&w, _mxarray27_);
    /*
     * y2=zeros(1101,21);
     */
    mlfAssign(&y2, mlfZeros(_mxarray29_, _mxarray30_, NULL));
    /*
     * for i=1:21
     */
    {
        int v_ = mclForIntStart(1);
        int e_ = 21;
        if (v_ > e_) {
            mlfAssign(&i, _mxarray31_);
        } else {
            /*
             * y2(:,i)=sin(w(i)*x');
             * set(h,'xdata',x,'ydata',y2(:,i));
             * drawnow;
             * pause(0.0005)
             * i=i+1;
             * %axis=theax;
             * end
             */
            for (; ; ) {
                mlfAssign(&i, mlfScalar(v_));
                mclArrayAssign2(
                  &y2,
                  mlfSin(
                    mlf_times_transpose(
                      mclArrayRef1(mclVv(w, "w"), mclVv(i, "i")),
                      mclVv(x, "x"),
                      _mxarray7_)),
                  mlfCreateColonIndex(),
                  mclVv(i, "i"));
                mclAssignAns(
                  &ans,
                  mlfNSet(
                    0,
                    mclVv(h, "h"),
                    _mxarray32_,
                    mclVv(x, "x"),
                    _mxarray34_,
                    mclArrayRef2(
                      mclVv(y2, "y2"), mlfCreateColonIndex(), mclVv(i, "i")),
                    NULL));
                mlfDrawnow(NULL);
                mlfPause(_mxarray36_);
                mlfAssign(&i, mclPlus(mclVv(i, "i"), _mxarray37_));
                if (v_ == e_) {
                    break;
                }
                ++v_;
            }
        }
    }
    /*
     * legend off
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray8_, NULL));
    /*
     * legend('sin(2x)')
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray38_, NULL));
    /*
     * pause(2)
     */
    mlfPause(_mxarray7_);
    /*
     * legend('将y=sin(2x)向左平移pi/6')
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray40_, NULL));
    /*
     * %hold on
     * pause
     */
    mlfPause(NULL);
    /*
     * %y=sin(wx)->y=sin(wx+pi/3)
     * h=line('Color',[1 0 1],'Marker','.','MarkerSize',4,'EraseMode','Xor');
     */
    mlfAssign(
      &h,
      mlfNLine(
        1,
        _mxarray12_,
        _mxarray42_,
        _mxarray16_,
        _mxarray18_,
        _mxarray20_,
        _mxarray22_,
        _mxarray23_,
        _mxarray44_,
        NULL));
    /*
     * w1=2;
     */
    mlfAssign(&w1, _mxarray7_);
    /*
     * y3=zeros(1101,20);
     */
    mlfAssign(&y3, mlfZeros(_mxarray29_, _mxarray46_, NULL));
    /*
     * for i=1:20
     */
    {
        int v_ = mclForIntStart(1);
        int e_ = 20;
        if (v_ > e_) {
            mlfAssign(&i, _mxarray31_);
        } else {
            /*
             * y3(:,i)=sin(w1*x'+i*pi/60);
             * set(h,'xdata',x,'ydata',y3(:,i));
             * drawnow;
             * pause(0.0005)
             * i=i+1;
             * %axis=theax;
             * end
             */
            for (; ; ) {
                mlfAssign(&i, mlfScalar(v_));
                mclArrayAssign2(
                  &y3,
                  mlfSin(
                    mclPlus(
                      mlf_times_transpose(
                        mclVv(w1, "w1"), mclVv(x, "x"), _mxarray7_),
                      mclMrdivide(
                        mclMtimes(mclVv(i, "i"), _mxarray47_), _mxarray48_))),
                  mlfCreateColonIndex(),
                  mclVv(i, "i"));
                mclAssignAns(
                  &ans,
                  mlfNSet(
                    0,
                    mclVv(h, "h"),
                    _mxarray32_,
                    mclVv(x, "x"),
                    _mxarray34_,
                    mclArrayRef2(
                      mclVv(y3, "y3"), mlfCreateColonIndex(), mclVv(i, "i")),
                    NULL));
                mlfDrawnow(NULL);
                mlfPause(_mxarray36_);
                mlfAssign(&i, mclPlus(mclVv(i, "i"), _mxarray37_));
                if (v_ == e_) {
                    break;
                }
                ++v_;
            }
        }
    }
    /*
     * legend off
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray8_, NULL));
    /*
     * legend('sin(2x+pi/3)')
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray49_, NULL));
    /*
     * pause(2)
     */
    mlfPause(_mxarray7_);
    /*
     * legend('将y=sin(2x)向左平移pi/6')
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray40_, NULL));
    /*
     * %y=sin(wx+pi/3)->y=2sin(wx+pi/3)
     * pause
     */
    mlfPause(NULL);
    /*
     * %y=sin(wx)->y=sin(wx+pi/3)
     * h=line('Color',[1 0 1],'Marker','.','MarkerSize',5,'EraseMode','Xor');
     */
    mlfAssign(
      &h,
      mlfNLine(
        1,
        _mxarray12_,
        _mxarray42_,
        _mxarray16_,
        _mxarray18_,
        _mxarray20_,
        _mxarray51_,
        _mxarray23_,
        _mxarray44_,
        NULL));
    /*
     * w3=1:0.05:2;
     */
    mlfAssign(&w3, _mxarray27_);
    /*
     * y4=zeros(1101,21);
     */
    mlfAssign(&y4, mlfZeros(_mxarray29_, _mxarray30_, NULL));
    /*
     * for i=1:21
     */
    {
        int v_ = mclForIntStart(1);
        int e_ = 21;
        if (v_ > e_) {
            mlfAssign(&i, _mxarray31_);
        } else {
            /*
             * y4(:,i)=w3(i)*sin(w1*x'+pi/3);
             * set(h,'xdata',x,'ydata',y4(:,i));
             * drawnow;
             * %pause(0.0005)
             * i=i+1;
             * %axis=theax;
             * end
             */
            for (; ; ) {
                mlfAssign(&i, mlfScalar(v_));
                mclArrayAssign2(
                  &y4,
                  mclMtimes(
                    mclArrayRef1(mclVv(w3, "w3"), mclVv(i, "i")),
                    mlfSin(
                      mclPlus(
                        mlf_times_transpose(
                          mclVv(w1, "w1"), mclVv(x, "x"), _mxarray7_),
                        _mxarray52_))),
                  mlfCreateColonIndex(),
                  mclVv(i, "i"));
                mclAssignAns(
                  &ans,
                  mlfNSet(
                    0,
                    mclVv(h, "h"),
                    _mxarray32_,
                    mclVv(x, "x"),
                    _mxarray34_,
                    mclArrayRef2(
                      mclVv(y4, "y4"), mlfCreateColonIndex(), mclVv(i, "i")),
                    NULL));
                mlfDrawnow(NULL);
                mlfAssign(&i, mclPlus(mclVv(i, "i"), _mxarray37_));
                if (v_ == e_) {
                    break;
                }
                ++v_;
            }
        }
    }
    /*
     * legend off
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray8_, NULL));
    /*
     * legend('2sin(2x+pi/3)')
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray53_, NULL));
    /*
     * legend off
     */
    mclPrintAns(&ans, mlfNLegend(0, NULL, NULL, NULL, _mxarray8_, NULL));
    mxDestroyArray(x);
    mxDestroyArray(y1);
    mxDestroyArray(ans);
    mxDestroyArray(h);
    mxDestroyArray(w);
    mxDestroyArray(y2);
    mxDestroyArray(i);
    mxDestroyArray(w1);
    mxDestroyArray(y3);
    mxDestroyArray(w3);
    mxDestroyArray(y4);
    mxDestroyArray(m);
    mclSetCurrentLocalFunctionTable(save_local_function_table_);
    /*
     * %pause(2)
     * %legend('将y=sin(2x)向左平移pi/6')
     */
}