shuxuexingtai.m

数学形态学进行ECG 滤波处理算法,完全原代码。不是MATLAB 自带的函数。

clear;
ECGBMP30=[0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0104,0105,0106,0108,0109,0109,0110,0111,0111,0110,0110,0110,0110,0109,0108,0106,0105,0104,0103,0101,0101,0100,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0100,0099,0098,0096,0096,0103,0114,0128,0144,0159,0167,0166,0155,0141,0125,0109,0095,0087,0083,0080,0078,0077,0077,0078,0078,0078,0078,0078,0079,0079,0079,0079,0079,0080,0080,0080,0080,0082,0082,0082,0083,0083,0084,0084,0085,0085,0087,0087,0088,0088,0089,0090,0091,0091,0093,0093,0094,0095,0095,0096,0096,0098,0099,0099,0100,0101,0103,0105,0106,0108,0110,0113,0114,0115,0116,0118,0118,0119,0118,0118,0118,0116,0116,0115,0114,0113,0111,0109,0108,0105,0104,0104,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0103,0103,0103,0101,0101,0103,0103,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0101,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0103,0103,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0104,0105,0106,0108,0109,0110,0110,0111,0111,0111,0111,0110,0110,0109,0108,0106,0105,0104,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0100,0099,0098,0096,0096,0103,0115,0128,0144,0159,0167,0166,0155,0140,0125,0108,0094,0087,0083,0080,0078,0078,0078,0078,0078,0078,0078,0078,0079,0079,0079,0079,0079,0079,0080,0080,0080,0082,0082,0082,0083,0083,0084,0084,0085,0085,0087,0087,0088,0089,0089,0090,0090,0091,0091,0093,0094,0094,0095,0096,0096,0098,0099,0099,0100,0101,0103,0104,0106,0108,0110,0113,0114,0116,0116,0118,0118,0119,0119,0118,0118,0116,0116,0115,0114,0113,0111,0109,0108,0106,0105,0104,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0106,0111,0118,0123,0128,0133,0139,0143,0146,0149,0150,0151,0153,0155,0156,0159,0164,0170,0176,0181,0182,0177,0171,0165,0160,0155,0151,0146,0141,0136,0130,0125,0118,0110,0103,0096,0093,0090,0090,0090,0090,0089,0089,0088,0087,0085,0084,0082,0080,0079,0078,0075,0074,0073,0072,0070,0069,0068,0067,0067,0065,0064,0063,0063,0062,0062,0060,0060,0060,0060,0060,0060,0062,0063,0064,0065,0067,0068,0069,0070,0073,0075,0078,0082,0085,0088,0091,0094,0095,0096,0098,0099,0099,0100,0100,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0101,0101,0101,0101,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0103,0103,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0103,0103,0101,0101,0103,0103,0101,0101,0101,0103,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0103,0101,0101,0103,0103,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0104,0105,0106,0108,0109,0109,0110,0110,0111,0111,0110,0110,0109,0109,0108,0106,0105,0104,0103,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0101,0099,0098,0095,0098,0104,0116,0130,0145,0160,0167,0165,0154,0139,0123,0106,0094,0085,0082,0080,0078,0078,0078,0078,0078,0078,0078,0079,0079,0079,0078,0079,0079,0080,0080,0080,0080,0082,0082,0083,0083,0083,0084,0084,0085,0085,0087,0088,0088,0089,0089,0090,0091,0091,0093,0093,0094,0095,0095,0096,0096,0098,0099,0099,0100,0101,0103,0104,0106,0108,0110,0113,0114,0115,0116,0118,0118,0119,0119,0118,0118,0116,0116,0115,0114,0113,0110,0109,0108,0106,0105,0104,0103,0103,0103,0103,0101,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0101,0101,0103,0103,0103,0103,0101,0101,0103,0103,0101,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0103,0103,0101,0101,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0103,0101,0101,0101,0101,0101,0101,0101,0101,0103,0103,0101,0101,0101,0101,0101,0101,0101,0103,0103,0103,0101,0101,0101,0103,0103,0103,0103
];
%x = ECGBMP30;

x=[0152,0135,0144,0152,0143,0154,0130,0140,0144,0149,0131,0134,0146,0120,0135,0130,0146,0149,0130,0139,0127,0136,0136,0124,0135,0143,0127,0136,0147,0127,0139,0119,0134,0136,0139,0125,0126,0132,0116,0137,0143,0167,0189,0189,0208,0196,0189,0174,0144,0138,0132,0119,0128,0135,0125,0141,0120,0136,0136,0138,0123,0126,0136,0119,0136,0130,0139,0144,0130,0146,0132,0135,0134,0124,0132,0144,0137,0132,0149,0132,0144,0128,0139,0138,0135,0128,0128,0144,0128,0140,0138,0146,0154,0135,0147,0139,0147,0141,0136,0146,0156,0136,0149,0158,0138,0151,0132,0150,0147,0147,0136,0137,0141,0124,0138,0132,0135,0141,0124,0136,0128,0141,0135,0122,0135,0146,0127,0135,0144,0124,0136,0120,0139,0141,0137,0123,0127,0143,0119,0134,0132,0135,0143,0127,0146,0134,0138,0134,0122,0131,0143,0128,0137,0144,0125,0137,0123,0135,0138,0134,0125,0128,0141,0121,0135,0127,0140,0146,0125,0138,0130,0138,0132,0124,0132,0144,0127,0136,0140,0127,0139,0120,0135,0143,0140,0132,0135,0149,0125,0140,0136,0146,0155,0137,0150,0138,0146,0146,0131,0135,0141,0127,0134,0143,0132,0131,0123,0136,0139,0141,0127,0126,0134,0120,0130,0130,0141,0146,0130,0139,0132,0136,0132,0122,0130,0132,0125,0147,0168,0158,0189,0186,0202,0196,0181,0158,0140,0138,0109,0125,0120,0132,0143,0126,0134,0126,0131,0131,0124,0138,0143,0127,0139,0149,0126,0138,0123,0132,0136,0137,0134,0136,0141,0120,0136,0137,0140,0149,0135,0140,0134,0139,0146,0134,0140,0147,0132,0140,0152,0134,0147,0130,0139,0144,0151,0137,0140,0150,0131,0147,0146,0152,0158,0138,0154,0144,0144,0143,0134,0143,0146,0134,0139,0147,0125,0137,0117,0132,0136,0134,0124,0131,0138,0121,0132,0131,0135,0141,0124,0137,0131,0135,0135,0124,0135,0147,0128,0134,0139,0124,0135,0120,0138,0128,0138,0126,0130,0144,0121,0131,0123,0136,0139,0127,0141,0131,0137,0134,0125,0132,0140,0128,0137,0140,0127,0139,0123,0127,0136,0138,0124,0127,0139,0123,0132,0128,0138,0146,0124,0137,0132,0138,0134,0127,0135,0147,0137,0154,0154,0135,0151,0135,0141,0144,0146,0130,0132,0143,0127,0140,0131,0135,0144,0132,0138,0131,0139,0131,0124,0132,0149,0134,0138,0143,0125,0137,0123,0135,0138,0139,0123,0124,0137,0114,0137,0143,0167,0189,0187,0212,0198,0186,0174,0149,0136,0134,0119,0130,0135,0124,0137,0122,0131,0135,0134,0124,0127,0137,0120,0138,0130,0139,0146,0130,0137,0130,0134,0134,0128,0134,0146,0132,0140,0152,0132,0139,0122,0136,0138,0141,0136,0126,0147,0126,0141,0141,0146,0150,0128,0146,0135,0144,0147,0136,0146,0154,0143,0150,0155,0140,0154,0130,0147,0151,0151,0130,0137,0149,0123,0137,0131,0140,0143,0126,0138,0128,0130,0130,0124,0131,0136,0125,0131,0139,0126,0144,0122,0132,0139,0141,0125,0127,0141,0117,0132,0128,0143,0151,0128,0136,0130,0139,0132,0124,0135,0138,0128,0136,0150,0131,0139,0121,0131,0136,0138,0126,0130,0141,0117,0134,0134,0138,0147,0125,0138,0130,0138,0135,0125,0130,0144,0131,0134,0141,0127,0140,0119,0135,0138,0143,0128,0134,0144,0126,0143,0137,0147,0159,0137,0151,0139,0144,0137,0128,0134,0143,0132,0136,0144,0127,0140,0126,0136,0136,0136,0124,0126,0139,0125,0127,0132,0139,0147,0134,0140,0128,0130,0132,0119,0127,0137,0125,0146,0166,0165,0189,0185,0205,0198,0180,0155,0144,0141,0104,0124,0125,0132,0141,0126,0140,0128,0134,0134,0124,0127,0138,0128,0136,0137,0124,0135,0119,0135,0146,0139,0126,0134,0149,0122,0137,0136,0138,0147,0131,0150,0140,0140,0136,0128,0141,0147,0135,0146,0147,0134,0146,0135,0146,0149,0146,0134,0138,0154,0132,0149,0146,0150,0156,0143,0152,0143,0144,0143,0134,0143,0149,0135,0137,0147,0128,0136,0119,0132,0137,0134,0124,0127,0140,0116,0131,0126,0136,0144,0125,0137,0135,0135,0135,0124,0134,0138,0127,0132,0140,0128,0137,0121,0134,0138,0141,0126,0126,0137,0117,0131,0130,0141,0137,0130,0139,0132,0140,0135,0122,0125,0139,0125,0135,0146,0127,0139,0121,0135,0136,0136,0126,0130,0137,0121,0136,0134,0131,0144,0130,0137,0131,0137,0139,0128,0137,0149,0137,0139,0150,0136,0147,0124,0139,0140,0140,0131,0140,0146,0121,0138,0136,0137,0144,0131,0137,0130,0136,0139,0131,0134,0139,0130,0140,0141,0127,0140,0119,0134,0137,0144,0128,0126,0134,0109,0135,0147,0167,0190,0187,0207,0196,0191,0171,0149,0137,0132,0120,0130,0137,0125,0135,0122,0135,0132,0135,0124,0130,0137,0119,0134,0132,0135,0143,0126,0138,0132,0135,0135,0130,0135,0146,0132,0140,0143,0130,0138,0123,0139,0140,0141,0131,0136,0152,0126,0138,0135,0141,0150,0136,0152,0132,0147,0146,0137,0150,0154,0139,0141,0154,0136,0151,0136,0147,0150,0149,0136,0137,0146,0125,0139,0130,0139,0147,0130,0132,0128,0136,0130,0122,0130,0141,0126,0134,0140,0126,0132,0117,0134,0136,0131,0123,0124,0137,0120,0140,0132,0136,0149,0132,0137,0130,0138,0131,0123,0131,0147,0136,0136,0139,0126,0141,0120,0134,0139,0134,0125,0127,0146,0123,0136,0130,0135,0143,0130,0139,0132,0137,0131,0124,0136,0141,0131,0134,0141,0127,0140,0121,0135,0136,0141,0132,0131,0144,0127,0144,0137,0146,0154,0139,0146,0137,0139,0139,0131,0135,0144,0136,0136,0144,0127,0140,0117,0134,0134,0136,0128,0128,0140,0121,0136,0136,0139,0144,0126,0136,0128,0136,0139,0116,0130,0136,0127,0150,0167,0160,0182,0183,0202,0196,0185,0156,0144,0139,0110,0124,0123,0134,0144,0124,0138,0131,0137,0126,0123,0132,0138,0127,0135,0144,0125,0137,0120,0135,0132,0135,0126,0130,0137,0120,0132,0131,0140,0155,0134,0141,0138,0146,0138,0128,0140,0144,0135,0141,0155,0140,0146,0125,0141,0151,0146,0136,0141,0150,0131,0147,0151,0154,0160,0140,0149,0141,0149,0144,0135,0141,0146,0134,0141,0144,0128,0135,0117,0132,0137,0135,0125,0124,0140,0120,0130,0127,0136,0146,0125,0137,0130,0137,0130,0123,0128,0139,0130,0132,0141,0131,0138,0122,0134,0138,0134,0125,0125,0138,0122,0134,0131,0137,0147,0134,0139,0128,0134,0131,0122,0132,0147,0122,0138,0144,0130,0143,0121,0131,0130,0135,0122,0127,0141,0121,0136,0131,0138,0144,0127,0139,0134,0136,0138,0131,0140,0143,0137,0146,0150,0135,0147,0131,0140,0143,0143,0132,0128,0140,0122,0135,0125,0135,0140,0125,0138,0138,0138,0134,0127,0137,0140,0130,0138,0140,0126,0138,0125,0140,0138,0134,0124,0128,0132,0111,0138,0141,0166,0187,0191,0214,0201,0187,0170,0146,0140,0135,0121,0130,0135,0124,0139,0119,0134,0134,0134,0125,0128,0138,0120,0131,0132,0138,0141,0126,0138,0134,0134,0135,0125,0136,0140,0131,0135,0144,0131,0139,0123,0141,0141,0141,0131,0135,0143,0124,0137,0135,0146,0152,0136,0147,0141,0151,0146,0134,0141,0151,0139,0149,0162,0132,0154,0134,0149,0154,0147,0134,0130,0144,0121,0137,0135,0138,0147,0127,0138,0128,0132,0132,0124,0127,0140,0130,0137,0136,0125,0138,0116,0132,0135,0137,0123,0126,0138,0120,0127,0126,0136,0144,0122,0135,0126,0132,0132,0131,0134,0139,0132,0139,0141,0126,0140,0117,0132,0136,0140,0131,0127,0136,0120,0138,0128,0137,0149,0125,0137,0130,0140,0138,0126,0131,0138,0128,0137,0143,0128,0140,0119,0134,0141,0139,0132,0131,0144,0126,0144,0138,0147,0152,0139,0150,0136,0140,0139,0131,0134,0143,0130,0138,0140,0126,0135,0119,0135,0135,0135,0131,0128,0141,0121,0136,0127,0136,0141,0126,0140,0131,0136,0134,0124,0132,0136,0124,0140,0162,0160,0189,0189,0198,0200,0183,0158,0149,0139,0106,0119,0122,0128,0141,0128,0138,0131,0134,0134,0123,0130,0140,0130,0132,0143,0127,0140,0114,0134,0138,0135,0126,0127,0143,0121,0135,0132,0140,0143,0128,0141,0135,0134,0135,0125,0134,0147,0143,0143,0149,0137,0151,0126,0141,0150,0143,0136,0138,0156,0138,0149,0140,0150,0164,0139,0151,0146,0143,0143,0134,0146,0152,0137,0137,0140,0126,0139,0119,0134,0137,0132,0123,0128,0137,0120,0128,0127,0135,0146,0126,0137,0127,0136,0135,0121,0128,0139,0131,0131,0140,0126,0140,0117,0132,0134,0134,0126,0125,0138,0124,0134,0131,0137,0146,0125,0137,0127,0132,0136,0124,0131,0140,0132,0140,0143,0125,0136,0117,0131,0137,0140,0119,0130,0140,0122,0138,0131,0134,0138,0125,0135,0131,0139,0138,0131,0137,0150,0137,0141,0151,0137,0144,0128,0143,0147,0136,0131,0134,0139,0121,0134,0132,0135,0143,0127,0139,0126,0131,0134,0125,0126,0138,0126,0132,0141,0135,0140,0120,0137,0141,0136,0125,0127,0132,0110,0134,0151,0171,0189,0182,0208,0202,0186,0171,0151,0136,0134,0120,0134,0140,0126,0137,0116,0131,0137,0135,0126,0128,0136,0119,0136,0128,0137,0141,0126,0137,0134,0135,0136,0123,0135,0144,0130,0135,0144,0132,0139,0122,0137,0143,0137,0130,0130,0143,0126,0137,0136,0149,0152,0136,0147,0140,0140,0143,0131,0140,0155,0141,0149,0156,0141,0156,0132,0144,0146,0143,0132,0137,0151,0119,0139,0134,0139,0150,0127,0135,0123,0131,0127,0122,0132,0140,0130,0134,0141,0125,0136,0119,0134,0132,0135,0126,0128,0131,0119,0134,0126,0135,0143,0128,0136,0128,0132,0134,0120,0127,0140,0130,0128,0138,0123,0135,0119,0140,0138,0135,0128,0132,0138,0120,0136,0127,0136,0144,0132,0144,0132,0131,0134,0128,0130,0140,0134,0132,0141,0126,0144,0123,0136,0138,0136,0128,0135,0146,0127,0144,0136,0144,0156,0136,0149,0136,0139,0139,0131,0135,0144,0128,0137,0144,0124,0136,0120,0136,0134,0135,0125,0131,0137,0120,0131,0128,0137,0143,0127,0143,0131,0136,0135,0124,0124,0134,0122,0139,0167,0164,0189,0185,0207,0202,0182,0154,0140,0136,0106,0124,0127,0125,0144,0126,0138,0135,0134,0130,0117,0127,0136,0128,0136,0141,0128,0138,0121,0134,0136,0136,0128,0126,0141,0124,0138,0126,0138,0150,0128,0140,0134,0140,0137,0127,0135,0149,0131,0137,0149,0135,0139,0124,0138,0143,0144,0143,0140,0151,0135,0152,0141,0150,0162,0137,0150,0141,0151,0151,0135,0136,0147,0138,0136,0143,0131,0135,0119,0132,0141,0138,0126,0126,0135,0117,0134,0128,0136,0146,0124,0136,0132,0132,0135,0121,0128,0139,0131,0132,0141,0124,0139,0121,0130,0135,0134,0127,0124,0137,0119,0136,0127,0135,0140,0125,0138,0128,0132,0138,0124,0132,0141,0132,0131,0140,0124,0136,0122,0134,0137,0135,0127,0134,0139,0119,0131,0127,0134,0144,0132,0131,0135,0137,0139,0135,0137,0146,0131,0140,0147,0135,0150,0128,0143,0144,0141,0131,0131,0141,0123,0131,0130,0137,0147,0122,0137,0132,0131,0134,0123,0132,0139,0128,0134,0143,0122,0135,0120,0134,0132,0132,0122,0122,0132,0120,0137,0144,0168,0192,0185,0208,0202,0185,0170,0146,0143,0140,0121,0123,0136,0128,0135,0119,0135,0132,0134,0124,0132,0141,0121,0134,0126,0135,0146,0126,0139,0134,0132,0134,0127,0132,0144,0130,0135,0144,0132,0140,0124,0135,0143,0140,0127,0130,0143,0127,0137,0135,0141,0154,0132,0144,0137,0140,0144,0132,0140,0158,0143,0149,0156,0140,0149,0131,0141,0147,0147,0135,0136,0147,0127,0143,0132,0136,0141,0124,0134,0130,0136,0125,0125,0130,0141,0134,0132,0138,0120,0135,0114,0131,0138,0135,0126,0126,0138,0119,0131,0128,0137,0140,0126,0139,0132,0126,0132,0125,0127,0139,0128,0135,0139,0125,0137,0120,0127,0134,0134,0126,0122,0135,0116,0131,0128,0143,0147,0126,0140,0136,0132,0134,0127,0127,0140,0130,0138,0149,0127,0135,0120,0138,0137,0137,0134,0130,0143,0126,0149,0141,0147,0154,0132,0146,0139,0140,0140,0131,0131,0141,0134,0134,0144,0124,0136,0120,0136,0136,0136,0124,0130,0139,0117,0131,0128,0138,0141,0126,0138,0135,0131,0132,0120,0124,0136,0123,0139,0170,0164,0189,0185,0206,0195,0180,0154,0140,0140,0109,0124,0123,0134,0149,0126,0136,0127,0128,0130,0123,0134,0134,0131,0134,0144,0132,0138,0117,0126,0135,0131,0125,0131,0140,0123,0136,0132,0138,0147,0130,0143,0132,0138,0139,0131,0128,0146,0137,0136,0147,0132,0147,0125,0139,0144,0146,0131,0135,0150,0132,0140,0139,0146,0154,0138,0158,0146,0144,0147,0138,0137,0147,0138,0134,0143,0126,0141,0124,0132,0131,0132,0128,0124,0137,0122,0128,0126,0134,0150,0128,0138,0130,0128,0132,0124,0128,0141,0131,0130,0139,0127,0136,0120,0130,0134,0134,0126,0126,0138,0117,0135,0130,0132,0141,0125,0139,0127,0131,0134,0126,0127,0139,0127,0131,0143,0124,0137,0124,0134,0137,0136,0127,0124,0137,0117,0132,0131,0136,0144,0127,0140,0137,0136,0136,0127,0132,0146,0138,0146,0144,0137,0149,0130,0147,0144,0139,0125,0128,0137,0121,0136,0130,0137,0144,0127,0137,0130,0132,0136,0122,0130,0143,0134,0126,0140,0128,0136,0119,0132,0138,0134,0124,0123,0134,0106,0131,0143,0166,0182,0181,0205,0196,0186,0180,0150,0137,0137,0125,0124,0137,0127,0134,0117,0131,0139,0139,0125,0122,0137,0123,0131,0128,0138,0141,0125,0137,0138,0137,0137,0125,0127,0141,0134,0135,0146,0132,0138,0122,0139,0139,0140,0127,0130,0143,0127,0138,0136,0139,0154,0135,0143,0137,0140,0146,0132,0139,0154,0146,0143,0154,0137,0149,0132,0143,0147,0151,0136,0137,0146,0126,0135,0130,0134,0143,0127,0137,0130,0131,0134,0127,0128,0138,0126,0128,0138,0125,0140,0110,0134,0136,0135,0130,0126,0135,0111,0130,0124,0134,0146,0126,0138,0130,0132,0132,0123,0127,0141,0127,0132,0143,0128,0131,0119,0134,0134,0132,0124,0127,0136,0119,0132,0130,0130,0140,0126,0138,0126,0130,0130,0122,0128,0149,0134,0132,0144,0132,0137,0120,0137,0135,0136,0130,0137,0151,0127,0138,0137,0149,0152,0135,0151,0136,0138,0139,0136,0137,0146,0131,0130,0140,0127,0137,0121,0135,0134,0134,0127,0126,0139,0117,0131,0128,0138,0144,0127,0136,0134,0135,0131,0121,0124,0137,0123,0138,0164,0166,0186,0182,0202,0196,0182,0154,0140,0143,0109,0125,0123,0132,0139,0124,0134,0127,0132,0132,0124,0128,0141,0136,0132,0140,0124,0134,0117,0134,0140,0127,0128,0128,0141,0126,0136,0130,0131,0144,0125,0139,0137,0137,0139,0128,0135,0147,0135,0137,0151,0131,0144,0128,0143,0138,0141,0137,0135,0149,0131,0147,0140,0147,0155,0139,0144,0139,0143,0144,0127,0134,0144,0132,0135,0151,0128,0137,0121,0137,0134,0132,0127,0123,0135,0119,0136,0134,0135,0138,0124,0140,0128,0131,0136,0121,0126,0139,0136,0135,0143,0125,0134,0117,0132,0135,0135,0127,0123,0136,0122,0132,0130,0132,0140,0125,0139,0130,0132,0131,0126,0130,0138,0128,0131,0144,0124,0136,0119,0136,0134,0134,0123,0124,0139,0117,0132,0135,0136,0144,0127,0139,0130,0134,0135,0127,0137,0149,0138,0140,0154,0140,0147,0126,0139,0139,0137,0131,0135,0134,0123,0135,0131,0139,0144,0124,0131,0128,0127,0132,0126,0130,0143,0131,0132,0141,0126,0136,0121,0130,0136,0136,0126,0117,0131,0114,0130,0141,0162,0187,0181,0206,0197,0186,0167,0143,0136,0136,0115,0121,0132,0121,0135,0125,0134,0135,0135,0130,0124,0136,0122,0130,0127,0134,0149,0132,0138,0127,0132,0138,0123,0128,0147,0130,0132,0144,0136,0144,0125,0136,0136,0136,0131,0130,0143,0127,0136,0134,0143,0150,0135,0143,0137,0139,0146,0134,0140,0151,0146,0146,0154,0137,0149,0134,0141,0146,0143,0138,0132,0143,0123,0135,0132,0134,0141,0130,0137,0131,0131,0134,0120,0126,0136,0127,0131,0140,0125,0136,0120,0137,0135,0132,0123,0122,0134,0120,0136,0121,0136,0143,0127,0140,0130,0128,0126,0121,0123,0139,0132,0131,0143,0126,0136,0119,0131,0134,0135,0123,0126,0138,0122,0126,0126,0136,0140,0124,0136,0132,0128,0132,0123,0130,0136,0127,0131,0143,0121,0132,0116,0130,0136,0143,0132,0131,0147,0132,0139,0137,0147,0150,0134,0146,0143,0144,0140,0125,0131,0149,0130,0131,0146,0124,0135,0119,0138,0139,0136,0125,0123,0137,0122,0132,0130,0137,0140,0125,0139,0130,0134,0131,0120,0123,0137,0125,0138,0160,0164,0187,0181,0202,0196,0182,0154,0139,0138,0114,0121,0122,0126,0139,0125,0134,0128,0135,0134,0124,0127,0141,0127,0128,0138,0124,0138,0120,0132,0136,0135,0131,0127,0138,0120,0131,0128,0138,0151,0122,0143,0135,0136,0143,0128,0131,0139,0134,0132,0147,0135,0144,0128,0140,0144,0141,0135,0135,0152,0128,0146,0143,0151,0151,0137,0151,0139,0140,0143,0135,0135,0146,0134,0135,0138,0123,0135,0120,0125,0130,0128,0122,0123,0143,0121,0131,0130,0138,0140,0124,0139,0127,0128,0132,0126,0132,0140,0126,0128,0144,0124,0135,0122,0128,0132,0132,0131,0128,0138,0120,0127,0127,0135,0141,0126,0139,0127,0130,0136,0123,0130,0138,0128,0130,0144,0125,0137,0117,0135,0136,0131,0123,0125,0139,0117,0131,0128,0138,0140,0125,0135,0130,0135,0135,0127,0139,0149,0139,0139,0154,0132,0144,0125,0138,0143,0139,0131,0130,0141,0126,0134,0128,0132,0139,0123,0137,0135,0124,0136,0124,0130,0146,0131,0128,0136,0124,0131,0119,0135,0136,0135,0125,0123,0131,0111,0131,0144,0159,0186,0185,0211,0191,0183,0174,0144,0135,0134,0122,0121,0135,0123,0136,0114,0127,0132,0134,0120,0121,0134,0116,0130,0136,0136,0143,0127,0141,0130,0131,0138,0122,0127,0143,0137,0138,0147,0125,0138,0126,0135,0138,0140,0127,0128,0141,0131,0141,0137,0140,0147,0132,0146,0138,0140,0146,0131,0138,0155,0141,0144,0154,0137,0149,0134,0144,0147,0141,0137,0136,0143,0123,0135,0134,0132,0140,0124,0139,0126,0128,0130,0121,0127,0137,0127,0135,0140,0125,0136,0120,0128,0132,0128,0123,0125,0137,0120,0131,0130,0139,0143,0124,0135,0126,0127,0134,0126,0121,0141,0131,0131,0147,0126,0134,0111,0130,0130,0131,0127,0125,0139,0120,0132,0128,0134,0141,0127,0135,0131,0132,0136,0119,0126,0143,0128,0130,0140,0127,0134,0119,0131,0138,0131,0126,0130,0146,0122,0136,0135,0140,0151,0141,0149,0138,0139,0146,0127,0131,0147,0130,0130,0141,0130,0141,0121,0128,0134,0137,0124,0124,0140,0119,0131,0128,0140,0147,0128,0137,0127,0130,0135,0121,0124,0137,0123,0136,0165,0160,0186,0181,0202,0196,0182,0156,0141,0136,0111,0124,0121,0127,0139,0126,0132,0127,0128,0136,0120,0125,0137,0128,0131,0139,0124,0140,0120,0134,0136,0135,0124,0125,0136,0121,0135,0132,0137,0146,0131,0146,0135,0135,0136,0124,0130,0147,0140,0128,0151,0134,0144,0131,0140,0140,0136,0132,0131,0147,0134,0144,0143,0149,0155,0137,0150,0141,0143,0141,0132,0137,0151,0128,0131,0144,0124,0135,0119,0132,0130,0130,0123,0125,0131,0115,0130,0128,0127,0137,0122,0132,0128,0137,0135,0123,0128,0146,0128,0128,0144,0123,0134,0119,0136,0139,0134,0121,0123,0139,0119,0130,0132,0131,0139,0125,0141,0134,0132,0134,0120,0125,0141,0130,0130,0144,0123,0135,0121,0131,0136,0131,0123,0124,0139,0120,0132,0127,0136,0143,0124,0135,0130,0135,0135,0126,0132,0152,0136,0137,0149,0132,0146,0126,0139,0147,0139,0132,0130,0141,0120,0131,0125,0132,0143,0126,0137,0130,0132,0139,0124,0126,0139,0127,0127,0141,0130,0128,0122,0134,0136,0138,0124,0120,0125,0109,0126,0140,0165,0185,0185,0207,0198,0185,0171,0146,0137,0132,0121,0123,0138,0117,0132,0120,0127,0131,0131,0125,0122,0135,0119,0132,0123,0130,0140,0127,0131,0126,0127,0132,0123,0136,0146,0134,0134,0150,0127,0138,0126,0134,0137,0136,0134,0135,0143,0122,0136,0138,0139,0149,0137,0141,0137,0138,0151,0137,0139,0154,0139,0140,0155,0138,0150,0134,0140,0144,0144,0135,0135,0141,0123,0135,0134,0135,0143,0123,0137,0130,0126,0130,0121,0127,0137,0127,0127,0144,0122,0134,0116,0128,0135,0128,0123,0128,0136,0120,0131,0128,0131,0139,0122,0134,0130,0130,0134,0122,0127,0146,0131,0128,0139,0122,0132,0120,0135,0126,0135,0126,0125,0140,0120,0128,0122,0131,0137,0124,0138,0130,0132,0134,0123,0126,0140,0130,0131,0139,0126,0137,0122,0126,0134,0137,0126,0128,0141,0128,0137,0136,0140,0154,0128,0141,0137,0138,0135,0124,0127,0140,0130,0137,0144,0125,0137,0124,0130,0134,0136,0123,0123,0137,0124,0137,0130,0130,0140,0128,0135,0130,0135,0131,0119,0122,0140,0128,0137,0164,0155,0183,0183,0202,0196,0182,0154,0139,0139,0109,0123,0120,0126,0138,0126,0135,0128,0126,0135,0123,0124,0138,0128,0131,0138,0123,0135,0123,0130,0134,0130,0124,0126,0136,0121,0138,0132,0137,0146,0131,0138,0132,0131,0136,0127,0131,0147,0136,0136,0154,0132,0141,0125,0136,0139,0141,0138,0127,0151,0132,0144,0147,0147,0154,0132,0146,0137,0139,0146,0132,0137,0149,0135,0131,0143,0125,0137,0116,0130,0134,0132,0117,0121,0137,0116,0127,0126,0134,0138,0123,0135,0130,0124,0128,0121,0127,0135,0125,0125,0138,0123,0141,0121,0131,0135,0136,0123,0123,0139,0117,0130,0128,0137,0147,0126,0132,0128,0134,0132,0122,0130,0138,0128,0128,0149,0128,0137,0120,0127,0132,0134,0124,0125,0139,0117,0130,0131,0134,0143,0124,0135,0130,0135,0137,0128,0131
];

x0=x;

k=[0,1,5,1,0];

n=length(x);

y1=zeros(1,n);
y2=zeros(1,n);
% Open,Close
y3=zeros(1,n);
y4=zeros(1,n);

y5=zeros(1,n);
y6=zeros(1,n);

yy3=zeros(1,n);
yy4=zeros(1,n);

y=zeros(1,5);
max=0;
min=0;

j=1;

%膨胀
for i=5:n-0    
    max=-10000;
    for kl=1:5
        %y(kl)=x(j+kl-1)+k(kl);
        y(kl)=x(j+kl-1)+k(kl);
        if(y(kl)>max)
            max=y(kl);
        end
    end
    y1(j)=max;
    j=j+1;
    %j=j+5;
end

%腐蚀
j=1;
for i=5:n-0    
    min=100000;
    for kl=1:5
        y(kl)=x(j+kl-1)-k(kl);
        if(y(kl)<min)
            min=y(kl);
        end
    end
    y2(j)=min;
    j=j+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%开运算 ,先腐蚀，再膨胀
x=y2;
j=1;
for i=5:n-0     
    max=-10000;
    for kl=1:5
        y(kl)=x(j+kl-1)+k(kl);
        if(y(kl)>max)
            max=y(kl);
        end
    end
    y3(j)=max;
    j=j+1;
end

%闭运算：先膨胀，再腐蚀
x=y1;
j=1;
for i=5:n-0        
    min=10000;
    for kl=1:5