2017年数学建模A组练习代码

clear,clc
%%表5处理图形
D=xlsread('D:\sparename\Documents\2017\A\A.xls',5)
% imagesc(D)
R=iradon(D,[0:179]+30,512)
imagesc(R)

clear,clc
%%表1处理图形
A=xlsread('D:\sparename\Documents\2017\A\A.xls',1)
imshow(A)
%%附件2
figure
B=xlsread('D:\sparename\Documents\2017\A\A.xls',2)
imagesc(B)

题目1

%% 求探测器距离d
B=xlsread('A题附件.xls',2);
imagesc(B)
W=sum(B>0);
d=80/max(W)
%% 求旋转中心偏移量
imin=min(W),jmin=find(W(:)==imin);%得投影最小长度imin探测器
i=jmin(1);
x=find(B(:,i)>0)%此时对应角度的所有接受的探测器
% 512/2得到近似中心探测器编号，而(max(x)+min(x))/2)得到椭圆中心对应的探测器编号。编号之差乘d得到距离。
nx=(max(x)+169)
idx=(512/2-(max(x)+169)/2)*d;
jmax=150-90;
y=find(B(:,jmax)>0);%此时对应角度的所有接受的探测器
idy=(512/2-(max(y)+min(y))/2)*d;
%%求偏移角度
alpha=180-jmin
j=jmin(2)
x=find(B(:,j)>0);%此时对应角度的所有接受的探测器
jdx=(512/2-(max(x)+169)/2)*d;
jmax=150-90
y=find(B(:,jmax)>0)%此时对应角度的所有接受的探测器
jdy=(512/2-(max(y)+min(y))/2)*d;
dx=-(idx+jdx)/2
dy=(idy+jdy)/2%%求偏移角度
alpha=180-jmin
%% 验证
B=[zeros(200,180);B;zeros(200,180)]theta=[0:1:179]+30
% 基于R-L滤波器的滤波反投影法
R= iradon(B,theta,'Hann')
figure
imshow(R)

题目2

clear,clc
%%表三处理
C=xlsread('D:\sparename\Documents\2017\A\A.xls',3)
C=[zeros(200,180);C;zeros(200,180)];
R=iradon(C,[0:179]+30,'Hann');
imagesc(R)
d=0.2768
dx=-9.2042
dy=5.8131
x=ceil(-dx/d)+size(R,1)/2+[ceil(-50/d):ceil(50/d)];
y=ceil(dy/d)+size(R,1)/2+[ceil(-50/d):ceil(50/d)];
% 调整大小
R=imrescale(R(y,x),[256,256])
T=flipud(R)
D=xlsread('D:\sparename\Documents\2017\A\A.xls',4)
d=100/256
c=ceil(D(:,1)/d)
r=ceil(D(:,2)/d)
ix=sub2ind(size(T),r,c);
rate=[D T(ix)];
imagesc([0,100],[0,100],T);
colorbar;
hold on
plot(D(T(ix)==0,1),D(T(ix)==0,2),'wx');
plot(D(T(ix)~=0,1),D(T(ix)~=0,2),'r+');
set(gca,'ydir','normal')
%
% figure
% imagesc([0,179]+30,[1,512],R)
%函数
function new=imrescale(old,newsize)
scale=numel(old)/prod(newsize);
new=imresize(old,newsize)*scale;
new(new<1e-1)=0;
end

2017年数学建模A组练习代码相关推荐

数学建模——灰色预测模型Python代码
数学建模--灰色预测模型Python代码 """ Spyder Editor This is a temporary script file. ""& ...
数学建模入门方法及代码
数学建模入门数学模型 (1)原型与模型模型是原型的替代物.(对原型提炼构造.) 模型不是原型,既简单于原型又高于原型. 模型的分类数学模型是对特定对象,特定目的,根据特有的内在规律,做出一些必要 ...
2022大湾区杯奥港金融数学建模竞赛思路及代码
占个位置吧,开始在本帖实时更新赛题思路代码,文章末尾获取! 持续为更新参考思路赛题思路会持续进行思路模型分析,下自行获取. A题思路: (比赛开始后第一时间更新) B题思路: (比赛开始后第一时间 ...
17届华为杯数学建模大赛B题代码
1 导入库: import numpy as np import pandas as pd import seaborn as sns import matplotlib.pyplot as plt ...
数学建模常用模型及代码
一.规划模型 1.线性规划线性规划与非线性规划问题一般都是求最大值和最小值,都是利用最小的有限资源来求最大利益等,一般都利用lingo工具进行求解. 点击进入传送门 2.整数规划求解方式类似于线性 ...
19年全国数学建模比赛A题代码（简单的迭代思想）
前言:怀念通宵刷题的日子,困得一批第一题 M = 10000; %循环次数 %需要用到的函数 %function E = EMPA(p) 计算弹性模量 p为当前时刻压强 %funtion QA = ...
2018国赛数学建模B题一道工序代码
问题B 智能RGV的动态调度策略图1是一个智能加工系统的示意图,由8台计算机数控机床(Computer Number Control ...
全国数学建模C题仿真代码
CUMCM已经结束了,我们最后的查重率是0%,但是我们搞错了论文,提交了一个18.3%的论文上去.给上我们队的排队论仿真代码,顺便给大家当一个警示!不要踩我们的坑! s=126; onetimefuw ...
2019数学建模国赛Apython代码
第二问: 凸轮角速度0.042,和优秀论文有差距,emmm,调了半天,大概率是前面拟合的问题. 想要图好看一点,可以自己悄悄改一下pb 代码有注释写了很多没用的列表,可以去掉 import math ...

2017年数学建模A组练习代码

2017年数学建模A组练习代码相关推荐

最新文章

热门文章