그리고 그래프로 원래값과 비교해 본다.
rom_iq = fopen("rom_iq.txt",'r');
rom_d = fopen("rom_d.txt",'r');
a = fscanf(rom_iq,'%f');
b = fscanf(rom_d,'%f');
k = 0;
for i = 1:8192
iq(i) = i^(4/3);
nf = i/128;
n = fix(nf);
if(nf == n)
n = n-1;
end
if(n<1)
k(i)=a(i);
else
if(n<4)
k(i)=b(n)+a(i-128*n)+(i-128*n)*b(n+64)/5.5;
else
k(i)=b(n)+a(i-128*n)+(i-128*n)*b(n+64);
end
end
end
subplot(2,1,1);
plot(iq);
subplot(2,1,2);
plot(k);
댓글 달기