z=4;
M=zeros(z,64);
for i= 1:1:64
filename = ['D:\iCentral\pictures\guti\3_blockId#' num2str(i), '.bin'];
fid = fopen(filename,'r');
[Vector, count] = fread(fid,'uint8');
I = reshape(Vector,64,64)';
I=I(20:21,33:34);
I=double(I)/255;
Vector=reshape(I,1,z)';
M(:,i) = Vector;
end
sum=[];
for j=1:2:63
sum=[sum M(:,j)+M(:,j+1)];
end
A=sum;%2ms像素矩陣，每一列為一幀影像
sum=[];
for j=1:2:31
sum=[sum A(:,j)+A(:,j+1)];
end
B=sum;%4ms像素矩陣
sum=[];
for j=1:2:15
sum=[sum B(:,j)+B(:,j+1)];
end
C=sum;%8ms像素矩陣
sum=[];
for j=1:2:7
sum=[sum C(:,j)+C(:,j+1)];
end
D=sum;%16ms像素矩陣
sum=[];
for j=1:2:3
sum=[sum D(:,j)+D(:,j+1)];
end
E=sum;%32ms像素矩陣
F=E(:,1)+E(:,2);%64ms像素矩陣
K0=mean(sqrt((std(M,0,2))./(mean(M,2))));
K1=mean(sqrt((std(A,0,2))./(mean(A,2))));
K2=mean(sqrt((std(B,0,2))./(mean(B,2))));
K3=mean(sqrt((std(C,0,2))./(mean(C,2))));
K4=mean(sqrt((std(D,0,2))./(mean(D,2))));
K5=mean(sqrt((std(E,0,2))./(mean(E,2))));
K6=mean(sqrt((std(E,0,2))./(mean(F,2))));
syms b p t v
x0=1/t;x1=2/t;x2=4/t;x3=8/t;x4=16/t;x5=32/t;x6=64/t;
s0=b*p^2*(exp(-2*x0)-1+2*x0)/(2*x0^2)+4*b*p*(1-p)*(exp(-x0)-1+x0)/x0^2+v-K0^2;
s1=b*p^2*(exp(-2*x1)-1+2*x1)/(2*x1^2)+4*b*p*(1-p)*(exp(-x1)-1+x1)/x1^2+v-K1^2;
s2=b*p^2*(exp(-2*x2)-1+2*x2)/(2*x2^2)+4*b*p*(1-p)*(exp(-x2)-1+x2)/x2^2+v-K2^2;
s3=b*p^2*(exp(-2*x3)-1+2*x3)/(2*x3^2)+4*b*p*(1-p)*(exp(-x3)-1+x3)/x3^2+v-K3^2;
s4=b*p^2*(exp(-2*x4)-1+2*x4)/(2*x4^2)+4*b*p*(1-p)*(exp(-x4)-1+x4)/x4^2+v-K4^2;
s5=b*p^2*(exp(-2*x5)-1+2*x5)/(2*x5^2)+4*b*p*(1-p)*(exp(-x5)-1+x5)/x5^2+v-K5^2;
s6=b*p^2*(exp(-2*x6)-1+2*x6)/(2*x6^2)+4*b*p*(1-p)*(exp(-x6)-1+x6)/x6^2+v-K6^2;
[b,p,t,v] = solve(s0,s1,s2,s3,s4,s5,s6,b,p,t,v)
警告: Cannot find explicit solution.
> In solve (line 316)
In hecheng (line 65)
b =
Empty sym: 0-by-1
p =
Empty sym: 0-by-1
t =
Empty sym: 0-by-1
v =
Empty sym: 0-by-1
matlab將.bin格式的影像檔案合成，計算對比值，在根據公式算其他的未知引數，結果一直無解，請問那個大佬能幫幫忙
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/284440.html
標籤：其他開發語言
上一篇：cmd命令執行"python test_login.py"，系統報錯：E Time Elapsed 0:00:00