作业帮 > 综合 > 作业

求matlab大神做循环

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/10 05:51:04
求matlab大神做循环

clear;
clc;
syms Tijicha a k w b r;
a=1.53;
k=0.20;
b=120;

w=300;

I=w*b/(4*pi);
r=1:10:5000;
Er=(5/36)*I*exp(-a*r)+(5/8)*(1+7*exp(-k*r))*I*k/pi/r*exp(-k*r);


f = @(Er, r)(5/36)*I*exp(-a*r)+(5/8)*(1+7*exp(-k*r))*I*k/pi/r*exp(-k*r)-Er;
Er1 = 10;
r1 = fsolve(@(r1)f(Er1, r1), 10)


Er2 = 1;
r2 = fsolve(@(r2)f(Er2, r2),1)
Tijicha=(4/3)*pi*(r2^3-r1^3)


大体模式如下:
while(w