clear,clc
x=[96 104 108 110 113 114 116 118 123 124 125 126 128 129 129 130 130 130 131 132 134 134 136 136 137 138 138 138 139 139 142 142 142 143 144 145 145 145 145 148 149 149 150 150 150 151 151 151 152 153 154 154 154 154 154 155 155 156 156 156 156 157 157 157 158 158 158 158 158 159 160 160 160 161 161 161 162 162 162 162 163 164 164 165 166 166 166 167 167 167 169 169 170 171 171 172 172 172 172 173 173 173 173 173 174 174 174 175 175 175 176 176 176 176 177 177 177 178 179 179 179 180 180 180 181 181 182 183 183 184 184 184 184 185 185 185 185 186 187 187 188 188 188 188 189 189 190 190 190 190 191 191 191 191 192 192 193 193 194 194 194 195 195 195 196 196 197 197 197 197 198 199 200 200 201 201 204 205 205 206 207 211 212 214 216 220 220 226 229];
[yi,xi] = ksdensity(x);
[normalfit,gof1] = fit(xi.', yi.', 'gauss1')
parmhat = wblfit(x);
wbl3 = fittype(@(a,b,c,x) c*b^(-c)*(x-a).^(c-1).*exp(-((x-a)/b).^c), 'coefficients', {'a','b','c'});
[wbl3fit,gof2] = fit(xi.', yi.', wbl3, 'Startpoint', [0, parmhat])
plot(xi, yi, 'b.', xi, normalfit(xi), 'r-', xi, wbl3fit(xi), 'b-');
legend('Data', 'Normal fit', 'Weibull-3 fit')
这个是用3参数韦伯分布拟合的,应为我想换数据所以想尝试用两参数韦伯分布去替换,但是替换后会出现错误!求大神解答!