| 查看: 837 | 回复: 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]) 求高手指点 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有255人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用matlab实现一个沿x轴运动的且半径随x坐标变化的圆的轨迹图
已经有4人回复
- 如何在matlab中计算对称矩阵的高精度本征函数？
已经有0人回复
- MATLAB程序总是出错?
已经有1人回复
- MATLAB问题 程序怎么会出错？Error: Missing variable or function. 怎么修改？
已经有1人回复
- matlab作出的图都缩到一块，如何调整？
已经有7人回复
- MATLAB程序设计与工程应用
已经有144人回复
- matlab提前设定语句运行时间的问题
已经有4人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 绝热毛细管（制冷）校核计算时壅塞流流量的计算
已经有17人回复
- 关于使用lsqcurvefit进行最小二乘曲线拟合时候出现的误差巨大问题求助！！！
已经有7人回复
- 使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响
已经有17人回复
- 分享matlab的一些代码（入门贴）
已经有19人回复
- 求高手指教用matlab拟合相变曲线，为什么不管怎么调整参数误差都这么大，有图
已经有2人回复
- 一种简单已用的热加工图画法
已经有22人回复
- 用MATLAB重新演绎PSCAD图线的方法
已经有6人回复
- 【求助】怎么使用MATLAB求PID与响应的调整时间关系
已经有2人回复
- 【beefly个人文集】画NBO轨道的最简单、最省钱方法
已经有65人回复
2楼2012-11-10 09:12:24
3楼2012-11-10 09:28:57
4楼2012-11-10 10:56:31
【答案】应助回帖
★ ★ ★ ★ ★ ★
小鱼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
10