| 查看: 675 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
小鱼118银虫 (小有名气)
|
[求助]
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]) 求高手指点 |
回复此楼» 猜你喜欢
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有6人回复
-
物理学I论文润色/翻译怎么收费?
已经有54人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有62人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用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人回复
4楼2012-11-10 10:56:31
2楼2012-11-10 09:12:24
小鱼118
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5821.5
- 散金: 51
- 红花: 4
- 帖子: 206
- 在线: 161.5小时
- 虫号: 1134404
- 注册: 2010-10-29
- 性别: MM
- 专业: 光学
3楼2012-11-10 09:28:57
【答案】应助回帖
★ ★ ★ ★ ★ ★
小鱼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
