谁会Matlab，帮我把这个程序转换成C语言，3ks

hxs880

clear all;
close all;
ts=0.001;
J=1/133;
q=25/133;
sys=tf(1,[J,q,0]);
dsys=c2d(sys,ts,'z');
[num,den]=tfdata(dsys,'v');
r_1=0;r_2=0;
u_1=0;u_2=0;
y_1=0;y_2=0;
error_1=0;
ei=0;
for k=1:1:1000
time(k)=k*ts;
s=2;
if s==1   %Step Signal
    rin(k)=1.0;
elseif s==2  %Sine Signal
    rin(k)=1.0*sin(1*2*pi*k*ts);
end
yout(k)=-den(2)*y_1-den(3)*y_2+num(2)*u_1+num(3)*u_2;
error(k)=rin(k)-yout(k);
derror(k)=(error(k)-error_1)/ts;
ap=22;bp=8.0;cp=0.8;
kp(k)=ap+bp*(1-sech(cp*error(k)));
ad=0.5;bd=2.5;cd=6.5;dd=0.30;
kd(k)=ad+bd/(1+cd*exp(dd*error(k)));
ai=1;ci=1;
ki(k)=ai*sech(ci*error(k));
ei=ei+error(k)*ts;
u(k)=kp(k)*error(k)+kd(k)*derror(k)+ki(k)*ei;
%Update Parameters
r_2=r_1;r_1=rin(k);
u_2=u_1;u_1=u(k);
y_2=y_1;y_1=yout(k);
error_1=error(k);
end
figure(1);
plot(time,rin,'b',time,yout,'r');
xlabel('time(s)');ylabel('error');
figure(2);
plot(time,derror,'r');
xlabel('time(s)');ylabel('derror');
M=1;
if M==1
    figure(3);
    plot(error,kp,'r');xlabel('error');ylabel('kp');
    figure(4);
    plot(error,kp,'r');xlabel('error');ylabel('kd');
    ad+bd/(1+cd)
    figure(5);
    plot(error,ki,'r');xlabel('error');ylabel('ki');
elseif M==2
    figure(4);
    plot(time,kp,'r');xlabel('time(s)');ylabel('kp');
    figure(5);
     plot(time,kd,'r');xlabel('time(s)');ylabel('kd');
     figure(6);
     plot(time,ki,'r');xlabel('time(s)');ylabel('ki');
end

hxs880

没有人知道么？

王紫豪

有人知道，也不会帮你做苦力

yikaikai

你这个程序 有些matlab自己的画图函数, 没法全部转换

yikaikai

你这个程序 有些matlab自己的画图函数, 没法全部转换

menb

可是俺不懂matlab

cjaizss

原帖由 yikaikai 于 2006-5-27 21:16 发表
你这个程序 有些matlab自己的画图函数, 没法全部转换

matlab里面应该有有编译工具

agan_1007

遗传算法做PID优化的？

scutan

我想除了你给钱, 几乎没有人愿意免费为你做, 除非他是技术狂人.
我去年暑假给别人将一个matlab程序转换成C代码过, 由于涉及到FFT的变换, 几乎所有的运算全是三维数组的转换. 而且matlab的一行代码, 用C要用很多行, 还是循环来完成.
不过有转换工具, 你可以找找.

Godbach

楼上的也做过Matlab转C啊，呵呵，是啊。一旦涉及到复杂的数学运算，Matlab一行代码，C里面要上百行啊。最常见的就是FFT，用C实现一遍FFT，老锻炼人了。