24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

曾全元

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 54.7
帖子: 92
在线: 36.6小时
虫号: 2639060
注册: 2013-09-06
性别: GG
专业: 计算数学与科学工程计算

[交流] 一个简单的编程，谢谢求解已有15人参与

编程求[100，200]之间第一个能被25整除的整数。

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有4人回复
论文终于录用啦！满足毕业条件了已经有27人回复
所感已经有3人回复
要不要辞职读博？已经有7人回复
不自信的我已经有11人回复
北核录用已经有3人回复
实验室接单子已经有3人回复
磺酰氟产物，毕不了业了！已经有8人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有10人回复
26申博（荧光探针方向，有机合成）已经有4人回复

吾生也有涯，而知也无涯

1楼 2013-12-13 15:11:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

超人不怕

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 57.8
帖子: 56
在线: 50.4小时
虫号: 2872303
注册: 2013-12-15
性别: GG
专业: 无机材料化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

4楼: Originally posted by lixiaofei001 at 2013-12-13 20:20:48
for(i=100;i<=200;i++)
{
  if(i%25==0)
  cout<<"我是第一个被整除的，我叫："<<i<<endl;
  continue;
}

超级菜鸟求助，能不能帮我把这道题转化为lsqcurvefit来解？急！
function KineticsEst1_Diff
clear all
clc

% 动力学数据
t = [0  20  40  60  120  180  300];
CAm = [10  8  6  5  3  2  1];

% 用最小二乘样条拟合法计算微分dCA/dt--使用不经过实验点的B样条插值函数
knots = 3;
K = 3;                % 三次B样条
sp = spap2(knots,K,t,CAm);
pp = fnder(sp);       % 计算B样条函数的导函数
dCAdt = fnval(pp,t) % 计算t处的导函数值
rAm = dCAdt;

% 绘制浓度拟合曲线
ti = linspace(t(1),t(end),200);
CAi = fnval(sp,ti);
plot(t,CAm,'ro',ti,CAi,'b-')
xlabel('t')
ylabel('C_A')
legend('实验值','B样条拟合')

% 非线性拟合
beta0 = [0.0053 1.39];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
   lsqnonlin(@OptObjFunc,beta0,[],[],[],rAm,CAm);
ci = nlparci(beta,residual,jacobian);

% 参数辨识结果
fprintf('Estimated Parameters:\n')
fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))
fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2))
fprintf('  The sum of the squares is: %.1e\n\n',sum(residual.^2))

% 绘制反应速率拟合曲线
figure
plot(t,rAm,'ro',t,Rate(CAm,beta),'b*')
xlabel('t')
ylabel('dC_Adt')
legend('Experiment','Kinetic Model')

% ------------------------------------------------------------------
function f = OptObjFunc(beta,rAm,CAm)
rAc = Rate(CAm,beta);
f = rAc - rAm;

% ------------------------------------------------------------------
function rA = Rate(CA,beta)
rA = -beta(1)*CA.^beta(2); % -rA = -dCA/dt = k*CA^n, 其中k=beta(1), n=beta(2)

赞一下

回复此楼

大菜鸟

19楼2013-12-16 12:49:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 24 个回答

fish.yfyh

铜虫 (小有名气)

应助: 42 (小学生)
金币: 607.9
红花: 4
帖子: 108
在线: 35.9小时
虫号: 2867094
注册: 2013-12-12
专业: 化学生物学与生物有机化学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

================Fortran代码===================
   program test
   implicit none
   integer :: i

   do i = 100,200
      if(mod(i,25) .eq. 0) then
            print*, i
      end if
   end do

   stop
   end program
=========================================
结果：
      100
      125
      150
      175
      200

赞一下

回复此楼

2楼2013-12-13 18:28:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖