clc;
clear;
t0 = clock;
%求解参数定义
format long;
Bz=5e-4;
Br=5e-4;
J1=155e6;
J2=155e6;
B0=1.5;
e=5e-6;
%meshgrid of the feasible coil space
%线圈区域1
z1_min=0;
z1_max=1.32/2;
r1_min=0.5;
r1_max=1;
n1_cell_r=100;
n1_cell_z=100;
s1_cell_r=(r1_max-r1_min)/n1_cell_r;
s1_cell_z=(z1_max-z1_min)/n1_cell_z;
area1_cell=s1_cell_r*s1_cell_z;
%
[R1,Z1]=meshgrid(r1_min+s1_cell_r/2:s1_cell_r:r1_max-s1_cell_r/2,z1_min+s1_cell_z/2:s1_cell_z:z1_max-s1_cell_z/2);
%线圈区域2
z2_min=0;
z2_max=1.32/2;
r2_min=0.5;
r2_max=1;
n2_cell_r=100;
n2_cell_z=100;
s2_cell_r=(r2_max-r2_min)/n2_cell_r;
s2_cell_z=(z2_max-z2_min)/n2_cell_z;
area2_cell=s2_cell_r*s2_cell_z;
[R2,Z2]=meshgrid(r2_min+s2_cell_r/2:s2_cell_r:r2_max-s2_cell_r/2,z2_min+s2_cell_z/2:s2_cell_z:z2_max-s2_cell_z/2);
%对称线圈
[R5,Z5]=meshgrid(r1_min+s1_cell_r/2:s1_cell_r:r1_max-s1_cell_r/2,-z1_min-s1_cell_z/2:-s1_cell_z:-z1_max+s1_cell_z/2);
[R6,Z6]=meshgrid(r2_min+s2_cell_r/2:s2_cell_r:r2_max-s2_cell_r/2,-z2_min-s2_cell_z/2:-s2_cell_z:-z2_max+s2_cell_z/2);
%将矩阵转换为向量
R1_v=reshape(R1,1,n1_cell_z*n1_cell_r);
Z1_v=reshape(Z1,1,n1_cell_z*n1_cell_r);
R2_v=reshape(R2,1,n2_cell_z*n2_cell_r);
Z2_v=reshape(Z2,1,n2_cell_z*n2_cell_r);
R5_v=reshape(R5,1,n1_cell_z*n1_cell_r);
Z5_v=reshape(Z5,1,n1_cell_z*n1_cell_r);
R6_v=reshape(R6,1,n2_cell_z*n2_cell_r);
Z6_v=reshape(Z6,1,n2_cell_z*n2_cell_r);
R=[R1_v,R2_v];
Z=[Z1_v,Z2_v];
%target points of DSV
center_DSV_R=0.0; %The center of DSV
center_DSV_Z=0.0;
r_DSV=0.25; %The radii of DSV
%l_DSV=0.045; %Te length of DSV
num1_targets=19;
the = linspace(0,pi/2,num1_targets);
%R1_targets=r_DSV*ones(1,num1_targets);
R1_targets=r_DSV.*sin(the);
%Z1_targets=0:pi/2/(num1_targets-1):pi/2;
%Z1_targets=r_DSV.*cos(Z1_targets);
Z1_targets=r_DSV.*cos(the);
%target points of stray field
R_stray=2.5;
Z_stray=4;
num2_targets=11;
the2 = linspace(0,pi/2,num2_targets);
Z2_targets = Z_stray.*cos(the2);
R2_targets=R_stray.*sin(the2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
A1=ones(num1_targets,n1_cell_r*n1_cell_z);
A11=ones(num1_targets,n1_cell_r*n1_cell_z);
A2=ones(num1_targets,n2_cell_r*n2_cell_z);
A22=ones(num1_targets,n2_cell_r*n2_cell_z);
for n=1:1:n1_cell_r*n1_cell_z
for m=1:1:num1_targets
A1(m,n)=A_mn(R1_v(1,n),Z1_v(1,n),R1_targets(1,m),Z1_targets(1,m));
A11(m,n)=A_mn(R5_v(1,n),Z5_v(1,n),R1_targets(1,m),Z1_targets(1,m));
end
end
for n=1:1:n2_cell_r*n2_cell_z
for m=1:1:num1_targets
A2(m,n)=A_mn(R2_v(1,n),Z2_v(1,n),R1_targets(1,m),Z1_targets(1,m));
A22(m,n)=A_mn(R6_v(1,n),Z6_v(1,n),R1_targets(1,m),Z1_targets(1,m));
end
end
A=[A1+A11,-A2-A22];
%%%%%%%%%
Cz1=ones(num2_targets,n1_cell_r*n1_cell_z);
Cz11=ones(num2_targets,n1_cell_r*n1_cell_z);
Cz2=ones(num2_targets,n2_cell_r*n2_cell_z);
Cz22=ones(num2_targets,n2_cell_r*n2_cell_z);
for n=1:1:n1_cell_r*n1_cell_z
for m=1:1:num2_targets
Cz1(m,n)=A_mn(R1_v(1,n),Z1_v(1,n),R2_targets(1,m),Z2_targets(1,m));
Cz11(m,n)=A_mn(R5_v(1,n),Z5_v(1,n),R2_targets(1,m),Z2_targets(1,m));
end
end
for n=1:1:n2_cell_r*n2_cell_z
for m=1:1:num2_targets
Cz2(m,n)=A_mn(R2_v(1,n),Z2_v(1,n),R2_targets(1,m),Z2_targets(1,m));
Cz22(m,n)=A_mn(R6_v(1,n),Z6_v(1,n),R2_targets(1,m),Z2_targets(1,m));
end
end
Cz=[Cz1+Cz11,-Cz2-Cz22];
%%%%%%%%%%%%%%
Cr1=ones(num2_targets,n1_cell_r*n1_cell_z);
Cr11=ones(num2_targets,n1_cell_r*n1_cell_z);
Cr2=ones(num2_targets,n2_cell_r*n2_cell_z);
Cr22=ones(num2_targets,n2_cell_r*n2_cell_z);
for n=1:1:n1_cell_r*n1_cell_z
for m=1:1:num2_targets
Cr1(m,n)=Ar_mn(R1_v(1,n),Z1_v(1,n),R2_targets(1,m),Z2_targets(1,m));
Cr11(m,n)=Ar_mn(R5_v(1,n),Z5_v(1,n),R2_targets(1,m),Z2_targets(1,m));
end
end
for n=1:1:n2_cell_r*n2_cell_z
for m=1:1:num2_targets
Cr2(m,n)=Ar_mn(R2_v(1,n),Z2_v(1,n),R2_targets(1,m),Z2_targets(1,m));
Cr22(m,n)=Ar_mn(R6_v(1,n),Z6_v(1,n),R2_targets(1,m),Z2_targets(1,m));
end
end
Cr=[Cr1+Cr11,-Cr2-Cr22];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xm=zeros(n1_cell_z*n1_cell_r+n2_cell_z*n2_cell_r,1);
%xm(1:1:n1_cell_z*n1_cell_r)=J1*0.999*area1_cell;
%f=[J1*R1_v,J2*R2_v,J3*R3_v,J4*R4_v]'./J1*area1_cell;
f=[R1_v,R2_v]'; %考虑到电流密度不同的话电流对线包体积贡献的不同
Az=[A;-A;Cz;-Cz;Cr;-Cr];
B=[B0*(1+e)*ones(num1_targets,1);-B0*(1-e)*ones(num1_targets,1);Bz*ones(num2_targets,1);Bz*ones(num2_targets,1);Br*ones(num2_targets,1);Br*ones(num2_targets,1)];
Aeq=[];
Beq=[];
xM=[J1.*area1_cell.*ones(n1_cell_z*n1_cell_r,1);J2.*area2_cell.*ones(n2_cell_z*n2_cell_r,1)];
ff=optimset;
ff.MaxIter=500;
ff.TolX=1e-10;ff.TolFun=1e-10;TolCon=1e-10;ff.Display='iter';
[x,f_opt,key,c]=linprog(f,Az,B,Aeq,Beq,xm,xM,[],ff);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
price=f'*x*2*pi/J1
etime(clock,t0);
curr1 = x(1:10000);
curr2 = x(10001:20000);
curr = curr1 - curr2;
curr = reshape(curr,n1_cell_z,n1_cell_r);
%curr = reshape(curr,n1_cell_z,n1_cell_r);
zmin_plt = 0;
zmax_plt = 1.32/2;
rmax_plt = 1;
rmin_plt = 0.5;
zb = linspace(zmin_plt,zmax_plt,n1_cell_z);
rb = linspace(rmin_plt,rmax_plt,n1_cell_r);
figure('color','w');
surfc(rb(1:end),zb(1:end),curr);
colormap jet
没有合适的资源?快使用搜索试试~ 我知道了~
基于MATLAB实现的电磁优化计算功能,进行线型规划优化电磁设计+使用说明文档.rar
![preview](https://csdnimg.cn/release/downloadcmsfe/public/img/white-bg.ca8570fa.png)
共2个文件
md:1个
m:1个
![preview-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/scale.ab9e0183.png)
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 171 浏览量
2024-05-22
17:32:22
上传
评论
收藏 6KB RAR 举报
温馨提示
CSDN IT狂飙上传的代码均可运行,功能ok的情况下才上传的,直接替换数据即可使用,小白也能轻松上手 【资源说明】 基于MATLAB实现的电磁优化计算功能,进行线型规划优化电磁设计+使用说明文档.rar 1、代码压缩包内容 主函数:main.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 Matlab 2020b;若运行有误,根据提示GPT修改;若不会,私信博主(问题描述要详细); 3、运行操作步骤 步骤一:将所有文件放到Matlab的当前文件夹中; 步骤二:双击打开main.m文件; 步骤三:点击运行,等程序运行完得到结果; 4、仿真咨询 如需其他服务,可后台私信博主; 4.1 期刊或参考文献复现 4.2 Matlab程序定制 4.3 科研合作 功率谱估计: 故障诊断分析: 雷达通信:雷达LFM、MIMO、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计:SOC估计 目标定位:WSN定位、滤波跟踪、目标定位 生物电信号:肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统:DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、DTMF、信号检测识别融合、LEACH协议、信号检测、水声通信 5、欢迎下载,沟通交流,互相学习,共同进步!
资源推荐
资源详情
资源评论
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
收起资源包目录
![package](https://csdnimg.cn/release/downloadcmsfe/public/img/package.f3fc750b.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
![file-type](https://csdnimg.cn/release/download/static_files/pc/images/minetype/UNKNOWN.png)
共 2 条
- 1
资源评论
![avatar-default](https://csdnimg.cn/release/downloadcmsfe/public/img/lazyLogo2.1882d7f4.png)
![avatar](https://profile-avatar.csdnimg.cn/default.jpg!1)
IT狂飙
- 粉丝: 4778
- 资源: 2640
![benefits](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-1.c8e153b4.png)
下载权益
![privilege](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-2.ec46750a.png)
C知道特权
![article](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-3.fc5e5fb6.png)
VIP文章
![course-privilege](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-4.320a6894.png)
课程特权
![rights](https://csdnimg.cn/release/downloadcmsfe/public/img/vip-rights-icon.fe0226a8.png)
开通VIP
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助
![voice](https://csdnimg.cn/release/downloadcmsfe/public/img/voice.245cc511.png)
![center-task](https://csdnimg.cn/release/downloadcmsfe/public/img/center-task.c2eda91a.png)
最新资源
- 高分项目,PID-电机类- 基于红外遥控的多功能智能电风扇控制器源码
- 高分项目,PID-电机类- 基于NRF24L01的PID电机转速控制器,源码+无线PID电机转速控制器+接收机
- 基于Matlab读取sif格式文件,并且画图
- Python基于Django图书管理系统(源码+数据库脚本+截图)高分期末大作业
- Screenshot_20240619_075859.jpg
- 210501230212
- JavaWeb停车收费管理系统的设计与实现
- matlab画图Some-matlab-code-master.zip
- python,python-sorts.rar
- python,python-searches.rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)
安全验证
文档复制为VIP权益,开通VIP直接复制
![dialog-icon](https://csdnimg.cn/release/downloadcmsfe/public/img/green-success.6a4acb44.png)