zmyjiqidaoshujifenzhisrm.m

样条函数方法

function[Zmy] = ZmyjiqidaoshujifenzhiSRM(BCtype,m,n)
disp('1 represent the BCtype is liang duan jian zhi liang ');
disp('2 represent the BCtype is liang duan gu ding liang');
disp('3 represent the BCtype is zuo duan jian zhi you duan gu ding liang');
disp('4 represent the BCtype is zuo duan   zi you you duan gu ding liang');
disp('please choose the type of the boundary condition');
type = input('type=');
switch type
case 1   %the BCtype is liang duan jian zhi liang
  format long g
    ZM = zeros(5,4);
  ZM =[ 0.5 4.9343 48.705 0.6366;
         0.5 19.739 779.28 0;
         0.5 44.413 3945.10 0.2122;
         0.5 78.955 12468 0;
         0.5 123.370 30440 0.1273];
  disp('ZM=');
  disp(ZM);
case 2   %the BCtype is liang duan gu ding liang
     format long g
  ZM = zeros(5,4);
  ZM =[1.0359 12.775 518.52 0.8445;
       0.9984 45.997 3797.1 0;
       1.000 98.920 14619 0.3637;
       1.000 171.58 39940 0;
       1.000 264.01 89138 0.2314;
       1.000 264.01 89138 0.2314];
   disp('ZM=');
   disp(ZM);
case 3   %the BCtype is zuo duan jian zhi you duan gu ding liang
     format long g
  ZM = zeros(5,4);
  ZM = [0.4996 5.5724 118.80 0.6147;
        0.5010 21.451 1250.4 -0.0586;
        0.5000 47.017 5433 0.2364;
        0.5000 82.462 15892 -0.0310;
        0.5000 127.79 36998 0.1464];
    disp('ZM=');
    disp(ZM);
case 4   %the BCtype is zuo duan zi you you duan gu ding liang
     format long g
  ZM = zeros(5,4);
  ZM = [1.8556 8.6299 22.6299 1.0667;
        0.9639 20.176 467.97 0.4252;
        1.0014 77.763 3808.5 0.2549;
        1.000 152.83 14619 0.1819;
        1.000 205.521 39940 0.1415];
    disp('ZM=');
    disp(ZM);
otherwise
  disp('Your choose is wrong');
end