查看: 2143|回复: 0

全国大学生数学建模竞赛

[复制链接]

字体大小: 正常放大

浅夏110

542 主题	15 听众	1万积分

TA的每日心情

	开心 2020-11-14 17:15

签到天数: 74 天

[LV.6]常住居民II

群组: 2019美赛冲刺课程

群组: 站长地区赛培训

群组: 2019考研数学桃子老师

群组: 2018教师培训（呼伦贝

群组: 2019考研数学站长系列

电梯直达

1^#

发表于 2018-11-12 09:24 |只看该作者 |倒序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

全国大学生数学建模竞赛

一、任务分析

二、题目一求解

题目分析：

模型建立：

模型求解：

3.1读数据，绘图

3.2探测器之间的距离计算

3.3小圆与发射—接收系统的相对运动

3.4提取小圆圆心位置坐标

三、题目二求解

一、任务分析

历年赛题

任务：CT标定，CT图像重建

知识：图像处理

标定模板为椭圆与球。CT测量所得数据为512*180的二维矩阵。

行数512:接收点个数，列数180:CT旋转次数。所成的二维矩阵并不是真实物体的成像，而是不断旋转的水平投影的叠加。

每个角度都有一个大小为512的投影向量，一共180个角度的叠加。

二、题目一求解

求解：CT系统旋转中心在正方形托盘中的位置、探测器单元之间的距离、CT系统使用的X射线的180个方向

题目分析：

标定 CT 系统的安装误差

模型建立：

（1）探测器之间的距离

模型求解：3.1读数据，绘图

表格一数据读取

9 y( c7 `7 ]5 q
, ?! e8 J' S2 K! H" J' r
A=xlsread('E:\math_model_data\data2017\A.xls');( |0 S6 G6 w. ]

# E5 e0 A' S) ^* d3 h- A; M" ^6 H
6 c6 V! h2 z/ N7 O: O
1 ]8 ^6 ]$ D& Y) i
imshow(A)+ ]( q9 z3 }6 ~- v% n l* v

2 L, S( V0 f( Z+ A- z! E1 O; @! n/ K- ~

表格二数据读取

: S$ M$ T3 c' x& ~7 l0 H8 T$ D) q) O: ]* q% ^
B=xlsread('E:\math_model_data\data2017\A.xls',2);
; ~+ e6 @/ \" {+ m. D; u b5 k& y! L$ X& _1 k
) `$ Y5 N6 a. a

" _, m; D) Z1 D) fX=1:size(B,2);4 E L, \' S5 a D; h J

8 W0 o$ g+ q4 b& N
& {: \! B7 l1 ^7 i
J0 Z |; J& O2 WY=1:size(B,1);2 l! P% O7 h" p- ^& G* ]
) G9 N2 Y. d3 N9 S
* l, s2 n4 m# D B/ g7 k9 X: `/ A# y7 X& a* Z! i
figure(1)
# Z" v Z Z3 w; w$ b* I1 t2 O. i- ^+ F; c; O' V
6 Q6 W0 p4 h' Y: R
& Z, ~3 z& F" }: q" ]
mesh(X,Y,B)
5 K* Q9 O& n& S( ?7 D1 Z' I1 W9 S3 p& f$ N

2 r; |6 b, f+ |

3.2探测器之间的距离计算

9 z& F5 w: V- }2 d$ S
1 _6 s7 J0 d$ Y a6 {B=xlsread('E:\math_model_data\data2017\A.xls',2);
7 h1 ~7 m) v1 V* k' H
) Q/ Y! P; [0 c3 ^2 i# S
& ?$ h$ I( {4 M4 r- k) y8 b* e

" M; {+ y H, x) D6 _- R& T
+ G1 _7 Y4 B$ i. ~
, q( t- W+ c# _5 i
1 y/ v. O/ @& \$ |0 E9 ]/ D* |: u% F* X0 `* {. e5 E5 o" O7 O+ a
figure- v2 T3 i" D, H- {0 J2 v; ]8 t

9 t0 g9 n6 k3 ]. ^) Z
! f5 ?# N" c+ {

, ?% L( f; X* R4 D: {lenb = (1:length(B));" _! J3 J3 _6 z& ]7 ^9 u+ ?

: `0 s1 N5 G9 o4 y& |) {* h
+ s q7 s7 S/ F! ~6 b+ v
5 G7 A, d& A6 C D) \( k4 X
# u( I1 F+ J/ K6 l9 S; v* j. K/ |3 d* @) m" Z- d
! {: y. S5 m# X2 q/ m1 R! Q" e7 U7 N7 W2 ?& @, p
subplot(221). V$ A! _. V9 r2 ?7 D: n. p
1 e2 J" |) @' p( {; z
& p# p0 G( ?. U u# H
, z0 A, b; b* S' X/ b9 uscatter(lenb( B(:,1)~=0 ),B( B(:,1) ~= 0 ,1),20,'filled')7 |4 B# e7 Z6 X, W7 F1 Q% w
8 o T7 Y& X/ T' D; r; K2 l& F, Q
) P6 r- [" }4 D. B: F6 Q- b% t, j3 t2 N/ q) p5 P
title('1')& Q M9 _+ [# \- G

& Y/ Y5 b8 n4 v) _& Y5 T
# B9 x5 i) T3 W& _0 S/ _) I
" D$ C1 N a' \! t8 W
0 U7 B8 _9 w- W4 f: ~: z1 I
* d9 C4 S% ~7 K
, I7 r$ E2 X( x+ `$ V$ H6 q. M4 W6 ?
subplot(222)0 J$ y6 {9 ?2 w1 D; a8 p2 A: Q

7 z/ I! A8 s5 @' W/ L
( o7 N1 j7 u* L* r9 P' H ~
" |6 G/ W1 k5 l$ S! z
scatter(lenb( B(:,6)~=0 ),B( B(:,6) ~= 0 ,6),20,'filled')5 Z2 k2 f" F0 r0 D l" O
* {" K* t3 g# C: M' y2 w
4 n7 j. Y% I2 ~! b
j8 D! v ?0 _; ]; M1 H: f
title('6'): l3 ]8 w+ V: a0 Q

# P; x0 [6 M; u# L1 Q! M
$ D: n3 [3 G$ O3 F# l( c9 K' ]
+ d( L4 D$ k) @# j" Y$ D
3 N/ Y0 L( G% y$ l) r1 [% n
! ?) m9 ]0 N* S0 P H* u6 [9 D
/ V% k' {4 {: A' E% Q; R# l6 k" c7 {/ F# h% }
subplot(223)4 q$ P' c) h% u# h2 ]

( m0 i6 k6 v5 A i, t2 E U
( K7 n! O4 ?1 T# X6 g7 m% ]* Q& u; D8 b' h% M
scatter(lenb( B(:,11)~=0 ),B( B(:,11) ~= 0 ,11),20,'filled')
9 @ B$ X( s0 R; A: N7 Y0 Q$ X
. U! J5 i8 L* B: X; Q3 h: z
$ l6 t$ R9 ?8 h) e R0 H
" a( i! W6 ]) a0 w, w( G$ r
title('11')
/ I! l$ p- T9 C, y3 p% v. c0 P" W
* X4 w" u! @: E7 j8 o$ R+ R9 a- I/ o& M" X* O
' u" N8 Z( h8 P& l3 d. V6 _
# B6 I) L- o, V5 R q. _
8 E3 h8 x; `* p ~8 X/ I7 H, h7 f$ O

6 r2 Y9 e( p# q' r4 {2 {! T; Y. e: lsubplot(224)! U7 C/ D9 D+ ]) ~
+ {3 C9 u/ q$ B) M
! W$ h+ t! S) J/ d/ L8 K
; W q; q0 S8 Qscatter(lenb( B(:,60)~=0 ),B( B(:,60) ~= 0 ,60),20,'filled')
, t p0 C- u% R9 O
2 e" s( s* l7 Y4 W
" F* s+ P, ]& @, Y. g9 K* l: J' R* X& f* C4 p; D
title('60重合')
* C- t7 C% T `; f8 ^) O- t% q0 J- \3 e% M, \4 H& b O4 Y

, l1 c! V4 V/ O0 Y* \- w