Guest User

Untitled

a guest
Jun 19th, 2019
115
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. function deconvfail
  2.  
  3. y = [0.068 0.050 0.049 0.087 0.088 0.056 0.063 0.049 0.010 0.039 0.093 0.105 0.103 0.130 0.151 0.126 0.097 0.098 0.095 0.066 0.061 0.072 0.053 0.047 0.042 0.033 0.053 0.054 0.041 0.055 0.070 0.078 0.106 0.127 0.097 0.077 0.103 0.086 0.050 0.075 0.100 0.092 0.123 0.158 0.166 0.159 0.130 0.132 0.125 0.120 0.152 0.139 0.113 0.128 0.150 0.141 0.109 0.103 0.094 0.044 0.027 0.084 0.092 0.073 0.103 0.101 0.098 0.116 0.111 0.118 0.105 0.070 0.057 0.032 0.011 0.037 0.063 0.091 0.177 0.253 0.296 0.349 0.419 0.484 0.530 0.610 0.710 0.786 0.874 0.977 1.083 1.198 1.316 1.408 1.498 1.593 1.644 1.717 1.828 1.902 1.983 2.115 2.213 2.284 2.401 2.506 2.596 2.703 2.814 2.919 3.008 3.100 3.201 3.316 3.454 3.587 3.690 3.779 3.883 3.985 4.065 4.181 4.325 4.410 4.504 4.659 4.801 4.950 5.092 5.177 5.283 5.358 5.390 5.514 5.657 5.755 5.895 6.039 6.114 6.197 6.317 6.419 6.512 6.595 6.646 6.647 6.612 6.549 6.437 6.287 6.085 5.859 5.630 5.393 5.190 4.986 4.776 4.594 4.402 4.175 3.957 3.740 3.495 3.326 3.212 3.047 2.914 2.838 2.724 2.570 2.455 2.354 2.233 2.140 2.082 2.043 1.964 1.861 1.773 1.662 1.570 1.527 1.480 1.412 1.363 1.316 1.253 1.219 1.176 1.143 1.166 1.145 1.102 1.124 1.143 1.114 1.109 1.139 1.156 1.122 1.071 1.075 1.067 1.028 1.019 0.982 0.889 0.774 0.657 0.580 0.548 0.525 0.531 0.531 0.496 0.500 0.514 0.487 0.495 0.518 0.493 0.487 0.503 0.490 0.468 0.481 0.492 0.476 0.509 0.583 0.629 0.678 0.751 0.786 0.791 0.793 0.791 0.792 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791 0.791];
  4. u = [0.087 0.076 0.076 0.077 0.080 0.084 0.083 0.087 0.091 0.085 0.088 0.095 0.088 0.086 0.089 0.083 0.081 0.087 0.086 0.082 0.092 0.100 0.097 0.101 0.110 0.108 0.097 0.090 0.086 0.081 0.074 0.070 0.072 0.077 0.075 0.077 0.090 0.093 0.087 0.088 0.089 0.079 0.079 0.086 0.088 0.092 0.087 0.082 0.074 0.068 0.077 0.085 0.101 0.105 0.113 0.113 0.103 0.104 0.085 0.104 0.081 0.080 0.094 0.058 0.124 -0.026 0.432 2.301 4.290 6.031 8.143 9.537 9.564 9.504 9.568 9.508 9.536 9.516 9.512 9.506 9.495 9.510 9.495 9.497 9.496 9.497 9.494 9.486 9.494 9.493 9.489 9.484 9.476 9.469 9.467 9.475 9.476 9.481 9.489 9.484 9.474 9.476 9.483 9.482 9.479 9.469 9.459 9.458 9.452 9.454 9.465 9.465 9.453 9.450 9.452 9.445 9.446 9.450 9.448 9.450 9.448 9.444 9.447 9.456 9.454 9.458 9.472 9.468 9.465 9.453 9.454 9.449 9.431 9.447 9.440 9.474 9.447 9.485 9.450 8.281 6.267 4.445 2.515 0.638 -0.068 0.009 -0.010 0.007 0.028 0.033 0.053 0.040 0.047 0.045 0.051 0.056 0.057 0.063 0.052 0.048 0.042 0.029 0.026 0.031 0.034 0.039 0.048 0.051 0.052 0.052 0.053 0.054 0.052 0.050 0.047 0.047 0.053 0.057 0.059 0.063 0.061 0.059 0.058 0.056 0.053 0.049 0.050 0.052 0.056 0.062 0.059 0.054 0.051 0.050 0.055 0.063 0.066 0.075 0.074 0.057 0.052 0.056 0.049 0.044 0.053 0.060 0.061 0.068 0.075 0.074 0.070 0.063 0.063 0.066 0.059 0.053 0.053 0.048 0.040 0.043 0.050 0.045 0.040 0.049 0.053 0.050 0.055 0.053 0.042 0.037 0.032 0.027 0.028 0.029 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028 0.028];
  5.  
  6. ir = real(ifft(fft(y)./fft(u)));
  7. Y = conv(ir, u);
  8. Y = Y(1:length(y));
  9.  
  10. clf;
  11. subplot(2,1,1);
  12. hold on
  13. plot(y, 'LineWidth', 2);
  14. plot(Y, 'LineWidth', 2);
  15. grid on
  16. legend('system output', 'reconstructed output');
  17. subplot(2,1,2)
  18. hold on
  19. plot(u, 'LineWidth', 2);
  20. plot(ir, 'LineWidth', 2);
  21. grid on
  22. legend('system input', 'impulse response');
RAW Paste Data