¿Qué método utiliza el sistema CT para dar la tasa de absorción de una sustancia desconocida específica?
El siguiente es el método para calcular la tasa de absorción de CT en matlab (combinado con el tema del Concurso Nacional de Modelado Matemático para Estudiantes Universitarios de 2017 y la calibración de los parámetros de CT):
El reflejo numérico de cada uno punto en la plantilla cuadrada La intensidad de absorción en este punto se llama tasa de absorción.
El medio bidimensional que se va a inspeccionar absorbe la energía del rayo atenuado y, después del procesamiento de ganancia, se obtienen 180 conjuntos de información recibida.
Cuando la intensidad Cuando un haz estrecho de rayos con I_0 atraviesa un objeto con un coeficiente de absorción de nu, su intensidad satisface la relación de reducción exponencial:
I=I_0*exp( -nu*t), donde t es el rayo que atraviesa el espesor de la capa de material. En situaciones reales, el material en estudio a menudo no está compuesto por un solo componente. de varios componentes diferentes, el nu en todas partes dentro del objeto
también puede ser diferente. En tal material, la intensidad del rayo después de pasar a través de este objeto es
I(L). )=I_0*exp(-int_L u(t)dt), donde u (t) es la tasa de absorción en t
El sistema CT analiza la distribución de u(t) dentro del objeto cambiando. un conjunto de trayectorias de rayos L y registra los cambios en la intensidad de emisión correspondiente I(L).
Aquí podemos suponer que la intensidad del haz incidente es 1, entonces el modelo de atenuación de intensidad es: I=ganancia *exp(-nu*t), la ganancia es el coeficiente de ganancia
Puedes saberlo según la plantilla. El coeficiente de absorción del material de la plantilla aquí es nu=1, o 0
Entonces, el modelo de atenuación de intensidad es: I=ganancia*exp(-t)
Usa solo la dirección de 90 grados Usa los datos en una proyección y una plantilla para realizar un ajuste de ganancia por mínimos cuadrados
y: La longitud media en la dirección vertical donde la tasa de absorción de la plantilla correspondiente a cada dirección de proyección en ct_radon90 no es 0
Analizando la proyección vertical de la plantilla con y sin desplazamiento del centro de rotación ,
Se puede saber que la información recibida rxv después de la transformación del radón es proporcional a la longitud y de los rayos X que pasan a través del medio, y el coeficiente proporcional es Coeficiente de ganancia ganancia
De acuerdo con la tasa de absorción del medio plantilla en el Apéndice 1, se puede calcular la ganancia de amplificación del sistema
p1=find(ct_radon90~=0);
p2=p1 (2:end);
p=find(p2-p1(1:end-1)~=1);
rxv es igual que la elipse en la plantilla Parcialmente valores de absorción de rayos X correspondientes después de la ganancia
rxv=ct_radon90(p1(1:find(p2-p1(1:end-1)~=1))); >
Calcular la longitud y del medio, calcular la longitud y del medio en la elipse -15lt;=xlt;=15
x=-15:30/(length(rxv)- 1):15 ;
Ecuación elíptica: x^2/225 y^2/1600=1
y=2*40*sqrt((1-x.^2/ 225));
Elimina elementos cuyas longitudes en ambos extremos de y sean relativamente cortas o 0, y las operaciones rxv correspondientes sean las mismas
y=y(10:length(y). )-10);
rxv=rxv(10:longitud(rxv)-10);
ganancia=rxv(10)/y(10);
[ganancia, res]= lsqnonlin(@(ganancia)myfun(ganancia, rxv, y'), ganancia);
xlswrite('coeficiente de ganancia del sistema CT.xls', ganancia);