| 查看: 685 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
MATLAB程序调整
|
|||
|
function gaosi() i=1; for N=1:0.1:6 [T,Y]=ode45(@last00,[0.001,1000],[0.0000001,0.0001]) JI(i)=max(Y(:,1)) JIm(i)=max(Y(:,2)); NN(i)=N; i=i+1; end plot(NN,JI) xlabel('N'); ylabel('Φ_{int}'); hold on; end function dy=last00(t,y) dy=zeros(2,1); dy(1)=dy(2); dy(2)=y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... exp(-y(1)*0.35^(1/2))+exp(-y(1)*0.45^(1/2))+exp(-y(1)*0.55^(1/2))+exp(-y(1)*0.65^(1/2))+... exp(-y(1)*0.75^(1/2))+exp(-y(1)*0.85^(1/2))+exp(-y(1)*0.95^(1/2)))+2*(exp(-y(1)*0.1^(1/2))+... exp(-y(1)*0.2^(1/2))+exp(-y(1)*0.3^(1/2))+exp(-y(1)*0.4^(1/2))+exp(-y(1)*0.5^(1/2))+... exp(-y(1)*0.6^(1/2))+exp(-y(1)*0.7^(1/2))+exp(-y(1)*0.8^(1/2))+exp(-y(1)*0.9^(1/2)))+exp(-y(1)))-... 0.807*log(1/(N*N))/(0.7846+log(1/(N*N)))*y(2)*(1-(1-exp(2.6*y(1)))/(2.6*y(1)))-0.4117*0.1238/(0.7846+log(1/(N*N)))*y(2)-0.5891*0.1238/(0.7846+log(1/(N*N)))*y(2)*((1-exp(-30*y(1)))/(30*y(1)))-... (1-0.1238/(0.7846+log(1/(N*N))))*y(2)-5.8e-16/(8*0.15)*7.444e15*0.1238/(0.7846+log(1/(N*N)))*30*y(2)*y(2) end 运行结果为 ??? Undefined function or variable 'N'. Error in ==> gaosi>last00 at 18 dy(2)=y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in ==> gaosi at 4 [T,Y]=ode45(@last00,[0.001,1000],[0.0000001,0.0001]) 求高手指点 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有233人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用matlab实现一个沿x轴运动的且半径随x坐标变化的圆的轨迹图
已经有4人回复
- 如何在matlab中计算对称矩阵的高精度本征函数？
已经有0人回复
- MATLAB程序总是出错?
已经有1人回复
- MATLAB问题 程序怎么会出错？Error: Missing variable or function. 怎么修改？
已经有1人回复
- matlab作出的图都缩到一块，如何调整？
已经有7人回复
- MATLAB程序设计与工程应用
已经有144人回复
- matlab提前设定语句运行时间的问题
已经有4人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 绝热毛细管（制冷）校核计算时壅塞流流量的计算
已经有17人回复
- MATLAB程序调整
已经有4人回复
- 关于使用lsqcurvefit进行最小二乘曲线拟合时候出现的误差巨大问题求助！！！
已经有7人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 分享matlab的一些代码（入门贴）
已经有19人回复
- 求高手指教用matlab拟合相变曲线，为什么不管怎么调整参数误差都这么大，有图
已经有2人回复
- 一种简单已用的热加工图画法
已经有22人回复
- 用MATLAB重新演绎PSCAD图线的方法
已经有6人回复
- 【求助】怎么使用MATLAB求PID与响应的调整时间关系
已经有2人回复
- 【beefly个人文集】画NBO轨道的最简单、最省钱方法
已经有65人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★
小鱼118: 金币+5, ★★★★★最佳答案, 非常感谢，我同学成功了！ 2012-11-10 13:56:25
dbb627: 金币+1, 感谢应助 2012-11-12 08:57:21
小鱼118: 金币+5, ★★★★★最佳答案, 非常感谢，我同学成功了！ 2012-11-10 13:56:25
dbb627: 金币+1, 感谢应助 2012-11-12 08:57:21
|
程序如下，,[0.001,1000],[0.0000001,0.0001]这两个参数可能不是很合理，可以自己调整一下。 function gaosi() NN=1:0.1:6; JI=zeros(size(NN)); JIm=zeros(size(NN)); for i=1:length(NN) N=NN(i); [T,Y]=ode45(@last00,[0.01,100],[0.001,0.01],[],N); JI(i)=max(Y(:,1)); JIm(i)=max(Y(:,2)); end plot(NN,JI) xlabel('NN'); ylabel('Φ_{int}'); hold on; function dy=last00(t,y,N) dy=[y(2);y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... exp(-y(1)*0.35^(1/2))+exp(-y(1)*0.45^(1/2))+exp(-y(1)*0.55^(1/2))+exp(-y(1)*0.65^(1/2))+... exp(-y(1)*0.75^(1/2))+exp(-y(1)*0.85^(1/2))+exp(-y(1)*0.95^(1/2)))+2*(exp(-y(1)*0.1^(1/2))+... exp(-y(1)*0.2^(1/2))+exp(-y(1)*0.3^(1/2))+exp(-y(1)*0.4^(1/2))+exp(-y(1)*0.5^(1/2))+... exp(-y(1)*0.6^(1/2))+exp(-y(1)*0.7^(1/2))+exp(-y(1)*0.8^(1/2))+exp(-y(1)*0.9^(1/2)))+exp(-y(1)))-... 0.807*log(1/(N*N))/(0.7846+log(1/(N*N)))*y(2)*(1-(1-exp(2.6*y(1)))/(2.6*y(1)))-0.4117*0.1238/(0.7846+log(1/(N*N)))*y(2)-0.5891*0.1238/(0.7846+log(1/(N*N)))*y(2)*((1-exp(-30*y(1)))/(30*y(1)))-... (1-0.1238/(0.7846+log(1/(N*N))))*y(2)-5.8e-16/(8*0.15)*7.444e15*0.1238/(0.7846+log(1/(N*N)))*30*y(2)*y(2)]; end end |
5楼2012-11-10 12:41:02
2楼2012-11-10 09:12:24
3楼2012-11-10 09:28:57
4楼2012-11-10 10:56:31
